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HAUPTAUFSÄTZE 


Kondensationserscheinungen in Überschalldüsen. 
Von Kl. Oswatitsch in Göttingen. 
(Aus (dem Kaiser-Wilhelm-Institut für Strömungsforsechung. ) 
Bei Laraldüsen. welche ron Wasserdampf oder feuchter Luft durchströnf 
werden, zeigen sich Kondensationserscheinungen erst bei sehr starken Uber- 
sättigungen. Bei den in Betracht kommenden Kondensationsgeschiwindig- 
keiten erweisen sich die möglichen Kondensationskernzahlen als riel zu ge- 
ring. um merkliche Nebelbildung zu verursachen. Bei starken Übersäftti- 
gungen tritt ein Zusammenbrechen der Kondensation ein, wie es schon ron 
«ınderen Übersättigungsversuchen her bekannt ist. Die theoretischen Re- 
sultate decken sich rollkommen zufriedenstellend mit Druckrerteilungs- 
messungen in Überschalldüsen, welche für Wasserdampf ron anderen Au- 
toren und für feuchte Luft rom Verfasser selbst durchgeführt wurden. 


l. Einleitung. 

Wenn man feuchte Luft oder Wasserdampf durch eine Lavaldüse strömen läßt, so zeigt 
sich erst etwas stromabwärts vom engsten Düsenquerschnitt Nebelbildung, obwohl die 
Sättigungsgrenze im allgemeinen bereits vor dem engsten Düsenquersehnitt erreicht wird. 
Bei feuchter Luft wird die Nebelbildung dabei vielfach erst bei besonders geeigneten Be- 
leuchtungsverhältnissen sichtbar. Man beobachtet aber bei großen mit feuchter Luft betriebenen 
Düsen verhältnismäßig leicht die sogenannten X-Stöße, deren Zusammenhang mit der Kon- 
(lensation zuerst von Wieselsberger entdeckt wurde. Man spricht daher auch vielfach von 
Kondensationsstößen. Sie treten am Ort des Kondensationsbeginnes in der Düse auf. 

Es ist der Zweck der vorliegenden Arbeit, die Frage nach dem Grund dieser Konilen- 
sationsverzögerung sowie die Frage nach dem Ablauf der schließlich eintretenden Kondensation 
zu beantworten. Zu diesem Zweck wurden Druck- und Dichtemessungen in einer von feuchter 
Luft durchströmten Lavaldüse gemacht, über die hier unter anderem berichtet werden soll. 
Für das Verhalten von Wasserdampf in Überschalldüsen werden Druckmessungen von 
Yellott [1]') und von Binnie und Woods [2] herangezogen. Die Beobachtungen von Nebel- 
bildung in Düsen werden hier nicht weiter berücksichtigt. In den Zusammenhang von Sicht- 
barkeit und Wassergehalt eines Nebels spielen zu viele Parameter herein, als daß es möglich 
wäre, auf Grund von Nebelbeobachtungen verläßliche Schlüsse ziehen zu können. Der Haupt- 
teil der Arbeit ist der theoretischen Behandlung des Problems gewidmet. 


D. Küchemannund F.Vandrey: Über den Eintluß 
der Düse anf Widerstandsmessungen im Freistrahl II 


©,Weber: Halbebene mit periodisch gewelltem Rand 


NH. Behrbohm: Kurze Bemerkung zur graphischen 
Lösung gewöhnlicher linearer Differentialgleiehun- 
gen 1. Ordnung . 


1) Die Zahlen in eckigen Klammern verweisen auf das Sehrifttum-Verzeiehnis am Schluß der Arbeit. 
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Historisch ist zu bemerken, daß als Erster R. v. Helmholtz, der Sohn des berühmten 
Naturforschers, auf die Unterkühlung von Wasserdampf, der durch eine Mündung in die Luft 
austritt, hingewiesen hat. Bei Lokomotiven läßt sich sehr leicht beobachten, daß der aus- 
tretende Wasserdampf unmittelbar hinter der Mündung vollkommen klar bleibt. Die Kon- 
densation so eines Wasserdampfstrahles steht aber mit der Kondensation von Wasserdampf 
in einer Lavaldüse in keinem Zusammenhang. Das einzig Gemeinschaftliche an beiden Er- 
scheinungen ist lediglich der Umstand, daß es sich stets um unterkühlten Wasserdampf 
handelt. Genaue Beobachtungen eines Freistrahles zeigen nämlich, daß die Kondensation in 
der turbulenten Vermischungszone auftritt, während die Kondensationserscheinungen in laval- 
düsen mit Turbulenz gar nichts zu tun haben. Eingehend hat sich Stodola [3] mit der Frage 
nach der Nebelbildung in strömendem Wasserdampf befaßt. Die ersten Beobachtungen über 
Kondensationserscheinungen in mit feuchter Luft durchströmten Lavaldüsen wurden im hiesigen 
Institut im Zusammenhang mit Schlierenbeobachtungen gemacht. Zur Veröffentlichung dieser 
Beobachtungen kam es aber nicht. In den letzten Jahren erschien außer den beiden früher 
zitierten Arbeiten [1] und [2] noch eine Arbeit von Yellott und Holland [#4] über Konden- 
sationserscheinungen in von Wasserdampf durehströmten Überschalldüsen. 


ll. Besprechung experimenteller Arbeiten. 


Sowohl Yellot als auch Binnie und Woods benützen Lavaldüsen von rechteckigem 
(Querschnitt. Die Düsen sind symmetrisch und werden von zwei profilierten Metallwänden 
und zwei Glaswänden begrenzt, welche es gestatten, die Nebelbildung zu beobachten. In der 
Düsenachse befindet sich ein Röhrchen mit einer Anbohrung zur Messung des statischen Druckes. 

Yellott (Bild 2 zeigt Versuch 1 und 3) mißt bei einer ungefähr gleichbleibenden „Kessel- 
temperatur“ von 421° abs. bei drei verschiedenen Drucken von 2040, 1525 und 1410 mm Hg. 
Seine Düse weist ein ganz beträchtliches Druckgefälle auf, welches zweifellos die Meßgenauigkeit 
herabsetzt, aber den Vorteil hat, daß Ablösungserscheinungen im wesentlichen verhindert 
worden sein dürften. Bei diesen Experimenten tritt der Kondensationsbeginn durch eine 
charakteristische plötzliche Verflachung der Druckkurven in Erscheinung. 

Binnie und Woods fBild 3 zeigt Versuch 9) messen bei nahezu gleichbleibendem 
Kesseldruck von 1080 mm Hg. Die Kesseltemperatur schwankt zwischen 435,5" abs. und 
386,7° abs. Dem Druckgefälle ist im Überschallgebiet ein stehender Wellenzug überlagert, 
welcher mit den Kondensationserscheinungen in keinem Zusammenhang steht. Die Verfasser 
führen diese Erscheinung auf den scharfkantigen Düseneinlauf zurück. Bei diesen Versuchen 
äußert sich der Kondensationsbeginn in einem manchmal recht beträchtlichen Druckanstieg. 
Es ist daher fraglich, ob dieser Druckanstieg lediglich einem Kondensationseffekt oder auch 
einer noch hinzukommenden Ablösung der Grenzschicht zuzuschreiben ist. Das große Versuchs- 
material macht diese Versuche besonders wertvoll. 

Ich selbst habe Luft in bekannter Weise durch eine Düse in einen Vakuunikessel ab- 
gesaugt. Die Düse hatte ebenfalls rechteckigen Querschnitt und stellte im wesentlichen die 
Hälfte einer symmetrischen Düse dar. Die Anordnung war so, 
daß in einem rechteckigen Kanal aus Metall ein verschiebbares 
Profil eingebaut war (Bild 1). Am Ort des Profils waren die seit- 
lichen Metallwände durch Glaswände ersetzt. In der dem Profil 
gegenüberliegenden Wand war eine Anbohrung zur Messung des 
statischen Druckes. Durch Verschieben des Profils war es mög- 
lich, den Druckverlauf längs der ganzen Düse zu verfolgen. Die 
jeweilige Stellung des Profils ließ sich mit Hilfe eines Nonius auf 
0,1 mm genau ablesen. Außer dem stationären Druck wurde auch 
die Dichte gemessen. Über die dabei angewandte Methode wird 
an anderer Stelle berichtet. 

Um den Einfluß der Feuchtigkeit auf die Strömung vom Ein- 
fluß der Reibungsschicht, welcher sich vorläufig noch nicht theo- 
retisch erfassen läßt, trennen zu können, wurden zunächst Mes- 
sungen mit trockener Luft vorgenommen. Dies wurde so gemacht, 
daß an die Düse ein weites Rohr angeschlossen wurde, welches 
mit einem Gasometer in Verbindung stand. Preßluftflaschen, 
welche mit trockener Luft von 150 bis 200 at Druck gefüllt waren, 
wurden nun in das Verbindungsrohr abgeblasen und diese Luft 
wurde durch die Düse in den Kessel abgesaugt. Das (rasometer 
diente hierbei nur zum Druckausgleich und nicht zum Auf- 
Bild 1. Einblick in die Düse mit pewahren von trockener Luft, da es sich zeigte, daß es für den 
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Luft wurde mit einem Aßmannschen Aspirations- Psyehrometer sowohl bei den Messungen 
mit troekener Luft wie auch bei den Messungen mit feuchter Luft kontrolliert. Die Messungen 
mit feuchter Luft wurden bei folgenden „Kesseltemperaturen“ #9, und relativen Feuchtig- 
keiten durchgeführt: 


Versuch Yo 
15,0 0,0 
14.0 0,75 
4 19.7 040, 


Da sich die Versuche vor allem durch die relative Feuchtigkeit der Luft unterschieden, wollen 
wir sie nach diesen Werten benennen. Beim Versuch 4,-=0.40 wurde einfach Laboratoriums- 
luft abgesaugt. Die höheren Feuchtigkeitsgehalte wurden dadurch hergestellt, daß vor der 
Düse ein Gerüst aufgestellt wurde, welches so mit Tüchern behängt war, daß die Luft durch 
die Tücher hindurch abgesaugt wurde. Durch verschiedenes Befeuchten der Tücher wurden 
die relativen Feuchtigkeiten von 0,75 und 0,90 erreicht. Diese Werte stellen einen Mittel- 
wert, der während des Durchmessens der Düse abgelesenen Feuchtigkeiten dar. Die Feuch- 
tigkeitsschwankungen betrugen wenige Prozent. (Gewisse Abweichungen der Meßwerte an 
Stellen, wo bereits Kondensation stattfindet. können «daher außer in der Meßungenauigkeit 
auch im gewissen Grad in der Änderung der Versuchsbedingungen ihren Grund haben. Der 
Kondensationsbeginn ist bei meinen Versuchen ebenfalls durch eine plötzliche Verflachung 
des Druckgefälles auf einer kurzen Wegstrecke charakterisiert, welches hier nun auch von 
einer Verflachung des Diehtegefälles begleitet ist (Bild 4). 


Ill. Theoretische Behandlung. 

1. Grundgleichungen. Wir rechnen im folgenden mit einer Düse, welche von feuchter 
Luft, in der sich Wasser- oder Eistropfen befinden, durchströmt wird. Es bezieht sich der 
Index 1 stets auf den Anteil trockener Luft, der Index 2 auf den Wasserdampfanteil und 
der Index 3 auf den Wasser- oder Eisgehalt der Luft. Größen ohne Index beziehen sich auf 
das Gemisch von Wasserdampf und trockener Luft, also auf die feuchte Luft. Wir erhalten 
dann die Lösung für Strömungen von Wasserdampf und Wasser allein dadurch, daß wir die 
Durehflußmenge @, des Anteiles trockener Luft gleich Null setzen. In diesem Falle sind die 
Größen ohne Index einfach gleich den Größen mit dem Index 2. 

Um das Problem nicht allzu stark zu komplizieren. wollen wir folgende vereinfachende 
Annahmen machen. Wir begnügen uns mit einer Stromfadentheorie für stationäre Strömung. 
Im allgemeinen verwenden wir dabei das vereinfachte Bild, daß wir uns die gesamte Wasser- 
ınenge, welche sich an irgendeiner Stelle befindet, in einen oder mehrere Stromfäden von 
Gesamtquersehnitt f, konzentriert denken. Der Querschnitt des Gasstromfadens sei f. Ledig- 
lich bei Betrachtungen über die Kräfte, welehe Gas und Wasser aufeinander ausüben, und 
über die Kondensationsgeschwindigkeit stellen wir uns das Wasser (oder das Eis) als aus 
zahlreichen Tropfen bestehenden Nebel vor. 

Wir müssen annehmen, daß Translationsgeschwindigkeit w und Temperatur T vom 
rer = und trockener Luft gleich sind. Es ist dann auch die Dichte des Gasgemenges 
gleich der Summe der Dichten der Einzelteile und der Druck des Gasgemenges gleich der 
Summe der Partialdrucke. Wir können also schreiben 


Da ferner feuchte Luft, Wasserdampf und trockene Luft dureh denselben Querschnitt hin- 
durehgehen, ist 


ı=h=f. 
Bezeichnen wir en ‘ das spezifische Volumen und mit @ die Durchtlußmenge, so erhalten 
wir wegen ( —=f-o-w die beiden Gleichungen 
Bezeichnen wir die Summe der Dur hflußmengen von Wasserdampf und Wasser mit @,, = @, 


+ G,, dann ist die gesamte Durchflußmenge 6, + @,,. Durch diese wollen wir die einzelnen G 
dimensionslos machen und erhalten so die Durchflußmengenverhältnisse 


G 


Anstatt der Gl. (2) können wir dann schreiben 
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Die Kontinuitätsbedingungen für unser Problem haben nun folgende Gestalt. 
Erstens muß die durch irgendeinen Düsenquerscehnitt F=f-+f, hindurchgehende Menge 
trockener Luft konstant und gegeben sein, zweitens muß die Summe der Durchflußmengen 
von Wasserdampf und Wasser an irgendeinem Düsenquerschnitt konstant und gegeben sein. 

G,=f:o, | 

G,= 9, Wr 9m, 
Gegeben muß ferner auch die Düsenform sein. Wählen wir die «-Riehtung als Strömungs- 
richtung, so ist also F=f+f, als Funktion von x bekannt. Der oberen Gleichung liegt ferner 
(lie Voraussetzung zugrunde, daß alle Wasserteilchen eine einheitliche lediglich vom Ort ab- 
hängige Geschwindigkeit w, haben. Diese Annahme ist schon deshalb nötig, da sich die 
Tropfengrößen und Tropfenzahlen sonst allein durch Zusammenfließen von Tropfen ändern würden. 

Wir werden sehen, daß es nicht möglich sein wird, anzunehmen, daß die durch einen 
(Juerschnitt hindurehfließende Tropfenzahl konstant ist. Nennen wir die auf die Volumeinheit 
(Gas-Wassergemisch kommende Tropfenzahl #» und die sich in der Volumeinheit des 
(Gases und in der Zeiteinheit bildende Tropfenanzahl J — diese Größe wollen wir auch 
Keimbildungshäufigkeit nennen -, so ist die durch einen Düsenquerschnitt F hindurchgehende 
Tropfenzahl gleich »- w,- F. Die im Volumen f-d. in der Zeiteinheit sich bildende Tropfen- 
zahl muß gleich sein der Vermehrung des durch den Düsenquerschnitt F hindurchgehenden 
Tropfentlusses, also: 


(3). 


order auch 
de m, F d.ır 
Diese Gleichung läßt sich gut übersehen. Zwei Umstände führen zu einer Änderung der 
Tropfenzahl in der Volumeinheit auf der Wegeinheit. Erstens der Umstand, daß sich neue 
Tropfen bilden und zweitens der Umstand, daß sich das Durchflußvolumen des Tropfen ent- 
haltenden Gases ändert. 
An einer Stelle, an welcher wir gleiche Tropfenradien & für alle Tropfen voraussetzen 
können, läßt sieh die Durchflußmenge des Wassers auch in folgender Form schreiben: 
G,=f, 3 0 
Wir wenden uns nun den Bewegungsgleichungen zu. Wir haben eine für die 
feuchte Luft und eine zweite für das Wasser abzuleiten. Auf die feuchte Luft wirken außer 
(len Druckkräften die Widerstandskräfte der Tropfen. Bezeichnen wir also den Widerstand 
eines Tropfens mit IW, so erhalten wir die Gleichung 
d + 
f-o-.d.r -f-dp F-.de.n-W. 
Die Beschleunigung der Luftmasse f-o-d. wird bewirkt durch die Abnahme der Druckkräfte 
fdp, vermindert um den Widerstand der dieser Luftmasse zukommenden F- d.-» Tropfen. 
Die Beschleunigung der Wassermasse wird dureh drei Ursachen bewirkt: 1. wirkt auf 
(las Wasser ein Druckgefälle, 2. übt die Luft auf die Tropfen eine Kraft aus und 3. kondensiert 
Wasserdampf, welcher eine größere Geschwindigkeit haben kann als das Wasser, und über- 
trägt «dabei Bewegungsgröße auf dieses. Die durch die dritte Ursache .. Beschleunigung 


ist gleich der Änderung der Bewegungsgröße in der Zeiteinheit, d. i. die Änderung der 
Wassermenge in der Zeiteinheit mal der Geschwindigkeitsabnahme: 
d 


Aus den beiden letzten Gleichungen erhalten wir leicht 


de idp 
| 
de, dp n 


Prinzipiell wäre zu bemerken, daß «diese beiden Gleichungen unter der Voraussetzung 
d@, 
dr 
aber im allgemeinen sehr klein und eine Beachtung der eben genannten Ungleichung ist da- 
her in der Regel überflüssig. 


1@ 
>0 abgeleitet wurden und exakt auch nur für 2 —0 Gültigkeit haben. Das Glied ist 


| 
<un 
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Hätten wir es mit einem Gas ohne Kondensation zu tun, so würden wir durch Heran- 
ziehen des ersten Hauptsatzes der Wärmelehre genügend Gleichungen für die Lösung der 
einfachen Stromfadentheorie von Lavaldüsen haben. Wir wollen daher auch bei unseren 
Problem das Gleichungssystem zuerst durch den ersten Hauptsatz vervollständigen. Zu diesem 
Zweck müssen wir uns vor allem überlegen, auf welche Arten bei der von uns behandelten 
Strömung Wärme produziert wird. 

Am meisten in das Gewicht fällt natürlich die bei der Wasserkondensation freiwerdende 
Verdampfungswärme r. In der Zeiteinheit wird auf der Strecke d.« die Wärmemenge rd@, 
auf diese Weise frei. Ferner verursacht der Reibungswiderstand der Tropfen Wärmeerzeugung. 
welche gleich der an den Tropfen in der Zeiteinheit geleisteten Arbeit sein muß. Für einen 
Tropfen ist diese Größe gleich W-(w Drittens nehmen wir an, daß die Wasserdampf: 
teilchen, welche auf das Wasser auftreffen, unmittelbar abgebremst werden. Dadurch ent- 
steht eine dem Carnotschen Stoßverlust entsprechende Wärmemenge. In der Zeiteinheit 
hat sie für die Strecke dx den Wert „(we w,)”d @,. Bezeichnen wir mit e die innere 

2 
Energie der Masseneinheit, so erhält der erste Hauptsatz folgende Gestalt: 
de de, dt, ‚de, 


de, 


=r:.dG, + F- den Wie gie 


Bedenken wir, daß wir für die zeitliche Änderung, die irgendeine einem Luftteilchen (oder 


Wasserteilehen) anhaftende Größe erfährt dj, \m stationären Zustand auch [bei Wasser 


d 
m, schreiben können. so erhalten wir: 


‚fde, dr, (de, de, ‚de, 


Nachdem aber die Verdampfungswärme r nach einem Resultat der klassischen Thermodynamik 
nichts anderes ist als die Differenz der Wärmeinhalte des gasförmigen und des flüssigen 
Aggregatzustandes 

Pasta), 


erhalten wir durch Multiplikation der ersten Bewegungsgleichung mit G und Addition der 
beiden Gl. (4) zu Gl. (5) unter Berücksichtigung der Kontjnuitätsbedingungen 


dG, dG, 
dx 0 und de dr 
und der Gl. (1) und (2) das Resultat 


Die Energiegleichung erscheint uns demnach wieder in sehr einfacher Form, die 
dadurch charakterisiert ist, daß sie Ausdrücke, welche Reibungsgrößen enthalten. wie etwa 
den Tropfenwiderstand, nicht aufweist. 

Unser Gleichungssystem ist noch nicht vollständig. Vor allem fehlen uns Angaben über 

die Größe der Kondensationsgeschwindigkeit =>. Wir müssen wissen, wieviel Tropfen sich 
bilden und wie rasch die Tropfen wachsen. Die Geschwindigkeit der Tropfen läßt sich dann 
aus den gemessenen Widerstandswerten von Kugeln [5] leicht errechnen. 

Die Gl. (3), (4) und (5), welche wir Grundgleichungen nennen wollen. bilden aber in- 
sofern ein einheitliches Gleiehungssystem, weil hier nichts über die Zusammensetzung des 
Gases, über das Verhältnis von Wassermenge zu (Gasmenge oder über die Kondensations- 
geschwindigkeit vorausgesetzt ist. Man erhält aus diesem System den Spezialfall für Düsen, 
welche von reinem Wasserdampf und von Wasser durchströmt werden. indem man o, = @, 
setzt. 


2. Kondensationsgeschwindigkeit. Beschäftigen wir uns zunächst mit der Kondensations- 
geschwindigkeit eines einzelnen Tropfens, so sehen wir sofort ein, daß zwei physikalische 


| | 
3 
) 
N 
| 
| 
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Vorgänge vor allem Einfluß auf das Tropfenwachstum haben müssen. Es muß erstens für 
den nötigen Massenstrom in Riehtung auf den Tropfen gesorgt sein und zweitens muß dafür 
gesorgt sein, daß die Kondensationswärme abgeführt wird. Geschieht letzteres z. B. nicht, 
so würde der Tropfen allein schon durch Temperaturerhöhung die Sättigungsgrenze erreichen 
und aufhören zu wachsen. In der Art, wie der Wärmetransport und der Massentransport 
vor sich geht, sind wieder zwei Gebiete zu unterscheiden. Ein Gebiet unmittelbar an der 
Tropfenoberfläche, in dem wir molekulare Vorgänge zu berücksichtigen haben, und ein Gebiet 
"im Gas. welches den Tropfen umgibt. Hier ist mit Wärmeleitungsvorgängen und allenfalls 
Diffusionsvorgängen zu rechnen. 

Diese Betrachtung allein schon zeigt, daß wir bei der Ableitung der Kondensations- 
geschwindigkeit den Fall, daß sich ein Tropfen in übersättigter feuchter Luft befindet, vom 
Fall, daß er sich im übersättigten Wasserdampf befindet, trennen müssen, wenn wir nicht 
allzu verwickelte Gleiehungen haben wollen. Im ersten Fall wird nämlich der Wärme- 
transport in überwiegendem Maß von den Luftmolekülen besorgt und für den Massentransport 
spielt die Diffusion eine maßgebliche Rolle. Iın zweiten Fall ist von Diffusion überhaupt 
keine Rede, der Wärmetransport wird aber von den Dampfmolekülen selbst besorgt. 

Wir wollen uns in dieser Arbeit für die ausführliche Behandlung des leichteren Falles 
entscheiden. Dies ist die Kondensation kleiner Wassertropfen in reinem Wasserdampf. Es 
wird sich später zeigen, daß noch ein zweiter Grund für unseren Entschluß maßgebend ist. 
Wir werden auch sehen, daß schon eine verhältnismäßig rohe Kenntnis des Kondensations- 
vorganges zur Lösung unseres Problems ausreicht, wovon wir im folgenden Gebrauch machen 
wollen. Da der Luftanteil bei den folgenden Überlegungen wegfällt, können wir anstatt der 
Größen mit dem Index 2 einfach die Größen ohne Index schreiben. 

Bei sehr kleinen Tropfen werden vor allem die molekularen Vorgänge an der Oberfläche 
maßgebend sein. Wir wollen uns daher zunächst mit dem vereinfachten Bild beschäftigen, 
daß ein kleiner Wassertropfen von der Temperatur 7: und dem dazugehörigen Sättigungs- 
druck pz sich in Wasserdampf von konstanter Temperatur T und konstantem Druck p befindet. 
Der Wärmetransport geht dann wie bei stark verdünnten Gasen vor sich [6] und [7, 8. 27]. Er 
beträgt pro Flächeneinheit und Zeiteinheit 


sl 
R ist die absolute Gaskenstante und m das Molekulargewicht. Nehmen wir an, daß von allen 
auf die Wasseroberfläche auftreffenden Dampfmolekülen y Teilchen kondensieren und nennen 
wir diese Größe den Kondensationskoeffizienten, so errechnet sich die auf der Wasser- 
oberflächeneinheit in der Zeiteinheit kondensierende Wassermenge zu [7, S. 9] 
JE). 
‘ + RI pP: 
Die kondensierende Wassermenge, multipliziert mit der Verdampfungswärme r, muß nun 
den Wärmestrom geben. Das führt uns auf folgende Gleichung 
17), 
Um nun die Größe ps zu eliminieren, linearisieren wir die Dampfdruckkurve an der Stelle 
des Druckes p und der dazugehörigen Sättigungstenperatur 7, und können dann init Hilfe 
(der Formel von Clausins und Clapevyron schreiben 
= (BE ’ , 
dar dabei die Ableitung des Sättigungsdruckes nach der Temperatur. Aus den letzten 
beiden Gleichungen lassen sich die Drucke entfernen und wir erhalten 


T,—T 1/RT, 
Für y müssen wir nach den jetzigen Messungen einen Wert von etwa 3° annehmen [8] 
und [9, S. 195]. 

Die letzte Gleichung zeigt aber, daß auch bei diesem Wert von y zwischen den beiden 
Temperaturdifferenzen T:z— T und 7,— T nur ein Unterschied von etwa 25°/» entsteht. 

Für die Massenzunahme eines Tropfens muß nun folgende Gleichung bestehen 


d Me: R 


XUM 
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Die an der Tropfenoberfläche freiwerdende Wärmemenge muß gleich der abgeführten Wärme- 


menge sein. Durch Ausdrücken von Mz durch 5 erhalten wir dann leicht 


de 3 p 
sr; 
Bedienen wir uns noch der Abkürzung 


dann können wir die Gleichung für die Kondensationsgeschwindigkeit eines Tropfens in 
Wasserdampf auch schreiben 
dE (Tt-D) dM: 
=—o —— oder =4n#io 
da ©,” m, dr 
Hätten wir uns gar nicht um die molekulare Struktur des Gases gekümmert und ein- 
fach angenommen, daß die Kondensationswärme durch Wärmeleitung abgeführt würde, so 
hätten wir bei konstant gehaltener Übersättigung T: - T folgendes Tropfenwachstumsgesetz 
erhalten 


/2k(T:—T 
ro; 

k ist das Wärmeleitvermögen des Dampfes, 7 die Temperatur in weiter Entfernung vom 
Tropfen. Die höchsten Unterkühlungen, die in Düsen auftreten, liegen bei etwa 30°. Wir 


erhalten bei Einsetzen von Zahlenwerten folgende Abhängigkeit der Tropfengröße £ von der 
Zeit, diese stets in Sekunden gemessen, 


und 


Diesen beiden Gleiehungen können wir entnehmen, daß die aus makroskopischen Über- 
legungen gewonnene Formel (9) das Tropfenwachstum in den ersten Zeitmomenten nach Ent- 
stehen des Tropfens viel zu groß angibt. Hingegen liefert für lange Zeiten die Formel (8) zu 
große Werte. Der Übergang beider Gesetze liegt bei der angenommenen Übersättigung bei 
einer Tropfengröße von etwa £x 10°* cm. Diese Tropfengröße wird in einer Zeit von un- 
gefähr # x 10°? see erreicht. Bei Überschallgeschwindigkeiten dürfen wir also ohne weiteres 
einige em nach dem Entstehen des Tropfens ınit Formel (8) rechnen. 

Nachdem der Wärmeleitungseinfluß im Gebiet des uns interessierenden Tropfenwachstums 
eine bedeutend größere Wachstumsgeschwindigkeit zuläßt, ist nun auch die Berechtigung der 
Annahme erwiesen, für unseren Fall im Tropfen umgebenden Gas die Temperatur als kon- 
stant gleich T anzusetzen. 

Den Einfluß der Oberflächenkrümmung des Tropfens auf die Dampfspannung lassen wir 
in unserer Formel unberücksichtigt. Wir erhalten dadurch für die eben entstehenden Tröpf- 
chen mit Radien von etwa £> 10°” cm wahrscheinlich ein zu großes Tropfenwachstum, doch 
werden wir sehen, daß wir uns diese Ungenauigkeit erlauben können. Außerdem können wir 
damit rechnen, daß sich der stationäre Temperatursprung an der Tropfenoberfläche T:  T so 
rasch einstellt, daß wir auch bei veränderlichen Werten von T in der Düse die Vorgänge als 
quasistationär betrachten können und die Formel (8) daher auch für ortsabhängige T:z — T 
verwenden dürfen, 

Nachdem wir uns ein beiläufiges Maß für die Geschwindigkeit, mit der ein Tropfen vom 

Radius £<<10°* em in Wasserdampf wächst, verschafft haben, müssen wir uns vor allem der 
Frage zuwenden, welche Tropfenzahlen wir in einer Düse annehmen dürfen. Bekanntlich 
sind bei den geringen Übersättigungen, wie sie etwa in der Atmosphäre anzutreffen sind, zur 
Tropfenbildung Staubteilchen, die sogenannten Kondensationskerne, nötig. Wir wollen zeigen, 
daß diese bei den Überlegungen über die Kondensationserscheinungen in Lavaldüsen ohne 
Bedeutung sind. 
Nehmen wir an, wir ließen ein Gasvolumen mit der starken Unterkühlung von 30° eine 
Wegstrecke von 10 em in einer Lavaldüse durcheilen. Für diesen Weg stünde dem Teilchen 
eine Zeit von etwa 3:10°* sec zur Verfügung. Die Tropfengröße kann nach Gl. (8a) dabei 
etwa £=3-10-° cm erreichen. Nehmen wir an, es befänden sich die große Zahl von 100000 
Kondensationskernen im em’, so würde in diesem Volumen auf der 10 cm langen Strecke 
eine Wassermenge von etwa 


| 
In, 
®o,n=10"* 
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kondensieren. Das heißt bei mittleren Dampfdichten von o= 1W°’U.G.S.E. würde auf einer 
Streeke von 10 em bei 30° Unterkühlung nur */ı0»o des Wasserdampfes an Kondensations- 
kernen ausfallen! Wir können es uns wohl ersparen, zu zeigen, daß die dabei freiwerdende 
Kondensationswärme keine Rolle spielen kann. 

Wir haben also bei der Ausdehnung von Wasserdampf (dasselbe läßt sich auch leicht 
für feuchte Luft zeigen) trotz Anwesenheit von Kondensationskernen dieselben Verhältnisse 
vor uns, wie bei der adiabatischen Ausdehnung kondensationskernfreien Dampfes in einer 
Nebelkammer. Hiermit ist die Frage nach «dem Grund der starken Unterkühlungen in Über- 
schalldüsen beantwortet und wir müssen uns nun die Frage vorlegen, wieso es schließlich 
trotzdem im Überschallteil der Düse zur Nebelbildung kommt. 

Diese Frage führt uns auf eine andere bereits bekannte Erscheinung. Belindet sich ein 
zesättigtes von Kondensationskernen freies Gas in einer Wilson-Kammer und übersättigen wir 
dieses dureh rasche Expansion, so tritt keine Kondensation ein, wenn die Übersättigung nicht 
allzu groß ist. Expandieren wir aber so stark, daß der Gasdruck etwa «den vierfachen Betrag 
(des Sättigungsdruckes erreicht, so brieht das übersättigte System in jedem Fall zusammen. 
Auch bei Dampfdüsen finden wir den Kondensationsbeginn bei etwa vierfacher Übersättigung 
und wir müssen demnach vermuten, daß wir die Erscheinung des Zusammenbrechens (des 
übersättigten Systems, welehe wir in der Nebelkammer in ihrer zeitlichen Aufeinanderfolge 
nur schwer erfassen können, bei der stationären Strömung von Wasserdampf durch eine Laval- 
«lüse in örtlicher Aufeinanderfolge vor uns abgerollt sehen. 

Der Schluß, daß Kondensation nur an Staubteilchen stattfinden kann, beruht auf der 
Anwendung der Formeln von Thomson und Gibbs für den Dampfdruck über gekrümmten 
Oberflächen 

m 290 
Ps 
(ps ist dabei der zu 7 gehörige Sättigungsdruck und # ist die Oberflächenspannung des 
Wassers) auf beliebig kleine Tropfen. Die Formel muß dort ihre Gültigkeit verlieren, wo es 
sich um molekulare Vorgänge handelt und wir das Gebiet der Kontinuitätsphysik verlassen. 
In einem Gas kommen nämlich dauernd Moleküle in das Gebiet gegenseitiger Molekularkräfte, 
wir sagen, sie stoßen zusammen, und es kommt auch immer wieder vor, daß drei und mehr 
Moleküle sich gleichzeitig im Bereich gegenseitiger Kräfte befinden. Wir können das auch 
so ausdrücken, dal wir sagen, es bilden sich dauernd Tropfen von der Masse mehrerer Mole- 
küle, im allgemeinen zerfallen diese Tropfen aber wieder. Die Zunahme eines Tropfens um 
ein Molekül in der Zeiteinheit ist dabei im Mittel dem Produkt aus seiner Oberfläche und 
dem Gasdruck proportional zu setzen, die Abnahme eines Tropfens um ein Molekül in der 
Zeiteinheit ist dabei wieder proportional seiner Oberfläche zu setzen und außerdem dem 
Dampfdruck des Tropfens ps. Daraus ergibt sich, daß ein Tropfen um so eher erhalten bleibt, 
je mehr er durch zufälliges Auftreffen von Molekülen bereits angewachsen ist und wenn 
schließlich ein Tropfen so groß geworden ist, daß sein Dampfdruck ps: kleiner als der Gas- 
druck ist, so wird er im allgemeinen nieht mehr zerfallen. Diese Tropfenbildung in einen 
übersättigten System wurde von Beeker und Döring [10] berechnet. Die Verfasser nennen 
Tropfen, welche so groß sind, daß sie nicht mehr zerfallen, Kondensationskeime. Die 
Bildungszahl solcher stabiler Tropfen in der Zeit- und Raumeinheit, «die sogenannte Keim- 
bildungshäufigkeit ./, berechnen sie zu 


4 a? 
Daher ist: 
. 
RTo, In? ! 7 
Ps 


\, ist die Zahl der Dampfmoleküle in der Volumeinheit Wasser, Z ist die Zahl der gas- 
kinetischen Zusammenstöße in der Raum- und Zeiteinheit, also die Bildungszahl von „Tropfen“ 
in der Raum- und Zeiteinheit, welche aus zwei Molekülen bestehen. £, ist der Radius eines 
Dampfmoleküls und © seine mittlere Geschwindigkeit. N, ist die Anzahl der Moleküle in 
einem Tropfen, der so groß ist, daß er bei der vorgegebenen Übersättigung mit dem Wasser- 
ddampf gerade im Gleichgewicht ist. 

Der charakteristische Teil dieser Formel ist die e-Potenz. Die Abhängigkeit von Z, 
« und N, von den Übersättigungsverhältnissen ist gegenüber der Abhängigkeit der e-Potenz 
von diesen Größen ganz bedeutungslos.. Wenn wir etwa «die Änderung der verschiedenen 
Größen bei Unterkühlung dureh adiabatische Ausdehnung untersuchen. so sehen wir, daß die 
Potenz bei einer Temperaturerniedrigung von 2" in dem uns bei Wasserdampf interessieren- 
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den Gebiet auf den fünf- bis zehnfachen Betrag anwächst, während sich die Größen Z. « 
und N, nur um einige Prozent ändern. Es kommt also gar nicht auf eine genaue Bestimmung 
der Größen vor der e-Potenz an. Eine Ungenauigkeit hier von der Größe des Faktors 2 wird 
dadurch kompensiert, daß die Erscheinung theoretisch bei einer um den Bruchteil eines 
Grades geänderten Temperatur beginnt. Gewisse Grenzen sind der Genauigkeit vor allen 
dadurch gesetzt, daß die Oberflächenspannung » in der dritten Potenz im Exponenten von v 
vorkommt. Eine genaue Berechnung verlangt also vor allem auch eine genaue Kenntnis der 
Größe o. Dies ist der gewichtig ‚ste Grund, weshalb wir zum Vergleich für die Theorie das 
Beispiel von Wasserdampfdüsen "heranziehen. Beim Absaugen von feuchter Luft durch Über- 
schalldüsen bricht die Übersättigung bei so tiefen Temperaturen zusammen, daß wir es, falls 
unterkühlte Wassertropfen entstehen, mit für diese tiefen Temperaturen unbekannten Ober- 
Hächenspannungen zu tun haben. Nehmen wir aber an, es entstünden Eisteilchen, so hätten 
wir für o die mittlere Oberflächenenergie der Eiskristalle zu setzen. über die man noch 
weniger weiß. Der Spielraum, den uns die Extrapolation aus den bekannten Werten der 
Oberflächenspannung des Wassers bei höheren Temperaturen gewährt, genügt vollauf, um den 
theoretischen Kondensationsbeginn mit «dem praktischen Kondensationsbeginn «durch ent- 
sprechende Wahl des o zur Deekung zu bringen. 

Die von Becker und Döring erree hnete Keimbildungshäufigkeit ist allerdings unter 
ler Voraussetzung stationärer Übersättigungsverhältnisse gewonnen. Wi ir werden diese Formel 
dann für veränderliche Verhältnisse anwenden können. wenn die Entstehungszeit eines sta- 
bilen Tropfens so klein ist, daß sich die Übersättigungsverhältnisse in dieser Zeitspanne nur 
unbedeutend ändern. Bei den Übersättigungsverhältnissen in Dampfdüsen ist die Molekül- 
anzahl eines Tropfens, welcher mit dem Dampf im Gleichgewicht ist, N, ungefähr 30. Das 
heißt, sein Radius ist rund dreimal so groß als der Molekülradius, also etwa gleich 5 - 10 ° em. 
Die Entstehungszeit so eines Tropfens beträgt nach Formel (Sb) etwa 5-10’ sec, d. h. er 
entsteht auf einer Wegstrecke von rd. 2,5: 10° em. Wenn die Genauigkeit der Formel (sb) 
für so kleine Tropfen auch durchaus anzweifelbar ist, so können wir doch annehmen, daß 
sich die Zustandsgrößen während der Entstehungszeit eines Tropfens in der Düse so wenig 
ändern, daß wir auch die Keimbildung als quasistationär behandeln können und daher die 
Formel (11) verwenden dürfen. Wegen der kurzen Entstehungszeit stabiler Tropfen wollen 
wir im folgenden auch keinen Unterschied zwischen Tropfen vom Radius 0 oder vom Radius 
eines Tropfens, welcher mit dem Dampf im Gleichgewicht ist, machen. Die Tropfenradien 
nehmen nach Formel (8) für << 10" * em proportional mit der Zeit, d. h. im wesentlichen mit 
dem Abstand von ihrem Entstehungsort zu. Ihre Masse nimmt also ungefähr mit der dritten 
Potenz des Abstandes vom Tropfenentstehungsort zu, und für die freiwerdende Wärme sind 
daher vor allem die großen Tropfen maßgebend, für welche «die besprochene Vernachlässigung 
ohne Bedeutung ist. 

Nachdem wir die Keimbildungshäufigkeit / als bekannt voraussetzen dürfen. ist es uns 
nun möglich, nach Formel (3a) die Zunalıme der Tropfenzahl in der Raumeinheit pro Weg- 
einheit, die das Volumenelement zurücklegt, anzugeben. Auch für die Kondensations- 
geschwindigkeit jedes einzelnen Tropfens haben wir ein Maß gefunden und es muß uns daher 


möglich sein, die noch fehlende Größe d ° zu bereehnen. Denken wir uns an irgendeine 


Stelle .» der Düse versetzt, an der schon wesentliche Kondensation statifinden soll. Hier 
werden wir verschiedene Tropfengrößen vorfinden entsprechend den verschiedenen Ent- 
stehungsorten .’ der Tropfen. Diejenigen Tropfen, welche an einem Ort entstanden sind, 
wo x’ einen kleineren Wert hat, werden an der Stelle » größer sein als diejenigen, welche 
an einer Stelle .’ entstanden sind, die sich näher an .r befindet. Die Tropfenmasse M: ist 
daher eine Funktion von x und x’. Greifen wir an der Stelle «’ ein Gebiet von der Breite d.’ 
heraus, so entstehen dort in der Zeiteinheit ./- f- d.r’ Tropfen, und da wir einen stationären 
Vorgang betrachten, müssen an der Stelle .r in der Zeiteinheit ebenfalls J:f-d.- Tropfen 
hindurehwandern, welche zwischen .’ und .’ +dzx’ entstanden sind. Der Beitrag. der von 
«diesem Ort herstammenden Tropfen zur Durchflußmenge @, an der Stelle .- ist demnach 
und daher die Durechflußmenge 
x 


Die Zunahme der in der Achsenriehtung ist daher, weil für M= 
und für J—=0 ist. 
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Diese Formel gilt ganz allgemein für alle Überschalldüsen, da bei diesen die Konden- 
sationskerne nie eine Rolle spielen. Erst durch Einsetzen der Kondensationsgeschwindigkeit 
für Wassertropfen in Dampf oder in feuchter Luft wird Gl. (12a) spezialisiert. Wir sehen 


3 


auch, daß für die Richtigkeit der Kondensationsgeschwindigkeit 3 E die Multiplikation der 


Größe d ” mit einem Faktor von der Größenordnung der Einheit recht bedeutungslos ist, da 


das auf dasselbe herauskommt wie eine ungenaue Bestimmung des Faktors vor der e-Potenz 
in der Keimbildungshäufigkeit J. Daher können wir uns auch beträchtliche Ungenauigkeiten 
bei der Ableitung der Kondensationsgeschwindigkeit erlauben. Wesentlich ist es für die 
Größe ro aber, in welche Abhängigkeit man den Tropfenradius von der Bildungszeit / 
stellt, ob man ihn also proportional f oder proportional y f setzt. 

Wir erhalten nun die Kondensationsgeschwindigkeit in Wasserdampf, indem wir die 
Gl. (Sa) in Gl. (12a) einsetzen. Dabei dürfen wir, wie wir gesehen haben, anstatt 7:  T die 
Größe T,— T einsetzen, da sich diese nur um rd. : 25° unterscheiden (Gl. (7)). Wir erhalten 


und 

| 

09 %; | 


Die Keimbildungshäufigkeit ./ ist dabei Gl. (11) und die Größe 4 Gl. (Sa) zu entnehmen. Die 
Frage nach der Kondensationsgeschwindigkeit wäre hiermit beantwortet. 


3. Gleichungssystem für Düsen mit Kondensation. Um einen möglichst handlichen Formel- 
apparat zu haben, wollen wir im folgenden noch einige sehr einleuchtende Vereinfachungen 
machen, die sowohl für Wasserdampf wie für feuchte Luft Gültigkeit haben. 

Wir sahen, daß die Kondensationskerne bei Überschalldüsen keine Rolle spielen können. 
Der Dampf dehnt sich daher zunächst trocken-adiabatisch aus, bis er eine so starke Über- 
sättigung erreicht, daß außerordentlich starke Keimbildung einsetzt, die zum Zusammenbrechen 
des unterkühlten Zustandes führt. Diese starke Abnahme der Unterkühlung unterbindet dann 
jede weitere Keimbildung. Im übrigen Teil der Düse wachsen dann die an der Stelle des 
Zusammenbrechens gebildeten Keime zu Tropfen steigenden Durehmessers an. Erst im Über- 
schallgebiet wird die Keimbildung wesentlich, wie sich leicht abschätzen läßt. Die Düse hat 
also dieselbe Durchtlußmenge wie wenn keine Kondensation stattfinden würde. 

Bei dem späten Kondensationsbeginn in Dampfdüsen sind die erreichten Tropfengrößen 
außerordentlich klein, so daß auch der im vorigen Abschnitt für bestimmte Bedingungen er- 
rechnete Wert von £&=3- 10°° cm als viel zu hoch gegriffen angesehen werden muß. So kleine 
Tropfen machen auch die starke Beschleunigung in Überschalldüsen praktisch vollkommen 
mit [11]. Wir können in unseren Gleichungen also stets setzen: 


wi 


Es kondensiert ferner bis zu hohen Machschen Zahlen stets nur ein Bruchteil des Wasser- 
dampfes. Wir finden z. B.. daß auch unter der Annahme, es fiele stets so viel Wasser aus, 
daß ständig Sättigung herrscht, an der engsten Stelle erst etwa 3"/o des Dampfes kondensiert 
ist. Da die Geschwindigkeiten von Wasser und Dampf gleichzusetzen sind, das Wasser aber 
rund die tausendfache Diehte des Dampfes hat, können wir ohne weiteres f, gegen f vernach- 
lässigen und daher schreiben 


f=F. 
Wir nehmen ferner für den Dampf die ideale Zustandsgleichung als gültig an und bekommen 
laher nach einigen Umformungen, indem wir die Durehtlußmengenverhältnisse (Gl. (2a)) ein- 
führen, folgende Gleichungen 


de I 


— 0 Impulsgleichung. 


Anstatt des ersten Hauptsatzes nehmen wir der Einfachheit halber die Energiegleichung 


1S, 
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Wir wollen nun die Machsche Zahl in bekannter Weise einführen, indem wir schreiben 


Ma’ = 


Die Schallgeschwindigkeit in Nebel bängt allerdings von der Frequenz ab [7]. Wir wollen 


die Machsche Zahl auf die Schallgeschwindigkeit für hohe Frequenzen \ »P- beziehen. 

Ferner bezeichnen wir mit spezifischer Wärme bei konstantem Druck die Größe von e,, bevor 

Kondensation eintritt 


Bei Besprechung der Keimbildungszahl / haben wir gesehen, daß es auf den Ausdruck vor 
der e-Potenz gar nicht so sehr ankommt. Der größte Teil der Keimbildung erstreckt sich auf 
ein kleines Gebiet im Beginn des Zusammenbreehens der Übersättigung, in welchem sich die 
Werte von f, »w und o nur wenig ändern und in dem mengenmäßig noch sehr wenig kondensiert 
ist. Wir können also setzen 
G,+6@,=0r 

wobei die Indizes k auf die charakteristische Stelle des Kondensationsbeginnes hinweisen 
sollen. In der Kondensationsgleichung fällt dann das f; und das f weg. Wir bekommen so 
das Gleiehungssystem 


= Is Kontinuitätsgleichung 
due I dp 
Ma? Impulsgleiehung 
wde 
ide Idpy Ido 1 
9. Kondensationsgleichung 
dır d.r 0% 


und speziell für Wasserdampfdüsen 


d.r 
x (13a). 
| 


Das Gleichungssystem (D gibt eine sehr genaue Beschreibung des Kondensationsvorganges. 
Wir können aber im allgemeinen noch zwei weitere Vernachlässigungen zulassen, indem wir 


die Verdampfungswärme als konstant anschen und 4, < ‚ setzen. Wir machen von dieser 


Vereinfachung Gebrauch, lösen das Gleiehungssystem (I) auf und erhalten 
wd.r  Ma’-i|fd. T d.r | 
Ido Ma Idf. +19: # 
od. Ma? "dx Ma? e„Tdr 


>)» ist das Verhältnis der spezifischen Wärmen. 
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Diese Gleichungen gelten mit ausreichender Genauigkeit für Düsen, welche von feuchter Luft 

oder von Wasserdampf durchströmt werden. Die Düsenform, also die Größe f Se, muß 


natürlich gegeben sein. Zusammen mit den entsprechenden Kondensationsgleichungen, für 
Wasserdampf also Gl. (13a), bilden sie eine vollständige Beschreibung der besprochenen 


Vorgänge. 
Wir wollen den Gl. (11) sofort ein qualitatives Resultat entnehmen. Die Kondensation 
im UÜberschallgebiet verringert das Druckgefälle und kann für | sogar 
c„T 


einen Druckanstieg verursachen, ein Resultat, dessen Bestätigung man sich aus den Versuchen 
leicht holen kann. 


4. Anwendung der Theorie und Vergleich mit Experimenten. Messungen mit trockener 
Luft haben gezeigt, daß man die Verdrängung durch die Grenzschicht auch in sehr kleinen 
Düsen noch vernachlässigen kann und daher die Druckkurven allein schon aus den Düsen- 
abmessungen mit ziemlich großer Genauigkeit bekommt. Anders ist es, wenn durch irgend- 
eine Einlaufstörung der glatte Strömungsverlauf in der Düse gestört ist, wie es bei den 
Versuchen von Binnie und Woods der Fall ist. Hier darf man nicht mehr die Düsen- 
abmessungen für die Durchrechnung einer Druckkurve mit Kondensation verwenden. Man 
muß vielmehr fiktive Düsenabmessungen aus Druckkurven von Versuchen, bei welchen die 
Kondensation erst an einer weiter stromabwärts gelegenen Stelle eintritt, erreehnen und mit 
diesen in die Gl. (ID) eingehen. Man muß also die Düse mit Dampf, bei welchem der 
Kondensationsbeginn möglichst herausgeschoben ist, „kalibrieren*. 

Wir wollen nun unser Gleichungssystem dazu verwenden, den Verlauf des Druckes in 
einer Düse mit Kondensation schrittweise zu bereehnen. Dabei kommt es uns sehr zustatten, 
daß die Strömung bis in das Überschallgebiet hinein genau so verläuft wie bei trocken- 
adiabatischer Ausdehnung. Wir werden die Stelle unseres Rechnungsbeginnes womöglich 
knapp vor dem Kondensationsbeginn auswählen, was sich durch Abschätzen der nötigen 
Kondensationskeimzahlen machen läßt. Zunächst müssen wir unter der Annahme trocken- 
adiabatischer Ausdehnung die Werte der Zustandsgrößen, der Geschwindigkeit, der Unter- 
kühlung usw. errechnen. Es ist gut, dies an einigen Stellen zu machen, bei denen die 
Kondensation noch unbedeutend ist. da man auf diese Weise dann das Entstehen der ersten 


"ropfen und hiermit die Größe 7 nach Gl. (13a) mit Berücksichtigung von (Sa) und (11) 


leicht ermitteln kann. Haben wir diese Größe an der Stelle .r;. so können wir die Zuwächse 
von Druck, Diehte usw. an dieser Stelle mit Hilfe der Gl. (II) gewinnen und so wieder die 
Zustandsgrößen an der Stelle .r; + d.r errechnen. Mit Hilfe der Größen an der Stelle „;,+d.r 


finden wir wieder die Kondensationsgeschwindigkeit d I» und können nun mit unserer Rech- 


nung fortfahren. Man kann das Verfahren dadurch verfeinern, dab man die Größe SR, 
die sich ziemlich rasch ändern kann, an den Zwischenstellen der Schritte berechnet. Hingegen 
dürfte es keinen Vorteil haben, eine kompliziertere Theorie der schrittweisen Berechnung zu 
verwenden, «da man wegen der starken Veränderlichkeit der Funktion J an sehr kleine 
Schritte gebunden ist. Erst an einer Stelle. an der die Keimbildung vollkommen aufgehört 
hat, kann man zu größeren Schritten übergehen. 

Über den Arbeitsaufwand ist zu sagen, daß für einen geübten Rechner etwa 1 bis 2 
Stunden zur Ausrechnung eines Schrittes notwendig sind. Die Schrittgröße betrug bei meinen 
Rechnungen bei dem großen Druckgefälle der Versuche von Yellott '/, mm, bei dem weniger 
starken Druckgefälle der Experimente von Binnie und Woods I mm am Ort der Keim- 
bildung und starker Änderung des Druckgradienten. Nach dem Abklingen der Keimbildung 
wurden an Stelle von wenig veränderlichen Druckgradienten größere Schritte gemacht. 
Bild 2 zeigt Druckkurven von Yellott und die errechneten theoretischen Kurven. Sie 
stimmen innerhalb der Meßgenauigkeit, welche bei diesem Verfasser nicht sehr groß ist, mit 
(len Experimenten überein. Es wurde hier absichtlich keine dimensionslose Auftragung ge- 
wählt. um die Kurven besser trennen zu können. In dimensionsloser Auftragung fallen die 
Druckkurven im Gebiet vor der Kondensation aufeinander, die Versuchspunkte zeigen aber 
beträchtliche Streuung. Der starke Druckanstieg auf der rechten Seite der Zeichnung hat 
mit Kondensation nichts zu tun, sondern rührt von einem Verdichtungsstoß her. Bild 2a 


zeigt die Tropfenzahl, bezogen auf die Volumeinheit des Entstehungsortes 5 0r, sowie die 


Änderung «dieser Größe mit der Fortsehreitungsrichtung für die beiden Versuche. Dieser 
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etwas komplizierte Ausdruck hat gegenüber der Zahl » den Vorteil, nur dort seinen Wert 
zu ändern, wo Keimbildung stattfindet. Ferner zeigt das Bild die Zunahme des Wasser- 
gehaltes g, als Verhältnis von Wasserdurchflußmenge zur gesamten Durchflußmenge. Man 
kann aus diesen Daten eine mittlere Tropfengröße von etwa S=0,5-10°% em berechnen. 
Allerdings kommen sehr verschiedene Tropfengrößen vor. Yellott und Holland [4] finden 
auch, daß die Tropfengröße klein gegen die Wellenlänge des blauen Lichtes sei. Während 
der Kurve von 9, beträchtliche Genauigkeit zukommen dürfte, es könnte sonst auch nieht 
die gute Übereinstimmung mit den Experimenten zustande kommen, dürfen wir an die anderen 
beiden Kurven keine allzu großen quantitativen Ansprüche stellen, weil ihre Richtigkeit mit 
dem Tropfenwachstumsgesetz in Zusammenhang steht. Bild 3 zeigt einen Versuch von 
Binnie und Woods. Die schwarzen Punkte sind Mittelwerte der Druckmessungen von 
Versuchen, bei denen noch keine Kondensation stattfindet. Diese wurden zum „Kalibrieren* 
der Düse verwendet. Die Versuche dieser Verfasser sind sehr genau und die Abweichungen 
der Theorie scheinen hier nicht innerhalb der Meßgenauigkeit zu liegen. 

Trotzdem muß die gute Übereinstimmung der Theorie mit den Experimenten über- 


raschen. Erstens ist bei den Versuchen von Binnie und Woods mit Ablösungserscheinungen 
am Orte des starken Druckanstieges zu rechnen. Dies hat eine Einengung «er Düse zur 


Bild 4 (rechts). 


Bild 3 (unten). 
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Folge. Es ist also wahrscheinlich, daß die von uns verwendeten Werte des Düsenquerschnitts 
nicht ganz richtig sind. Zweitens ist es von vornherein gar nicht sicher, daß die an ruhiger 
Wasseroberfläche gemessene Oberflächenspannung mit derjenigen an so rasch entstehenden 
Tropfen übereinstimmt und wir hätten uns nicht wundern dürfen, wenn die theoretischen 
Kurven systematische Abweichungen nach einer Seite gezeigt hätten, welche bereits «dureh 
einen geringen Fehler in » verursacht hätten sein können. 


Außer den hier gezeigten drei theoretischen Kurven wurden bisher keine Experimente 
nachgerechnet. Die drei Versuche sind aber voneinander wesentlich verschieden und es ist 
daher nicht zu erwarten, daß noch Diskrepanzen von Theorie und Experiment auftreten. 
Bild 4 zeigt meine eigenen Versuche mit trockener und fenehter Luft. Auch hier wurde 
ein theoretischer Druck- und Dichteverlauf errechnet, indem ich ein Wachsen des Tropfen- 
radius proportional zur Entstehungszeit und zur Dampfdiehte annahm. Dieses Gesetz dürfte 
beim Entstehen der Tropfen erfüllt sein, da die Sättigungsdiehte gegenüber der Dampfdichte 
an diesen Stellen tiefer Temperatur bei feuchter Luft zu vernachlässigen ist. Es wurde das 
Entstehen von Eiskristallen vorausgesetzt und die mittlere Obertlächenenergie dieser Kristalle 
so gewählt, daß der Kondensationsbeginn ungefähr richtig herauskam. Die mittlere Ober- 
tHächenenergie liegt dabei durchaus in dem zu erwartenden Bereich. Von einem weiteren 
Ausbau der Theorie wurde aber hier vorläufig Abstand genommen, da man bei einigen 
Arbeitsaufwand durch entsprechende Wahl des unbekannten Tropfenwachstumsgesetzes und 
der unbekannten mittleren Oberflächenenergie von Eiskristallen ohne weiteres Überein- 
stimmung von Theorie und Experiment erzielen kann. Für die Prüfung der Theorie sind 
diese Versuche also wenig geeignet. Sie sind nach der Prüfung der Theorie durch Experimente 
an Dampfdüsen eher dadurch wertvoll, daß man durch sie vielleicht das Wachstumsgesetz 
kleiner Tropfen und deren mittlere Oberflächenenergie finden könnte. Allerdings wäre dann 
zuerst wichtig festzustellen, ob man es mit Wasser- oder Eisbildung zu tun hat. Ferner 
kann mit Hilfe dieser Experimente die Gesetzmäßigkeit des Kondensationsbeginnes bei von 
feuchter Luft durchströmten Lavaldüsen gefunden werden, worüber an anderer Stelle berichtet 
wird [11]. Die striehpunktierten Linien geben die theoretischen Werte von - und 2 . die 

0 
man erhält, wenn man Sättigung in jedem Augenblick voraussetzt. Sie wurden für eine 
Kesseltemperatur von 20°C und einer relativen Feuchtigkeit von y, = 1,00 errechnet. Man 
sieht, daß diese Kurven mit den Experimenten gar nichts gemeinsam haben, die Dichtekurve 
liegt sogar auf der falschen Seite. Bis zur engsten Stelle decken sich die Versuchswerte 
der drei Versuche mit feuchter Luft und des Versuches mit troekener Luft innerhalb der 
Meßgenauigkeit vollkommen. 


Herr Professor Prandt! hat mir die Anregung und zahlreiche wertvolle Ratschläge zu 
dieser Arbeit gegeben, die Deutsche Forschungsgemeinschaft ermöglichte es mir, mit der 
Durchführung der Arbeit zu beginnen. Sowohl Herrn Professor Prandtl als auch der 
Deutschen Forschungsgemeinschaft bin ich daher zu aufriehtigstem Dank verpflichtet. 

An dieser Stelle möchte ich auch Herrn Dr. Hellwig danken, der mich unermüdlich 
bei langwierigen Reehnungen unterstützte. 
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Über den Einfluß der Düse auf Widerstandsmessungen 
im Freistrahl Il. 


Von D. Küchemann und F. Vandrey in Göttingen. 

Die Untersuchungen einer früheren Arbeit desselben Titels ‘Z. angew. Math. 
Mech. Bd. 21 (1941), S. 17 bis 31 werden fortgesetzt und ergänzt. unter be- 
sonderer Berücksichtigung des kurzen Widerstandskörpers dargestellt 
durch einen Dipol,. Die durch die Kanalgrenzen herrorgernfenen Zusatz- 
geschwindigkeiten werden berechnet und” diejenige Lage des Körpers be- 
stimmt, in welcher die Zusatzgeschwindigkeit am Orte des Körpers ver- 
schwindet, in welcher sich also die Einflüsse von Düse und Freistrahl gerade 
aufheben. Die Rechnungen sind für ebene und für rotationssymmetrische 
Strömung durchgeführt. An die Behandlung der ebenen Aufgabe werden 
einige allgemeine Betrachtungen über die bei Randıwertaufgaben mit ge- 
mischten Grenzen auftretenden mathematischen Fragen ungeschlossen. die 
qualitatir auch für die räumliche Aufgabe gelten. 


I. Aufgabenstellung. 

In einer ersten Mitteilung’) wurden Halbkörper und Quellsenkenkörper im Freistrahl 
von kreisförmigem Querschnitt untersucht, mit dem Ergebnis, daß die Düse praktisch keinen 
Einfluß mehr auf die Zusatzgeschwindigkeit am Orte des Widerstandskörpers hat, wenn der 
Körper etwa einen Kanalradius oder mehr von der Düse entfernt ist. Wenn das Problem 
nun nochmals aufgegriffen wird, so hat das vor allem zwei Gründe: 

Zunächst sollte die Reihe der Beispiele dahin vervollständigt werden, daß auch kurze 
Körper ‘(dargestellt durch einen Dipol) behandelt werden, da für diese ein deutlicherer Ein- 
Huß der Düse zu erwarten ist. Weiterhin sollte genauer festgestellt werden, wie die Düse 
wirkt, wenn sich der Körper sehr nahe an der Düse befindet. 

Schließlich sollten einige grundsätzliche Fragen an Hand einer Reehnung für die ent- 
sprechende ebene Strömung aufgeklärt werden. 

Im ersten Teil wird die ebene Strömung untersucht, welche sich mit Hilfe der kon- 
formen Abbildung berechnen läßt, und im zweiten Teil wird die rotationssymmetrische 
Strömung mit Hilfe von Singularitätenbelegungen auf den Kanalgrenzen behandelt. In beiden 
Fällen wird dabei noch die Voraussetzung gemacht, daß sich die freie Strahlgrenze nicht de- 
formiert. Wieweit diese Annahme berechtigt ist, müßte eine besondere Untersuchung zeigen. 
Jedoch ist anzunehmen, «daß bei den kleinen Entfernungen zwischen Düse und Körper, wie 
sie hier untersucht werden, der mögliche Einfluß der Strahldeformation nur klein sein wird, 
so daß der Einfluß der Düse. auf den es uns hier ankommt. doch im wesentlichen richtig 
erfaßt wird. 

ll. Die ebene Strömung. 

Für den Fall einer ebenen Strömung läßt sieh nach der Methode der konformen Ab- 
bildung eine verhältnismäßig einfache Lösung der vorliegenden Aufgabe in geschlossener Form 
angeben. Die praktische Bedeutung des ebenen Falles ist zwar nur gering, da seine Voraus- 
setzungen meist nicht erfüllt sind, dafür gestattet aber seine Lösung einen Einblick auch in 
(die mathematischen Fragen, die bei der räumlichen Aufgabe auftreten. 

Wir betrachten die folgende Anordnung 


IR. 
(Bild 1a): Die beiden Halbgeraden y= +5 D, 
der z-Ebene bilden die festen Wände — 
des Kanals, die Halbgeraden „= _ ö 
die Grenzen des Freistrahles. Die Kanalbreite a) 2= Ebene 8) z=Ebene 
ist dann gleich D, die Geschwindigkeit der i 
ungestörten Strömung im Kanal sei Als 
Störung nehmen wir einen Dipol in der Kanal- | | 
achse bei a mit dem Moment 
Parallel c) 2, =fbene d) 23=£bene 
Dieser würde in einer ungesti rten a E . Bild 1. Aufban der Abbildungsfunktion zur Lösung der 
strömung von der Geschwindigkeit einen 


1) D. Kücheman nu. F. Vandrey: Über den Einfluß der Düse (oder des Auffangtricehters) auf Widerstands 
messungen im Freistrahl. Z. angew. Math. Mech. Bd. 21 (1941), S. 17 bis 31. 
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Kreis vom Durchmesser d als Grenzstromlinie liefern: wenn wir d- D voraussetzen, so er- 
setzt der die Störung darstellende Dipol einen Körper, der von einem Kreiszylinder mit dem 
Durchmesser d nur wenig abweicht. 

Zu der Dipolströmung soll nun eine Zusatzströmung bestimmt werden, so daß für die 
durch Überlagerung beider entstehende Strömung die festen Wände Stromlinien und die Frei- 
strahlgrenzen Potentiallinien werden. Die erste ‚dieser Forderungen ist genau, die zweite da- 
gegen eine Näherung. Eigentlich wäre eine neue Freistrahlgrenze so zu bestimmen, daß auf 
ihr die Geschwindigkeit (und damit der Druck) der Gesamtströmung konstant ist. Die 
Näherung ist zulässig, solange «die Geschwindigkeit der Zusatzströmung auf der alten Frei- 
strahlgrenze klein bleibt gegen die Hauptgeschwindigkeit 7. 

Zur Bestimmung der Zusatzströmung bilden wir den Kanal konform auf ein Gebiet ab. 
bei dem sich die Randbedingungen durch Spiegelungen erfüllen lassen. Hierzu setzen wir 
zuerst 


Dem Kanal in der z-Ebene entspricht dann in der z,-Ebene ein Streifen zwischen den Ge- 
raden 9, =0 und y, = (Bild 1b). Diesen Streifen bilden wir weiter so auf die obere Hälfte 
einer 2,Ebene ab, daß die untere feste Wand in die Strecke O<.r,<|1, die obere in die 
Strecke 1<a#,<0, die untere Freistrahlgrenze in die Halbgerade 1<x,<» und die obere 
in die Halbgerade »<.#,<1 übergeht (Bild Te). Dies leistet die Abbildung 


Schließlich bilden wir die obere Hälfte der 2,-Ebene noeh mittels der Funktion 


auf den Halbstreifen = einer 2,Ebene ab (Bild Id) °). Faßt ıman diese Ab- 
bildungen zusammen. so erhält man den Zusammenhang zwischen der =-Ebene ımd der z,-Ebene 
in der Form 

Der Dipol, der in der z-Ebene bei z2= «a lag, liegt in der z,-Ebene bei 


Die Strecke soll Stromlinie, die beiden Halbgeraden +5. m 


sollen Potentiallinien werden 


Diese Randbedingungen lassen sich nun leicht wu” 
«durch Spiegelung erfüllen (Bild 2). Hierzu spiegeln | 
wir zunächst den Dipol fortgesetzt unter jeweiliger | 
Umkehr seines Vorzeichens an den (reraden = + , 
7 
und erhalten so eine Dipolreihe längs der Geraden 
7, = Ar Ziner«!D in den Punkten „,=0, +2n, 
+47, ... (nach unten gerichtete Dipole) und 
+37, +57, .... (nach oben gerichtete 
Dipole). Die ganze Dipolreihe spiegeln wir noch- Bild 2. Lösung der ebenen Aufgabe durch 
Spiegelungen in der 23-Ebene. 
mals unter Umkehr ihrer Richtung an der reellen 
Achse „, 0 und erhalten eine zweite Reihe von Dipolen auf der Geraden y, = - Ar Sin er «ID, 


Für beide Reihen läßt sich das komplexe Potential in geschlossener Form angeben. Man 


hat für die obere Reihe 
. . 2 2 2 20202 


und für die untere 

PB, —=—kisin (ea, +aresinieteD) . . 2 222. ..(7b, 
wo k so zu bestimmen ist, daß sich in der z-Ebene für den Dipol die riehtige Stärke ergibt. 
Insgesamt erhält man in der z-Ebene 


P=k [1/sin (are sin id e:!D — are 
I /sin (are sin are sin 
oder | 
| m. 


2) Statt dessen könnte man aueh die längs der Strecke -I1<xa<-+ 1 aufgeschnittene z293-Ebene auf das Äußere 
des Einheitskreises einer 23-Ebene abbilden nnd die 


tandbedingungen durch Spiegelung an diesem Kreise befriedigen. 


ere 
en. 
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Die Konstante k erhält man dabei durch Entwicklung von ® für z=> a und Vergleich mit (1). 


Es wird 
Aus dem Potential erhalten wir die komplexe Geschwindigkeit 
dd 1 | 
dz D’ 4 


° 


(Yl+e-:7aD +y1+e-?7:1D) 


Das durch (9) gegebene Potential erfüllt die Bedingung, daß die festen Wände des 
Kanals Stromlinien und die Freistrahlgrenzen Potentiallinien sind. Trotzdem kann es noch 
nicht als endgültige Lösung der Aufgabe angesehen werden, denn die Forderung konstanten 
Potentials auf der Freistrahlgrenze ist nur dann eine Näherung für die Forderung konstanten 
Druckes, wenn die Normalgeschwindigkeiten klein sind gegenüber der Geschwindigkeit der 
Parallelströmung. Dies ist aber in der Nähe der Düsenmündung (22/D=+i) auch für 
kleine d/D nicht der Fall, denn eine Entwicklung der Geschwindigkeit (11) nach Potenzen von 


ergibt 


1 

Die Geschwindigkeit geht dort also wie [ * gegen Unendlich. Es gibt nun Aufgaben mit ge- 
mischten Grenzen, bei denen es grundsätzlich unvermeidbar ist, daß die Geschwindigkeit an 
den Punkten unendlich wird, in denen die Randbedingung wechselt’). Ein Beispiel hierfür 
bietet die Frage nach der Korrektur des Anstellwinkels eines Flügels in einem teilweise ge- 
schlossenen Kanal. Physikalisch würde dieses Unendlichwerden der Geschwindigkeit eine 
starke Verformung des Strahlrandes an den Enden der freien Grenzen bedeuten, die in Wirk- 
lichkeit aber gar nicht eintritt. In dem hier vorliegenden Falle läßt sich nun auch theoretisch 
eine endliche Geschwindigkeit auf der ganzen Grenze erreichen, wenn man dem Potential ® 
noch ein weiteres Potential hinzufügt, das wir als Düsenpotential ®) bezeichnen wollen und 
das folgende Eigenschaften besitzt: 


1. D) hat die festen Wände als Stromlinien und die Freistrahlgrenzen als Potential- 


linien. 


3. er ist sonst im Kanal und im Freistrahl einschließlich der Grenzen endlich. 


Dies Potential entspricht einer Strömung, die senkrecht durch die Freistrahlgrenze eintritt 
und dabei die Düsenmündung in entgegengesetzter Richtung wie die vom Potential ® her- 
rührende Strömung umfließt. Weiter innen im Kanal geht die von ®) herrührende Strömung 
dann in eine der Anströmgeschwindigkeit !’ entgegengerichtete Parallelströmung über. Physi- 
kalisch kann man die Düsenströmung dahin deuten, daß die Düse durch den Körper zum 
Teil versperrt wird und deshalb der Druck im Kanalinnern höher und die Geschwindigkeit 
dort geringer wird als ohne den Körper. 

Aus der Existenz des Potentials ®) ergibt sich, daß die gestellte Aufgabe — im ebenen 
wie im räumlichen Falle — mathematisch nicht eindeutig lösbar ist, da das Potential O+4®) 
bei beliebigem 4 die zunächst gestellten Bedingungen erfüllt. Unter allen diesen Lösungen 
gibt es jedoch eine, für welche die Geschwindigkeit auch an der Düsenmündung noch endlich 
bleibt. Wenn wir nun endliche Geschwindigkeit an der Düsenmündung zusätzlich fordern 
und so diese Lösung vor den übrigen auszeichnen, so stützen wir uns hierbei auf das von 
Prandtl?) in die Hydrodynamik eingeführte Axiom, nach dem von mehreren mathematisch 
möglichen Lösungen einer Potentialaufgabe diejenige das beste Bild von den Verhältnissen 
in wirklichen Flüssigkeitsströmungen gibt, bei der unendliche Geschwindigkeiten möglichst 
weitgehend vermieden sind. 

Um nun für ®, einen analytischen Ausdruck zu finden, gehen wir davon aus, daß es 
weit im Kanal einer Parallelströmung entspricht. Wir erhalten daher nach Abbildung in die 
2,-Ebene die Strömung einer Senke im Nullpunkt. Um die Randbedingungen des senkrechten 
Durchströmens durch die Freistrahlgrenzen zu erfüllen, müssen wir diese Senke fortgesetzt 


3) Ein Beweis dieser Tatsache wird im Anhang der Arbeit gegeben. 4 
4) L. Prandtl u. A. Betz: Vier Abhandlungen zur Hydrodynamik und Aerodynamik. Göttingen 1927, S. 15. 
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unter Umkehr ihres Vorzeichens an den Geraden «,—= + 5 spiegeln. Wir erhalten so eine 


Reihe von Senken in den Punkten @,=+2nz und von Quellen in den Punkten 
+2 n+Na(n=0,1,2,...), deren Potential sich leicht in geschlossener Form angeben 


läßt. Es wird j 1 
oder in der z- Ebene 
(15) 
mit 
Zur Bestimmung von A entwickeln wir (16) in der Nähe der Düsenmündung 
d 
um. 
Da das Glied mit [ “: in der gesamten Zusatzgeschwindigkeit fehlen soll, erhalten wir durch 
Vergleich mit (13) 
Insgesamt erhalten wir also für die Zusatzgeschwindigkeit einschließlich der Dipolströmung 
1 
dz D4 
+1 2 (19). 
Wir berechnen nun zunächst die Zusatzgeschwindigkeit weit im Kanal (= —»), die von der 
Versperrung der Düse herrührt. 
x) — (20), 
dz 
Ud?n? 


In Bild 3 ist die dimensionslose Zusatzgeschwindigkeit (u (— »)/ U)(D?/d’) in Abhängigkeit 
von der Lage 2a/D des störenden Dipols aufgetragen. Sie ist erwartungsgemäß dann am 
größten, wenn der Dipol ungefähr in der Düsenmündung liegt. Weiter berechnen wir die 
Zusatzgeschwindigkeit am Orte des Dipols z2=a (ohne den Einfluß des Dipols selbst). 


U 
(23). 
Für a=-- x erhalten wir als Grenzwert die Zusatzgeschwindigkeit am Orte eines Dipols, 
der sich in einem geschlossenen Kanal befindet. 
2 


Für a=+»w erhalten wir entsprechend die Geschwindigkeitskorrektur für einen Dipol im 
Freistrahl 

Zwischen diesen beiden Grenzlagen gibt es eine Stelle a, für welche die Geschwindigkeits- 
korrektur «(a) gleich Null wird, sie ergibt sich nach (23) aus der Gleichung 


zu f 
2a 


Die Stelle verschwindender Korrektur liegt also noch innerbalb der Düse. In Bild 4 ist die 
dimensionslose Zusatzgeschwindigkeit am Orte des Dipols (uw (a)/ U) (D*/d”) in Abhängigkeit 
von der Lage des Dipols 2a/D aufgetragen. An dem Bilde fällt auf, daß der Übergang von 
der Korrektur für den festen Kanal zu der für den Freistrahl nicht monoton erfolgt, sondern 


BER. 
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Bild 3. Ebene Strömung. Zusatzgeschwindigkeit weit — 4 
im !Innern des festen Kanalteiles in Abhängigkeit chen | 
von der Lage des Dipols. | 
-15 -10 -05 [7 05 70 15 20 25 
2z 


Bild 5. Ebene Strömung. Zusatzgeschwindigkeit in der 
Kanalachse, wenn der Dipol in der Düsenmündung liegt (1). 
Zum Vergleich sind die Zusatzgeschwindigkeiten für einen 
20 r being Dipol in einem geschlossenen Kanal (2) und in einem 


Freistrahl (5) eingezeichnet. 
| 
057 1 
| | zZ | 
| | 7 | 
S | SS =) | 
=2 7 7 2 
| | | 
Be io 05 pi 05 10 15 Bild 6. Ebene Strömung. Zusatzgeschwindigkeit in der 
Kanalachse, wenn der Dipol einen Düsenhalbmesser vor 
der Düsenmündung liegt (1). Zum Vergleich ist die 
Bild 4, Ebene Strömung. Zusatzgeschwindigkeit am Zusatzgeschwindigkeit für einen Dipol im Freistrahl (2) 
Orte des Dipols in Abhängigkeit von seiner Lage. eingezeichnet. 


daß der Wert der Freistrahlkorrektur zunächst in der Nähe der Düse beträchtlich über- 
schritten wird und dann der Betrag der Korrektur wieder gegen den für den reinen Frei- 
strahl abnimmt. Dies Verhalten kommt durch den ziemlich starken Einfluß des Düsenpotentials 
in der Nähe der Düse zustande. 

In Bild 5 ist noch die Geschwindigkeitskorrektur längs der Achse (ohne den Einfluß 
des Dipols selbst) für den Fall gezeichnet, daß der Dipol in der Düsenmündung liegt (a=0) 
(Kurve 1). Zum Vergleich sind die entsprechenden Kurven für den festen Kanal (Kurve 2) 
und den Freistrahl (Kurve 3) mit eingetragen. Der Verlauf ist in diesem Falle erheblich 
anders als für den Freistrahl oder den festen Kanal, auch sind die auftretenden Zusatz- 
geschwindigkeiten wesentlich größer als in diesen beiden Fällen. Der Grund hierfür liegt 
wiederum hauptsächlich im Einfluß des Düsenpotentials. Geht man mit dem Körper vor die 
Düse, so klingt dieser Einfluß ziemlich rasch ab. Als Beispiel einer solchen Lage ist in Bild 6 
die Geschwindigkeitskorrektur in der Achse für eine Entfernung des Dipols von der Düse 
gleich dem halben Düsendurchmesser (Kurve 1) mit der Korrektur des reinen Freistrahles 
(Kurve 3) verglichen. 

Bei den vorstehenden Überlegungen ist als Bezugsgeschwindigkeit stets die Geschwin- 
digkeit U der ungestörten Parallelströmung genommen worden. Diese herrscht in der ge- 
störten Strömung auf dem Rande und weit von der Düse auch im Innern des Freistrahles, 
nicht aber weit im Innern des festen Kanalteiles.. Man mißt nun die Strahlgeschwindigkeit 
zumeist durch Überdruck Ap in der Düsenvorkammer des Kanals nach der Beziehung 


2 
wo F, der Düsenquerschnitt, F, der Querschnitt der Düsenvorkammer ist. Im Falle ebener 


Strömung kann man F,/F, durch das Verhältnis der Durchmesser D,/D, ersetzen. Nun hat 
man in der Düse die mittlere Geschwindigkeit 


in der Düsenvorkammer ist daher die Geschwindigkeit aus Kontinuitätsgründen 


en | | | | 
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Der Druck in der Düsenvorkammer ist dann 


Der Druck im Freistrahl ist i 


Aus diesen beiden Gleichungen folgt für die Geschwindigkeit 7’ 


Dies Ergebnis läßt sich so aussprechen, daß wegen der Versperrung der Düse durch den 
Körper der Endquerschnitt «des Freistrahles kleiner wird als der Düsenquerschnitt. Der Be- 
trag dieser Verkleinerung kann aus den in Bild 3 gegebenen Werten für «. (-») leicht be- 
rechnet werden. Da bei der Auswertung von Messungen ein Staudruck zugrunde gelegt wird, 
der einfach nach (28) berechnet ist, so erhält man Beiwerte, die im Verhältnis 
zu klein sind. 


III. Die rotationssymmetrische Strömung. 


Den Untersuchungen dieses Abschnittes wird eine kreisförmige Düse°) mit anschließen- 
dem zylindrischen Freistrahl von Kreisquerschnitt zugrunde gelegt. Wir benutzen ein Ko- 
ordinatensystem mit dem Nullpunkt auf der Symmetrieachse in der Düsenmündung und be- 
zeichnen mit « die axiale, mit r die radiale Richtung. Der Durchmesser des Kanals sei D, 
die ungestörte Anströmgeschwindigkeit sei U’. Als Störung wählen wir einen Dipol auf der 
Achse im Abstand #—=a von der Düsenmündung mit dem Moment 


. 


welcher in einer ungestörten Parallelströmung eine Kugel mit dem Durchmesser d als Grenz- 
stromlinie liefern würde. h 

Die Rechnung ist nach der Methode der in Fußnote 1 angegebenen Arbeit durchgeführt, 
so daß darauf nicht näher eingegangen werden soll. Sie beruht auf dem Singularitäten- 
verfahren, wonach zur Erfüllung der durch den Dipol gestörten Randbedingungen die feste 
Wand mit Quellringen und die freien Strahlgrenzen mit Wirbelringen belegt werden. Wir 
erhalten so eine Integralgleichung für die Dichte « der Belegung, welche numerisch durch 
Iteration gelöst werden kann. Die Konvergenz des Iterationsverfährens machte bei den an- 
gegebenen Beispielen gewisse Schwierigkeiten, so daß teilweise bis zur 15. Näherung gereclmet 
werden mußte. Jedoch ist die Konvergenz als gesichert anzusehen. — Die Rechnung wurde 
durchgeführt für 3 Lagen des Dipols: 2a/D=0; 0,5; 1,0. In Bild 7 sind die Belegungs- 
dichten u aufgetragen. Für negative « (feste Wand) bedeutet «u die Ergiebigkeit von Quell- 
ringen, für positive « (Freistrahl) die Stärke von Wirbelringen. Die Ausgangsnäherung u‘ 
(doppelte Normalgeschwindigkeit des störenden Dipols auf der festen Wand, Tangential- 
geschwindigkeit auf den Freistrahlgrenzen) ist zum Vergleich gestrichelt eingetragen. 

Bereits in der ersten Mitteilung wurde darauf hingewiesen (S. 26), daß unsere numerische 
Methode nicht geeignet ist, die Verhältnisse an der Übergangsstelle genau zu erfassen. Ein 
Unendlichwerden der Geschwindigkeit an dieser Stelle wird bei der Integralgleichungsmethode 
dadurch vermieden, daß die Belegung an der Kanalmündung nach Null geht. Damit erhalten 
wir, im Gegensatz zur ebenen Potentialaufgabe des vorigen Abschnittes, wie in der ersten 
Mitteilung sofort eine eindeutige Lösung, welche im Innern der festen Wände einen Stau 
hervorruft (entsprechend der Gesamtlösung des vorigen Abschnitts, welche das Düsenpotential 
mit enthält). Die aus den Singularitätenbelegungen berechneten Zusatzgeschwindigkeiten 
dürften trotz der Unsicherheit an der Düsenmündung ein gutes Bild von den wirklichen Ver- 
hältnissen geben, besonders in der Kanalachse. 

Die axialen Zusatzgeschwindigkeiten u am Orte des Dipols sind in Bild 8 aufgetragen. 
Die Zusatzgeschwindigkeit für eine feste Lage des Dipols in dessen Nachbarschaft zeigt 
Bild 9. Dabei ist « die Differenz zwischen der wahren durch die Kanalgrenzen veränderten 
Geschwindigkeit an dem betreffenden Punkt und der Geschwindigkeit U auf dem Strahlrand 
vermehrt um die vom störenden Körper herrührende Geschwindigkeit. In Bild 8 ist außer- 


6) Die Form der festen Wand ist aus theoretischen Gründen, welche in der ersten Mitteilung (Fußnote ]) aus- 
führlich dargestellt sind, bis etwa zum 1,5-fachen Kanalradius von der Düsenmündung kreiszylindrisch und erweitert 
sieh dann, um schließlich in eine zur Achse senkrechte Wand überzugehen (Düse 1 der angeführten Arbeit). 
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dem Kanalrand für verschiedene Lagen des Dipols zur Düse. 


Bild 7. Rotationssymmetrische Strömung. Singularitätenbelegung auf x=0 | 
| 


Bild 8 (rechts oben). Rotationssymmetrische Strömung. Zusatzge- 05 
schwindigkeiten im Mittelpunkt für die Kugel in verschiedenen Ab- | | 
ständen von der Düse. | | 


Bild 9 (rechts unten). Rotationssymmetrische Strömung: Oben: Zusatz- | | | ) 
geschwindigkeit in der Achse für verschiedene Lagen des Dipols. n) 27 . 22 


Unten: Zusatzgeschwindigkeit in Abhängigkeit von der Entfernung von u 
der Achse für die Kugel in der Düsenmündung. LEREN VIEIITER 


@ 


dem noch eine Geschwindigkeit u" — u — ü eingetragen, wobei @ den Mittelwert der Axial- 
komponente der Zusatzgeschwindigkeit im zylindrischen Teil des geschlossenen Kanals be- 
deutet®). Hat man also den statischen Druck in der Düsenvorkammer gemessen und daraus 
mittels des Düsenfaktors nach (32) die Geschwindigkeit auf dem Strahlrand bestimmt, so 
gibt w die Größe der Geschwindigkeitskorrektur an. Aus Bild 8 ersieht man nun die be- 
merkenswerte Tatsache, daß für « der Einfluß der Düse auf Körper im Freistrahl äußerst 
gering ist. Noch für solche Körper, welche sich einen Kanalradius von der Düse entfernt be- 
finden, haben wir für « praktisch nur Freistrahleinfluß. Erst bei noch größerer Annäherung 
an die Düse macht sich deren Wirkung bemerkbar. Die Stelle verschwindender Korrektur 
liegt etwa '/ı Kanalradius von der Düsenmündung entfernt. Dieses Ergebnis ist in der Haupt- 
sache dadurch bedingt, daß bei Annäherung an die Düse wohl die Quellbelegungen der Düse 
(Bild 7) stark anwachsen, aber ebenfalls die Wirbelbelegungen des Freistrahls; beide rufen 
entgegengesetzt gerichtete Zusatzgeschwindigkeiten hervor, so daß im Endergebnis der Düsen- 
einfluß wieder herabgesetzt wird. 

Die in Bild 8 eingetragene Zusatzgeschwindigkeit «* würde Bedeutung erlangen, wenn 
experimentell direkt die mittlere Geschwindigkeit in der Düsenvorkammer bestimmt würde. 
In diesem Falle liegt die korrekturfreie Stelle (u*=0) wesentlich weiter, nämlich einen Kanal- 
radius, von der Düsenmündung entfernt. — Aus Bild 8 kann man ebenfalls & = u — u* ent- 
nehmen. 

Bild 9 zeigt vor allem, daß sich die Zusatzgeschwindigkeiten außerhalb der Achse für 
Körper mäßiger Dicke nur wenig von den Werten auf der Kanalachse unterscheiden, so daß 
es stets genügt, die letzteren allein zu berechnen. 


6) u entspricht u (— ©) im ebenen Fall und ist stets der Anströmung entgegen gerichtet (<0). 
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IV. Zusammenfassung. 


Zur Bestimmung der Geschwindigkeitskorrektur in einem Windkanal mit freier Meß- 
strecke, in welchen sich ein kurzer Widerstandskörper, dargestellt durch einen Dipol, befindet, 
wird zunächst der ebene Fall als Potentialaufgabe mittels konformer Abbildung behandelt. 
Es zeigt sich zunächst, daß es eine Lösung gibt, welche zwar die Randbedingungen erfüllt, 
jedoch an der Ü bergangsstelle Düse —Freistrahl eine unendlich große Geschwindigkeit liefert 
Um dies zu vermeiden, wird eine zweite Lösung überlagert, welche die Düsenmündung eben- 
falls mit unendlich großer Geschwindigkeit jedoch in der entgegengesetzten Richtung um- 
strömen läßt, die Randbedingungen erfüllt und im Innern des Kanals singularitätenfrei ist. 
Einige allgemeine Betrachtungen über die bei Randwertaufgaben mit gemischten Grenzen 
auftretenden mathematischen Fragen werden angeschlossen, welche qualitativ auch die Vor- 
gänge bei dem entsprechenden räumlichen Problem klären. In einem Anhang wird gezeigt, 
daß sich bei der Behandlung der Potentialaufgabe, die Anstellwinkelkorrektur in einem 
Windkanal mit teilweise offener und teilweise geschlossener Meßstrecke zu bestimmen, das 
Unendlichwerden der Geschwindigkeit an den Übergangsstellen grundsätzlich nicht ver- 
meiden läßt. 

Der räumliche Fall wird nach der in der ersten Mitteilung (siehe Fußnote I) angegebenen 
Methode behandelt; dabei wird die Integralgleichung für die Singularitätenbelegung des 
Randes numerisch durch Iteration gelöst. 


Bestimmt man die Anströmgeschwindigkeit experimentell aus dem Düsenvorkammer- 
druck, so muß mit einem etwas veränderten Düsenfaktor gerechnet werden, falls der Wider- 
standskörper sich sehr nahe an der Düse befindet. 

Für einen Körper im Freistrahl ändert sich dabei die Geschwindigkeitskorrektur infolge 
Düseneinflusses erst dann merklich, wenn er sich näher als einen Kanalradius an der Düsen- 
mündung befindet. Die Stelle verschwindender Korrektur für u liegt etwa "ji Kanalradius 
von der Düse entfernt, eine Lage, die bei Widerstandsmessungen praktisch unbequem ist. 
Jedoch könnte man in manchen Fällen den Freistrahleinfluß durch Annäherung des Körpers 
an die Düse wenigstens herabsetzen. Weiter eıgibt sich, daß sich die Zusatzgeschwindigkeiten 
außerhalb der Achse für mäßige Dicken nicht erheblich von den in der Kanalachse herrschen- 
den Zusatzgeschwindigkeiten unterscheiden. 


Anhang. 
Eine Bemerkung zu Randwertaufgaben mit gemischten Grenzen. 


In dem in der vorliegenden Arbeit behandelten Falle gelang es, durch Hinzufügen einer 
im Inneren des betrachteten Gebietes regulären Potentialfunktion zu erreichen, daß die Ge- 
schwindigkeit in allen Randpunkten endlich bleibt, auch in solchen, wo die Randbedingung 
wechselt. Dies ist jedoch nicht bei allen Aufgaben mit gemischten Grenzen möglich, wie im 
folgenden kurz gezeigt werden soll. 

Wir betrachten als einfachstes Beispiel dieser Art die Aufgabe, eine Strömung zu be- 
stimmen, die die Strecke -1<xr<+1 als Stromlinie, den übrigen Teil der reellen Achse 
als Potentiallinie hat. Eine erste Lösung dieser Aufgabe ist offenbar die symmetrische 
ke senkrecht gegen eine Platte, eine zweite die zirkulatorische Strömung um die 

Platte, die allgemeine Lösung ist eine Linear- 
kombination dieser beiden. Beide Strömungen 
tv haben für sich in den Punkten = + 1 unend- 
liche Geschwindigkeit, und zwar ist die Um- 


strömung der Streckenenden bei der ersten 
symmetrisch, bei der zweiten antisymmetrisch 
Mm zur y-Achse (Bild 10). Eine passende Linear- 


kombination beider Lösungen kann also wohl 
in einem der beiden Endpunkte eine unend- 
liche Geschwindigkeit vermeiden, nicht aber 
in beiden zugleich, da die Umströmungsrich- 
tung immer in einem der beiden Enden gleich, 
im anderen aber entgegengesetzt ist. Fordert 
man, wie es bei der Bestimmung der Korrekturen eines teilweise abgedeckten Kanals nötig 
ist, daß die Potentialwerte auf den beiden Teilen der x-Achse von — bis — 1 und von 
+1 bis + gleich sein sollen, so bleibt nur die senkrechte Anströmung der Platte als mög- 
liche Lösung, bei der in beiden Endpunkten unendliche Geschwindigkeiten auftreten. 324 


Bild 10. Senkreehte Anströmung und zirkulatorische 
Umströmung einer ebenen Platte. 
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Untersuchung über die Grenzschicht einer Flüssigkeit über 
einer rotierenden Scheibe bei kleiner | 
Winkelgeschwindigkeitsänderung. 
Von Karl-Hans Thiriot in Friedrichshafen. 
(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


Zur Untersuchung der stationären Grenzschichtströmung nach einer kleinen 
Winkelgeschwindigkeitsänderung wurden die Navier-Stokesschen 
Differentialgleichungen aufgestellt und für den vorliegenden Fall integriert. 
Im Ergebnis sind die Neigungen der sich als logarithmische Spiralen 
ergebenden Stromlinien verglichen mit den Messungen des Verfassers und 
Rechnungen von Cochran und Bödewadt. 


Wenn eine Luftmasse über der rotierenden Erdoberfläche zunächst in Ruhe ist und 
nun unter Verminderung ihrer Höhe auseinandergezogen bzw. unter Zunahme ihrer Höhe 
seitlich zusammengeschoben wird, so nimmt sie dadurch relativ zur Erdoberfläche gleichzeitig 
eine Rotation an, die nach Beendigung der Umlagerung konstant bleibt. Der erste Fall stellt 
den Typus eines Hochdruckgebiets (Antizyklons) dar, der letztere denjenigen eines Tiefdruck- 
gebiets (Zyklons). Bei diesen beiden Strömungsvorgängen ist die Verteilung der Bodenwinde, 
die im ersteren Falle eine Komponente nach außen, im letzteren Fall eine nach innen haben, 
von besonderem Interesse. Diese Vorgänge sollen nun in folgender Weise idealisiert werden. 

Vorhanden sei eine mit einer Scheibe in gemeinsamer gleichförmiger Rotation befind- 
liche volumbeständige Flüssigkeit von der Winkelgeschwindigkeit w,. Die Flüssigkeit sei 
unendlich ausgedehnt und unendlich hoch. Zur Zeit =0 werde die Scheibe plötzlich um 
die Winkelgeschwindigkeit &’ beschleunigt oder verzögert. Die endgültige Winkelgeschwin- 
digkeit sei Vorausgesetzt sei 


o'—=@,— klein gegen 


Es soll nun die stationäre Strömung nach der Winkelgeschwindigkeitsänderung unter- 
sucht werden, wobei das Hauptinteresse sich auf die hier als laminar angenommenen Grenz- 
schiehten konzentrieren wird (der meteorologische Vorgang verläuft zwar mit turbulenten 
Reibungsschichten, die charakteristischen Vorgänge werden aber qualitativ auch bei der 
laminaren Strömung in die Erscheinung treten). 

Wir betrachten die Strömung von dem mit @, rotierenden System aus und benutzen 
die Navier-Stokesschen Differentialgleichungen in Polarkoordinaten (vgl. etwa Thiriot, [1]*). 
Der Rotationssymmetrie wegen werden alle Ableitungen nach g zu Null. Wenden wir noch 
die Vernachlässigungen der Prandtlschen Grenzschichttheorie [2] an, so erhalten wir folgendes 
System partieller Differentialgleichungen zweiter Ordnung. 


(l), 
0 2° 
"r d + = +v, +2 (2), 
mit den Randbedingungen 
für z=0 v=0, 
für 2=o w=—u’r, 


wobei = — klein gegen 


*) Das Schrifttumsverzeichnis befindet sich am Schluß der Arbeit. 
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Erste Näherung. 
Zur Integration der Differentialgleiehungen machen wir den Ansatz 


a 


/ 
wobei 

Gehen wir mit diesen Ansätzen in die Differentialgleichungen ein, so erhalten wir aus 
Gl. (1) die Geschwindigkeitskomponente v,, aus Gl. (2) die Geschwindigkeitskomponente vy 
und aus der Kontinuitätsgleichung ”;. Berücksichtigen wir noch, daß wegen des klein vor- 
ausgesetzten &’ die quadratischen und höheren Glieder von ©’ zur Berechnung der ersten 
Näherung außer acht gelassen werden können, so erhalten wir 


0° 


Differenzieren wir Gl. (5) zweimal nach z und multiplizieren Gl. (4) mit » und Gl. (5) 
mit 2, und substituieren, so erhalten wir 


Diese Gleichung läßt sieh integrieren mit Hilfe des Ansatzes 


v 


Wir setzen den Ansatz (7) in die Differentialgleichung (6) ein und erhalten für i die 
Gleichung 4. Ordnung 


ja E 


Diese Gleichung hat 4 Lösungen: 


+i 
v v 


/ 


Für die Differentialgleichung (6) erhalten wir demnach die allgemeine l.ösung 
Zur Befriedigung der Randbedingungen muß 
für z=x; fi 
für z=0; sein. 
Führen wir die nunmehr erhaltene speziellere Lösung in Gl. (4) ein, so erhalten wir 


dr, 
= 20, (fr) e fıv (r) ) (+a). 


v v 
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Für diese Differentialgleichung erhält man die Lösung 


wo e, und e, die noch von r abhängigen Integrationskonstanten sind. 
Auf Grund der Randbedingungen werden aber diese zu Null. Außerdem wird deshalb 


fu 


fıv 
=. ) 
+ ar ( 


und wenn man die erhaltenen Werte für 4,’ und 4,? einsetzt, bekommt man 


fıv) 

y 
Daraus folgt: 


fn=fiv: 

also 
or 


Damit wird 


und aus der Kontinuitätsgleichung 


cos | | De + sin | 


/ 


Zweite Näherung. 


Um die Differentialgleichungen (1), (2) und (3) genauer integrieren zu können, wenden 
wir ein Iterationsverfahren an, indem wir mit den Lösungen der 1. Näherung in die Diffe- 
rentialgleichungen (1), (2) und (3) eingehen. Wir erhalten dann folgende Differentialgleichungen: 


[073 / \ 
" (sm cos] . . (la) 
2 / 2 
' "z+sin]/ "z-e ' . . (ib). 
y 02° 


Durch Elimination wird dann 


„2 Up 


+40, (ty sin cos|/ "ze 


v 


Die Differentialgleichung hat als partikuläres Integral 


und die Differentialgleichung (la) 


| 

2yro, 


sin e 
v 


m 
942 
us 
IT- 
en 
4), 
5) 
_ 

: 

N 

2) 

/ 

v 

’ 
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Das vy der 2. Näherung wird dann unter Berücksichtigung der Randbedingungen 


-— "2008 ze - 2 
v 
[22] 
v v 


r 
sın — 2 

10, v 


wo die unterstrichenen Glieder die Bestandteile der ersten Näherung sind, und das v, der 
zweiten Näherung 


10, 


wo das unterstrichene Glied der Bestandteil der ersten Näherung ist. Die zweite Näherung 
v, wird aus der Kontinuitätsgleichung zu 


008 2 ze 
=0% 
e + \eos ze 
/ 
a - [47] 
v v 


/ 


sin|/ —ze + 


v v. v v 


Zur Darstellung der Funktionen in Bild 1 machen wir », und vy mit ©’ r und v, mit Yv|o'| 
dimensionslos. Wir erhalten dann 


2 


’ 
= — | — ine 'sin? — 
or 20, 1 10, 
100, 1 
== '— = + cosı 
f: or + 20, 7 1 1 
[2 
|cosne ' sınye 
= = |nceosne"— — sinne ' 
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wobei n=2 \ ®o ist und die oberen Vorzeichen für positives o’ und die unteren für negatives 
v 


1 | 
Gültigkeit haben. Für die beiden Fälle “+7 und m’ u: wurde das Ergebnis auf- 


getragen (Bild 1 bis 3). 


05 
ANSEHEN 
0% 
02 
1/ 
N 
; | 
l l 
or 
| | | 
| — 7 7 + 
| hfü | | 
4 Ale" | | | Bild 1 (links oben). Radiale Geschwindig- 
& 05 T T N T | Me keitsfunktion in Abhängigkeit von der Höhe. 
4 | | 
| | | | | Bild 2 (links). Transversale Geschwindig- 
as | keitsfunktion in Abhängigkeit von der Höhe. 
| | | | 
- 2 5 Bild 3 (rechts oben). Vertikale Geschwindig- 
& keitsfunktion in Abhängigkeit von der Höhe. 
Stromlinien. 
\ Zur Ermittlung der Stromlinien bilden wir aus 
r=e't (logaritlimische Spirale), 
tg rdo 


wobei = unabhängig von r ist. Wir untersuchen die Grenzschicht am Boden, wo z und 


damit auch 7=0 ist. Da sich dann für — = ergibt, muß Zähler und Nenner differenziert 
und dann der Grenzübergang durchgeführt werden. 

j Der Exponent der e-Funktion wird dann gleich 


=—11. 


wobei für o' >0 positives und @’< 0 negatives Vorzeichen gilt. In Bild 4 sind die Ergeb- 

nisse dieser Rechnung dargestellt. Es sind zum Vergleich auch noch die folgenden Punkte 

miteinandergetragen: 1. die Strömung einer nichtrotierenden Flüssigkeit über einer rotieren- 

den Scheibe Fr -1) nach der Theorie von Cochran (3) und nach Messungen des Verfassers 


im rotierenden Laboratorium des Göttinger Instituts, 2. die Strömung einer rotierenden 
Flüssigkeit über der nichtrotierenden Scheibe nach den Rechnungen von U. T. Bödewadt (4) 


(=- x). Die beiden theoretischen Zahlenwerte sind — m =0,83 und 1,26, die gemesse- 


nen 0,80 und 1.2. 


Um bei größeren (positiven und negativen) = -Werten etwa beim stationären Aus- und 
Do 


Anlauf größere Übereinstimmung zu schaffen, müßten noch weitere Näherungen gerechnet 
werden, was sich jedoch als äußerst mühsam erwies. 


Z. angew. Math. Mech. 
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die Grenzschicht einer Flüssigkeit 


11 911 
Cochran gerechnet Ehman-Spiroie 
Thiriof gemessen I — —- 1 
| Cochran \ 
= as 
| j 
oc 
035 03 05 06 07 08 
w' | 
2.2 


Bild 4 (links oben). Spiralensteigung in Abhängigkeit 
Winkelgeschwindigkeitsänderung. 
C=f | | Bild 5 (rechts oben). Ekmanspirale und deren Abweichungen bei 
| Winkelgeschwindigkeitsänderung. 
Bild 6 (links unten). Beiwert der vertikalen Geschwindigkeitsfunktion 
für große Höhe bei verschiedenen Winkelgeschwindigkeitsänderungen. 


von der 


x 


Um ein Bild über die Abweichung von der 
Ekman-Spirale (Geschwindigkeitsvektor abhängig von 
der Höhe) zu erhalten, wurde noch f, als Funktion 
von f, zusammen mit der Ekman-Spirale für die beiden 


1 
=:zr und die entsprechenden Werte von 


Cochran und Bödewadt aufgetragen (siehe Bild 5). 


| [O7 5 5 5 5 
0.1 0,0914 — 0,0892 — 0,1036 0.0958 — 0,0042 + 0,0041 
0.2 0,1633 — 0,1620 — 0.2051 0,1901 — 0,0158 + 0.0155 
0.3 0.2180 — 0,2199 — 0,3028 0.2818 — 0.0335 + 0,0331 
05 0.2846 — 0,2970 — 0,4819 0.4535 — 0,0786 + 0.0796 
1.0 0.2897 — 0,3294 — 0,8124 0,7901 — 0,2134 + 0.2264 
15 0,1988 — 0,2464 — 0,9835 0.9849 — 0,3239 + 0,3573 
2,0 0,1049 — 0,1412 — 1.0456 1.0671 — 0,3910 + 0,4437 
25 0.0402 — 0,0580 — 1,0512 1.0803 — 0,4222 + 0,4871 
3.0 0.0062 — 0,0078 — 1.0364 1.0622 — 0,4316 + 0,5006 
35 — 0,0067 + 0,0145 — 1.0198 1.0367 — 0,4310 + 0,4982 
4.0 + 0,0191 1,0080 1.0160 — 0,4273 + 0.4903 
45 — 0,0064 + 0.0154 — 1,0016 1.0031 042338 + 0,4824 
5.0 — 0,0035 + 0.0094 — 0,9991 0,9971 | +0,47 

Zusammenfassung. 


Um den Einfluß der Winkelgeschwindigkeitsänderung auf den Grenzschichtzustand 
über der rotierenden Scheibe zu übersehen, wurden die Grenzschichtgleichungen für den 
stationären Fall integriert und die bodennahen Stromlinien, die sich als logarithmische 
Spiralen ergeben, berechnet. Weiterhin wurden zusammenfassend die bisherigen Ergebnisse 
im Problem der rotierenden Scheibe in bezug auf das m im Exponent der logarithmischen 
Spiralen für die verschiedenen Fälle dargestellt. 

Nicht verfehlen möchte ich, meinem verehrten Lehrer, Herrn Professor Prandtl, für 
die vielfache Anregung bei der Abfassung der Arbeit meinen ergebensten Dank zu sagen. 
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Halbebene mit periodisch gewelltem Rand. 


Von Constantin Weber in Dresden. 


Unter Benutzung einer konformen Abbildung der Halbebene mit geradem 
Rande auf eine Halbebene mit periodisch gewelltem Rande wird eine ge- 
schlossene Lösung für die Airysche Spannungsfunktion gebracht und die 
Anwendung an einem Beispiele gezeigt. 


Wir bilden die z-Ebene konform auf die ©-Ebene durch 
. . 2 2 22.2.0 
ab, wobei a,,@,,..., «a, ganze reelle Zahlen sind. 


Dann ist 


2 22222.) 


Die obere [-Ebene 7 >0 geht hierbei in den 
oberen Teil der z-Ebene über, der durch eine perio- 
disch gewellte Linie begrenzt ist, Bild 1. Wir wollen 
diesen Teil der z-Ebene als Halbebene mit pe- 
riodisch gewelltem Rand bezeichnen. — 

Dieser Teil der z-Ebene werde in «-Richtung x 
gezogen, so daß für große Werte von y eine gleich- =" 
mäßige Zugspannung Ox,n = 1 entsteht. Bild 1. Halbebene mit periodisch gewelltem Rand. 

Gesucht wird die Spannungserhöhung im 
Grunde der periodischen Kerben. Die Bestimmung der Airyschen Spannungsfunktion werde 
erst für den einfachsten Fall z=[-+ia, ei» gezeigt. Dann ist 


y=nta,e"cosE. 


Die Airysche Spannungsfunktion lautet 
1 


worin =Rf, und Potentialfunktionen von &, »; und damit von #, y sind. 
Wir setzen 


cosE. 


Für große Werte von y geht 9, >b,, und hiermit F>b, +5 y(y-+ec,). Die Spann- 


nung o,„ geht über in = 
0x 


—() werden. 


Für „=0 muß bei unbelastetem Rand sowohl F=0, als auch 5 


Man erhält: 


+ la, +e,)e="cos £) 


1 1 1 


und 


1 


1 1 . 
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Für „=0 müssen die einzelnen Glieder der eos-Entwicklungen gleich null werden. 


Aus F=0 für „—=0 erhält man: 


b,+ Ta, (a, +e)=0, 
1 


1 
bt ga, (a, +te)=0. 
OF 
Aus für erhält man: 


1 1 
+ e)=®, 


1 1 


b, +7 


Die sechste Bedingungsgleichung ist identisch mit der dritten, so daß zur Bestimmung 
der fünf Koeffizienten b,, b,, b, und e,, e, auch nur fünf lineare Gleichungen zur Verfügung 
stehen. 

Hiermit kann die Airysche Spannungsfunktion bestimmt werden, womit auch die 
Spannungen gegeben sind. 

Wir wollen nunmehr zum allgemeinen Fall übergehen. 

Die Halbebene mit periodisch gewelltem Rand wird nunmehr durch die Gl. (2) und (3) 
auf die Halbebene 7 >0 abgebildet. Die Reihen mögen mit den Gliedern sinn und cosn 
abbrechen. Alle Koeffizienten a; mit k>n und k=0 sind gleich null. 


Für die Potentialfunktionen p, und , der Airyschen Spannungsfunktion 


1 
(4) 


machen wir die Ansätze: 


cos nd. (6). 


Wesentlich ist, daß die Reihe von 9, mit cos2nS& und die Reihe von g, mit cos n £ endet. 
Alle nieht angeführten Koeffizienten 5b, und ce, sind gleich null. 


Wir erhalten 2» +1 Koeffizienten für und n+1 Koeffizienten für Aus den Be- 


oF 


dingungen —=0 für dürfen also höchstens +2) Gleichungen zur Bestim- 


mung dieser Koeffizienten folgen. 
Die Spannungsfunktion wird wie folgt umgeformt: 


k=2n hen 


1 Y 
F= > > a), eK" cos (a, + e„)e=""cosh :) 
k=u h=V 


(mit a,—=0), 


ken 

1 


1 
+7 aplayn + 


h=u 


N 


| 
| 
_ 
N 
k=n h=n 
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1 kn k=!n 
1 k=n 
> > (ayn +h)& 
h=0 
k=n h=n 
> > ax (a, —h)E. 
k=0 


Für n=0 muß F=0 sein, hieraus erhält man die Bedingungen 


für die Glieder mit cosO £: 


1 
a, a; (a, + cr) = 0 


k=u 
für die Glieder mit cos(+m&d, 


kn 


1 1 
+7 > ax (ak + Ck- m) 


k=m 


1 ) 


h=m 


für die Glieder mit n<m<2n: 
k=n 


1 ) 
k =0 


Die Ableitung wird: 


F 1 


k=2!n kn 
= 1 
k=0 


1 | 


(k + 


F 
Für „=0 muß sein; hieraus erhält man 


für die Glieder mit cos O&£: 


(2a,+c,) 2 karla, +c)=0 


k=u 


(S), 


(do), 


| 31 
_ 
| 
| 
) 
). 
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h=n 
1 k=n hen 
h=V0 
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für die Glieder mit cos (+ m&), 0<m<zn: 


11), 
k—n hen 
v 
4 (2k- m) a; (a; mt ck m) \ (2h 
4 Bi 4 
k=m h=m 
für die Glieder mit n<m=?2n: 
1 
— mbn— gm N 


Wie man sieht, stimmen die Bedingungsgleichungen (12) mit den Bedingungsgleichungen (9) 
überein. 

Auf diese Weise reduziert sich das System auf (3n +2) Bedingungsgleichungen. 

Das Besondere des Ansatzes besteht darin, daß es hierdurch möglich ist, die Koeffi- 
zienten zu bestimmen und daß bei einer endlichen Anzahl von Gliedern in Gl. (2) und (3) 
auch die Potentialfunktionen y, und 9, nach Gl. (5) und (6) nur endlich viele Glieder besitzen. 

Für die später vorzunehmende Bestimmung der Spannungen im Kerbgrunde werden 
nur die Koeffizienten e, bis ce, benötigt. Zu ihrer Bestimmung addieren wir Gl. (11) und das 
m-fache der Gl. (8) und erhalten für Oo<m=n die n-Gleichungen mit k-m=h—m=t: 


t=n—m 


... . . (IB). 


Für n=5 lauten diese Gleichungen: 
m=1: 2a,+te— (a,+e,) +a,(a,+e,)] 
-2[a, (a, +e,) +a,(a,+ c,)] 
-3la, (a,+e,) +a,(a,+ 
—4[a, (a,+ec,) 
m—2: 2a,+0,— 
2la,(a,+e)+a, (a,+c,)] 
Bla, (a, + e,) +a, (a, + e,)] = 0, 
m=3: 2a,+c,— 
- 2la, (a, +, (a, + 0,)]=0, 
m—4: 2a,+e,— J[a,(a,+e,) +a,(a,+ec)]=0, 


m=5: 2a,+e,—V 


2 
Die Spannung im Kerbgrunde wird da dort o,,=0 ist, so 
folgt: 
1 
(15). 


(Diese einfache Bezieliung gilt natürlich nur für die Stelle 2=0, =.) 


Beispiel. Für eine Halbebene mit periodi- 
schen Halbkreiskerben nach Bild 2 wollen wir die 
Spannungserhöhung im Kerbgrunde finden. Die, 
Periodenlänge ist 1-27, der Halbmesser der Kerben 


% 
77 


Bild2. Halbebene mit periodischen Halbkreiskerben. 


a 
sei o0=7. Wir ersetzen die Begrenzung der Halb- 


ebene durch eine Begrenzung nach Gl. (2) und (3) mit „=0 und nehmen n=4, so daß wir 
vier Koeffizienten a,, a,, a, und a, wählen können. Zu ihrer Bestimmung nehmen wir 
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des Krümmungshalbmessers im Kerbgrunde beim Übergang zum Nachbarpunkte gleich Null. 
Eine vierte Bedingung wählen wir nicht, sondern suchen aus der Schar der so erhaltenen 
Begrenzungslinien gefühlsmäßig die günstigste aus. 


folgende Bedingungen: Kerbtiefe t— Krümmungshalbmesser 0, = 5: zweite Ableitung 


Die drei gewählten Bedingungen geben folgende Gleichungen für den Rand 7=%: 


Aus der dritten Bedingung folgt: 


Setzt man hierin (3 .) u; 7 nach Gl. (17) ein, so vereinfacht sich die letzte Glei- 


chung in: 


-3[; I\=-- +16a, 


"br 
+8a,+ 27a, +6ta,)(l-a, -2a, -3a, —ta,) 
3la, +ta,+Y9a, + 16a,” 


(18). 


Die Gleichungen sind naeh Wahl eines Koeffizienten durch Probieren gelöst. 


Es wurden folgende Werte gewählt: 


a,=--0224, 
a,= 0,003. 
Bild 3 zeigt die sich hieraus ergebende Begrenzungslinie. 
Nunmehr berechnen wir nach Gleichungssystem (14) die 
Koeffizienten e,.c,, c,,c,. Vorübergehend schreiben wir a, +c, =d,, 
a,+06,=d,, a,+e,—d,, a+ec,=d, und erhalten die Glei- 
chungen : 


(A323 3) 


Bild 3. Gewählte Begrenzung. 


1224d, 083 d,+-0,457d, 0,135d, +0,70 0, 
0,008 + d, +0,05 — 0, 
d, - 0,003 — 0, 


nit der Lösung d, = — 0,474, d,— 0,170, d,— 0,042, d, 0,003 und hieraus ec, — - 1,214, 
6,034, 0, ec, 0,006 . 


Die Spannung im Kerbgrunde wird nach Gl. (15): 


0,663 
’ | 


Man erhält eine 2,13-fache Spannungserhöhung, während bei einer einmaligen Halb- 
kreiskerbe eine 3,07-fache Spannungserhöhung entsteht. 323 
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Von den freien Schwingungen eines Kreiselpendels 
bei endlichen Ausschlägen. Il. 
Von !!. T. Bödewadt in Göttingen. 


In der ersten Mitteilung war das Gleichungssystem des Kreiselpendels durch 
ein einfacheres ersetzt worden, um daran die Methode darzustellen: Entwick- 
lung der periodischen Eigenlösungen in Fourierreihen nach der Zeit, Ent- 
wicklung der Fourierkoeffizienten in Potenzreihen nach der Schwingungs- 
weite. In dieser Mitteilung wird das strenge Gleichungssystem aufgestellt 
und nach derselben Methode behandelt: ein Zahlenbeispiel ist hinzugefügt. 


Die erste Veröffentlichung des Verfassers über diesen Gegenstand ') enthielt insofern 
noch eine Unrichtigkeit, als die dort angenommenen Gleichungen I (1.1) nicht eigentlich die 
Gleichungen des Schiffskreisels waren’), sondern aus ihnen erst durch Symmetrisierung ent- 
standen sind. Dadurch wurde die Darstellung der mathematischen Behandlung vereinfacht; 
doch ging andererseits die Beziehung zum physikalischen Gegenstand halbwegs verloren, ins- 
besondere hatte die Durchrechnung des Zahlenbeispiels in I, $ 9 nicht viel Sinn. In dieser 
zweiten Veröffentlichung soll darum ein besser zutreffendes Gleichungssystem abgeleitet und 
nach der in I dargestellten Methode behandelt werden. Da die grundsätzlichen Überlegungen 
ungeändert bleiben *), sollen hier nur die neuen Entwicklungen formelmäßig angegeben und 
das Zahlenbeispiel neu durchgerechnet werden. 


Wir legen unserer Betrachtung ein mechanisches Modell zugrunde, das aus drei in sich 
starren Teilsystemen besteht: aus Schiff, Kreiselrahmen und Kreisel. Das Schiff soll eine 
Symmetrieebene haben. In dieser liegt waagerecht die Rollachse, um welche das Schiff 
pendeln („rollen“) kann. Diese soll raumfest sein, d. h. Kursänderungen und Stampfbewe- 
gungen des Schiffes bleiben außer Betracht. Im Schiff liegt senkrecht zur Symmetrieebene 
die Lagerachse. Um diese ist der Kreiselrahmen drehbar, der, wieder die Kreiselachse trägt, 
um die sich der Kreisel drehen kann. Die Kreiselachse soll stets in der Symmetrieebene des 
Schiffes liegen und die Lagerachse schneiden. — Wir führen mehrere kartesische Koordinaten- 
systeme ein: ein im Raume festes System # yz, ein schiffsfestes System «vw, ein rahmen- 
festes System £7{ und ein kreiselfestes System U VW. In der Ruhelage des Gesamtsystems 
fallen alle Koordinatensysteme zusammen, der gemeinsame Nullpunkt ist der Schnittpunkt 
von Lagerachse und Kreiselachse; das Steuerbordende der Lagerachse gibt die Richtungen 
+0, +u, +U, waagerecht voraus laufen parallel zur Rollachse im Abstande unter 
ihr die V-Achsen, und lotrecht nach oben gehen die +2, 
Achsen. Außerhalb der Ruhelage fallen stets die + W und die + £-Richtungen in der Kreisel- 
achse zusammen, die +u und die +$-Richtung in der Lagerachse, während die +, und 
+ r-Riehtungen immer parallel zu der raum- und schiffsfesten Rollachse bleiben. Zur Be- 
schreibung der gegenseitigen Lage der drei Modellteile zueinander benutzen wir drei Winkel: 


y=<(+z, +ır) im Drehsinne (+2...+%): 

(+-w,+£) im Drebsinne (+v... +1): 

—+U) im Drehsinne (+8...+ 9). 
_ Der Kreisel soll symmetrisch sein, sein Schwerpunkt also in der Kreiselachse liegen und sein 
Trägheitsmoment um die Lagerachse unabhängig von sein. Der Rahmen soll symmetrisch 


zur Ebene „7=0 sein; und der Schwerpunkt des Gesamtsystems soll in der Symmetrieebene 
des Schiffes liegen. 


Dieses Modell hat dann drei Freiheitsgrade y, 9, g (Translationen kommen nicht in Be- 
tracht). Seine Lagrangesche Funktion 


1) Z. angew. Math. Mech. Bd. 20 (1940), S. 218 bis 234. — Im folgenden mit „I“ zitiert. 
2) Herrn E. Schmid habe ich dafür zu danken, daß er mich hierauf wieder aufmerksam machte. 


3) Lediglich die Bemerkungen in I, $ 3 sind dahin zu ergänzen, daß die Konvergenzgrenze aus physikalischen 
Gründen wahrscheinlich bei einer Schwingungsweite von M=90®° liegen wird, da darüber hinaus bereits zyklische 
(nieht mehr periodische) Lösungen anftreten. 


4 


n 
h 


en 
he 
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hat unter Fortlassung einer unwesentlichen Konstanten bei Umrechnung auf die Veränder- 
lichen 9, das Aussehen 


= 5 + Esin?d — cos 


CO sindyg + — qgh)ecos y-+ gheos y cos 


Hierin bedeutet J das Trägheitsmoment des Gesantsystems um die Rollachse bei 9—=0; 
E ist das Moment von Kreisel.mit Kreiselrahmen; m die Masse, my das 
Gewicht von Kreisel und Kreiselrahmen; — h die S-Koordinate des Schwerpunktes von Kreisel 
mit Kreiselrahmen; ! der Abstand zwischen Lagerachse und Rollachse; j das Trägheitsmoment 
von Kreisel mit Kreiselrahmen um die Lagerachse, € das Trägheitsmoment des Kreisels um 
seine Achse, X das Auslenkungsmoment \Fld m des Kreiselrahmens, 9 das Gewicht des 
Gesamtsystems, // die metrazentrische Höhe des Schiffes, d. h. der Abstand des Schwer- 
punktes des Gesamtsystems (im Zustande 90) von der (durch des Metazentrum gehenden) 
Rollachse. 

Für den Drehimpuls um die Kreiselachse ergibt sich Konstanz: 


= sindy)—=N. 


Man kann also die Veränderliche g und ihre Ableitung eliminieren und muß dann die Funktion 


+(E -C)sin’d mhl(l - cos#)} 


+5 Nsindu +(QH  qgh)cos y + qh cos y cos 


verwenden. Läßt man hier bei y” die von # abhängigen Glieder fort, vernachlässigt also 
—E (das ist im wesentlichen das Trägheitsmoment des Kreisels um die Lagerachse) und 
m hl. gegen .J, ferner K gegen .J und j, so bleibt 
Somit werden die Bewegungsgleichungen, wenn man in der ersten noch qh gegen QH ver- 
nachlässigt: 


+ OH sin y 


jÖ-+gh 


Die Gleichungen I (1.1) waren hieraus entstanden, indem in der zweiten Gleichung cos y 
in 1 und cos # in cosy geändert wurden; so wiesen die geänderten Gleichungen dasselbe 
Maß an Symmetrie auf wie im linearen Grenzfall. Ohne diese Änderung ist es zwar auch 
noch möglich, die Gl. (1) in y zu linearisieren und dadurch auf eine Gleichung vierter Ord- 
nung für 9 allein zu kommen, doch wird diese dann so viel verwickelter gegenüber der 
Gleichung I (3.8), daß sich die Elimination nicht lohnt. Wir gehen also nach der in Is 10 
angedeuteten Methode vor und brauchen uns dann auch nicht in y auf kleine Winkel zu 
beschränken. 


Schreibt man in den Gl. (1) » für y und y für 9, so gelangt man durch Division mit J 
und j zu der Form 


+ heosyy—d, | 
| 

Wir suchen die reinperiodischen Lösungen der Form 

3) 

die die Bedingung er 


erfüllen und für M>0 in Form und Periode in die Lösungen des linearen Grenzfalles über- 
gehen. Für die in (2) auftretenden Funktionen von «= und y ergeben sich Entwicklungen 
der Gestalt 


h. 
h 
| 
lt 
4 5 
e 
N 
4 
n | 
ff 
» 
t, € 
IS 
n 
& 
]- 
d 
“ 
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sin.e= D’MC,sin (2r+1)owt, 


sin y= D,cos +l)wt, 

cos siny= ME,cos (27 +1) mt, 


Einsetzen von (3) und (5) in (2) liefert für r—=0,1,2,...: 
A,— a (2r+l)ohD,=0, 
@r+l)okF,—0 | 


Unter der i. A. zutreffenden Voraussetzung, daß keine „Grenzresonanz* stattfindet (vgl. 
I $ 4, haben die für hinreichend kleine Werte von M möglichen Potenzreihenentwicklungen 
der Größen A,, B,. C,. D,, E,, F, die Form 


Dabei kann man setzen 
Da»=Buyt | 8) 
(9), 
Bart us; Far = Aut fm 
und die Zusatzgrößen Ayy: Eny: fu, können aus den An. und Bmu mit m <n berechnet 
werden. Die Formeln dafür bis n=4 sind unten in Tafeln zusammengestellt. — Auch die 
Frequenz wird als Potenzreihe in A angesetzt: 
Es ist also 4 0,0, @,n,, insbesondere am Anfang 9, Wir setzen 
ferner 
n—| 
Durch Einführung der Potenzreihen (7), (9) in (6) entsteht : 
n— 
(il). 
n—v 
2 — (2r+1)k N | 
- m m,» - m En-m,»r — 
=@ 


Für = 0 bleiben die Gleichungen*) des linearen Grenzfalles: 


* 


aus denen bekanntlich 
+h: k) =0 


und damit < min (a?, bzw. max (a?, folgt, sowie 


Für >0 sind zunächst in den Fällen die Bu, aus dem 
Gleichungspaar (11) zu bestimmen. Nach diesen Unbekannten geordnet heißen die Gleichungen: 


— a} + Gyr | 


wobei @,„, und H,„, Abkürzungen sind für?) 


(12), 


4) Dabei ist benutzt worden, daß aus (4) Boo=1 folgt: ferner wurde zur Abkürzung AJoo = .I geschrieben. 
5) Die in (13) rechts stehenden Summen sind leer fürv—n. 


| 

XUM 


(). 


(6). 


gl. 


zen 
(S), 
net 
die 


zen 


10). 


il). 


N 
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n—ı 
(21 +1)? Au-m,», +0’ Ca du» Om D.-m. 


n—v n—v (13). 
= (2 +1)? Im Bu- mvt b’en-+ (2 1) k fn „+ Fn- m, | 
Die Determinante des Gleichungssystems (12) ist 
1 


Unsere Voraussetzung, daß keine Grenzresonanz stattlindet, sagt das Nichtverschwinden der 
Determinante (14) aus. 


Setzen wir also für v>0O 


a,={@r +1? 
(15) 
„= | 

so ist 


Bıv=ß, 9, 


die Auflösung von (12). Sind die Größen A,„,. B,, für v»—=1,2,..., » berechnet, dann er- 
gibt sich aus (4) 
n 


Für »=0 enthalten die Gl. (11) die beiden Unbekannten A,, und ,„: 
Aus — a) + A- h)=@G,+heo, Bus: | 


(18); 
Aus (— ko,) kA)= - Fo) 
darin bedeuten @,„o. folgende Ausdrücke: 
n—1 
—=1 
m (19) 
n—1 n 
— 9 „ tbe.tk fn.t+ Om Fa-m, | 
Die Determinante des Systems (18) ist 


verschwindet also nicht; denn sonst müßte &} als geometrisches Mittel zwischen «a? und b? 
liegen, während beide Werte von »? außerhalb des Intervalles (a?...b?) liegen. Wenn man 
nun setzt 


el), 
so ist die Lösung von (18): 


Der Frequenzkoeffizient w,„ hängt also nieht von B,„, und damit nicht von den Amplituden- 
koeffizienten der gleichen Stufe ab. 

Es sind jetzt alle Formeln angegeben, um schrittweise die Koeffizienten der Reihen (3), 
(7), (9) für jede Ordnung berechnen zu können; nur die Ausdrücke für ey», du», Enns fu» 
fehlen noch. Bis n—=4 sind die ce, e, f in den nachstehenden Tafeln 1 bis 3 angegeben. Für 
die d„,„ waren die Ausdrücke bereits in I (Tafel 1°)) mitgeteilt worden; man hat jetzt nur 


6) Leider enthält jene Tafel einen Fehler: in Ego muß der Zähler des Faktors von C;o eine 7 sein (nicht 1). 


| 
| 
B..= — 
| 
2 
| 
m 
lem 
en: 
12), 
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noch die dortigen Bezeichnungen C',y, Ey», in B,,. d„, abzuändern. Jedoch lassen sich die 
Formeln für die d,, auch aus der untenstehenden Tafel I für die e„, gewinnen, indem man 
dort A,„, dureh (- 1” B,, und e„, durch (—1)”d,, ersetzt: die so erhaltenen Ausdrücke 
können dann durch Ausnützung der Beziehungen B,,=1 und (17) vereinfacht werden. — 
h,.h,,...h,, und H', H®,... H' sind Hilfsgrößen, die eingeführt wurden, um die Berechnung 
zu vereinfachen. Aus dem gleichen Grunde wurden in Tafel 2 auch die D,, und d,, ver- 
wendet. 
3. 

Nach den oben abgeleiteten Formeln wurde das in I $ 9 gewählte Zahlenbeispiel erneut 
durehgerechnet. Es sind also für die in (2) vorkommenden Konstanten folgende Werte an- 
genommen worden: 

#056; k=221. 


Tafel 1. Formeln für die durch (4), (5). (7). (8) erklärten c,,, 


n 
n=r n=r 
1 1 
> 4. 
1 1 
7 


1 1 


46080 ° ° 
1 1 
4 4: Au 34° 1+%Aı + Mm Aıt 


1 
(Ay — + 4 A{Ay Ass + Ayı — + N; — 4 
1 
+ 8 Ayo + (2 Ayo — 3 + 6.4,,) 
1 
4 (Ay - 2 Ad) + Au (An — + (Ay — 4 Ay + 6 


1 


1 
34° — 8 + A) + 12 Ad An — 3 + 21 — Au) — 84 h,. 


1 
3 An + N Ant Ad) + Au 2 + (10 Ay, — 5 Ası + 


1 
+ 8 Ayo Ay —2A},) + hz\ 10 — 19 Ayo Au 6 A) 21 Ay) 120 4, 


XUN 
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ie Tafel 4 zeigt die zunächst erhaltenen Koeffizienten für die Potenzreihen und läßt erkennen, 
ın daß die Konvergenz bei nicht zu großen Ausschlägen recht gut ist. Nimmt man als Aus- 
ke schlagweite zum Beispiel M=60° an, so findet man in Tafel 5 die erhaltenen Werte für die 
Be Frequenz und die Amplituden der Fourier-Teilschwingungen. Ein Vergleich” dieser Zahlen 
18 mit denen, die man aus dem linearen Ansatz 


hy=0; 


Tafel 2. Formeln für die durch (4). (5). (7). (8) erklärten e,„, 


ut 
Zr 
1 1 
13 6 > 4 22 | 8 
1 1 
4 8 4 
ho): 
+ h, Du + hu: 
= da + (Du —2D,) Dot Mm: 
Dot Mıs- 


1 
— 2 Aı- + Ar) + Ay Ass + Au (Ast Az} - 
+ + An An) +2 Ad An An Arı)+ {A An + (An + Ası))- 


— 23 + (9 Ara + (4 Ada 6 Ayo) + — 3 — + {4 (Ayo — + — Az) — 


2D,) Da Dat Ms Dot 


du +h (Do —2 Da + Pa) + Au Dat Dat Dit Dot 


angew. Math. Mech. 


40  bBödewadt, Schwingungen eines Kreiselpendels bei endlichen Ausschlägen »2 Nr. 1 Febr. 1943 
Tafel 3. Formeln für die durch (4). (5). (7), (8) erklärten /,» 
dr 
nr 


1 3 
In 24 A. lo = 8 4 
5 

H=+Botgg- H? = 48° 

3 
Al. 

1 1 
In = Au + Al. 
Pl 4 \ 2 « 

1 a 8 1 

1 5 1 

1 3 3 
5 (Aut An) +3H? Au+H® An — H® AN. 
= H® — H® — H4—=H, — 

720° 180 
—B, + Bu Ba +13 
24 1920 072580480 ° 


— 


—?2 B;; B;. + Bo 


+ Bio (Ba 


5 1 , 
35 But 5 Br Bic+ 


+5M A, 3IP A + 
i J > Loy» N 13 | 
Tdu-y +35 +3H' 4, +3H® An — Au +H 
1 7 3 \ 


5 5 


XUM 
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Tafel 4. 
Die Potenzentwicklungskoeftizienten A». B,» und o, für die Frequenz und die 
Amplituden der Teilschwingungen bei a? 0.62. b? 0,56, h= 0.06, 221. 


Präzession Nutation 

0.155192 —+ 3.7830 

©, = + 0.026404 ©, = — 0.883400 

©, —= + 0.004126 + 0.024720 

0, + 0.000678 = — 0.005758 

= 0.000160 — 0.001311 
Aoo = — 0.015626 Bao = + 1.000000 Aoo = + 0.016513 Bao = + 1.000000 
Ayo = — 0.000136 Bio = — 0.050215 40 = + 0.002452 Bio = + 0.019942 
Aso = —- 0.000153 B.0 = — 0.008062 Aso = + 0.000613 Bao = + 0.003894 
Ayo = — 0.000043 Bao = — 0.000546 Azo = + 0.000179 Bao = + 0.000937 
Ayo = — 0.000018 = + 0.000275 Age = + 0.000057 + 0.000256 
= — 0.000592 = + 0.050215 Au: 0.000326 0.019942 
Ası = -- 0.000113 Baı = + 0.002284 Ası 0.000074 Ba, 0.004614 
Azı = 0.000018 0.001117 Ası 0.000018 Bzı 0.001181 
= 00000005 = — 0.000428 An 0.000005 0.000331 
Ass = — 0.000055 Bas = + 0.005778 Ass = + 0.000012 Bas = + 0.000720 
Asa : 0.000017 By» + 0.000804 = + 0.000004 = + 0.000274 
Ars = — 0.000003 Bis 0.000200 = + 0.000001 Bas + 0.000090 
= 0.0000059 + 0.000859 = 0.0000005 — 0.000031 
= 0.0000028 Bis = + 0.000202 = — 0.0000003 Bis = — 0.000017 
Aya = 0.0000005 = + 0.000150 Aya = + 0.0000000 Bir = + 0.0000015 


Tafel 5. 
Frequenz und Teilamplituden für die Schwingungsweite M = 60° 
bei @=0.62, 5=056, A=0.M, 221. 


Präzession Nutation 
0.1902 sek”! o= 2.85 sek! 
T=2n|/o = 33,0 sek =2n/o=2.21sek 


Ao = — 0.01604 
A, = — 0.00081 


MA, = —0.% 
MA, = —0.05° 
MA, = — 
MA; = — 0.00 
MA, = —0.00° 


+ 0.02026 
Aı = — 0.00048 
As — + 0.00002 
As 0.00000 
Ay = + 0.00000 


MA, — + 1.22 
MA, — — 0.039 
— + 0.00 9 
MA, — 0.000 
MA, — + 0.00 0 


As 0.00009 
= — 0.00001 
Ay = — 0.00000 
Bo = + 0.9349 


B, = + 0.0557 
— + 0.0077 
— + 0.0014 
= + 0.0002 


MB, = + 56.1° 
MB, = + 33° 
MB, =+ 0.5° 
MB, =+ 0.1° 
MB, = + 0.00° 


Bo = + 1.0282 
B, = — 0.0295 


— + 0.0014 
B; = — 0.0001 
B, = + 0.0000 


MB, = + 61.70 


MB, = — 
MB=+ 01° 
MB, = — 0.0° 
MB, = + 0.0° 


für dieselbe Schwingungsweite = M = 60° erhalten würde, führt zu folgenden Ergebnissen: 


Die Frequenz der reinperiodischen Lösung ist bei der Präzession in Wirklichkeit um 
22.6 v.H. größer, bei der Nutation um 25.0 v.H. geringer als nach der linearen Theorie; man 
hätte dort die Umdrehungszahl des Kreisels, mithin die Größen h und k nur mit 799 v.H. 
bzw. 72.6 v.H. ihres tatsächlichen Wertes einsetzen dürfen, um bei gleichen Werten der 
Konstanten a, b die Frequenz richtig zu erhalten. Der Schlingerwinkel des Schiffes 


max 


ist mit 0.92° bzw. 1.25° um 2.0 v.H. geringer bzw. um 25.7 v.H. größer, als die lineare 
Theorie ergibt; das ist freilich bei der Kleinheit dieses Winkels belanglos. — Die Beschleu- 
nigung der Kreiselschwingung im Umkehrpunkt, — (0), ist nach I (9.1) bei der Präzession 
2.59 mal, bei der Nutation 0.45 mal so groß wie der Wert, den der lineare Ansatz liefert; der 
Umkehrruck beansprucht das Schiff also weit stärker, als die lineare Rechnung erwarten - 
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Beiträge zu den Interpolationsverfahren der numerischen 
Integration von Differentialgleichungen 1. und 2. Ordnung. 
Von Lothar Collatz in Karlsruhe und Rudolf Zurmühl in Darmstadt. 


Aus bekannten Abwandlungen des Ada ms- Verfahrens werden zwei Inter- 
polationsformeln für Differentialgleichungen 1. und 2. Ordnung ausgewählt, 
welche sich durch Einfachheit und Genauigkeit auszeichnen. Für diese 
Formeln, ron denen die für Gleichungen 1. Ordnung mit der Simpson- 
schen Regel übereinstimmt, werden Fragen der Handhabung. der Schritt- 
bemessung und der Fehlerfortpflanzung besprochen. 


Für die numerische Integration gewöhnlicher Differentialgleichungen haben die nach 
Adams!) und Störmer?) benannten Verfahren’) in steigendem Maße Eingang gefunden, 
wobei in neuerer Zeit die Interpolationsverfahren in den Vordergrund getreten sind, welche 
dureh iterative Verbesserung die Genauigkeit der Näherung wesentlich steigern. Im folgenden 
greifen wir aus den verschiedenen im Schrifttum bekannt gewordenen Abwandlungen?) zwei 
Interpolationsformeln für Differentialgleichungen 1. und 2. Ordnung heraus, welche die ersten 
Glieder der ein- bzw. zweimal integrierten Stirlingschen Interpolationsformel benutzen und 
der Form nach an Verfahren von Nyström und Störmer anknüpfen. Die Formeln sind 
nieht neu, doch bisher zu numerischen Rechnungen anscheinend nur wenig herangezogen 
worden, obgleich sie einfache Handhabung mit guter Genauigkeit verbinden. Für diese 
Formeln werden Fragen der Handhabung und der Fehlerfortpflanzung untersucht’). 


Teil I (Differentialgleichungen 1. Ordnung) behandelt zunächst Handhabung der Inter- 
polationsformel, stufenweisen Aufbau der Ausgangswerte, Bemessung und Abstufung der 
Schrittweite, Beschleunigung der Konvergenz und andere praktische Fragen. Auf die Fehler- 
abschätzung braucht hier nur kurz eingegangen zu werden, da diese für Differentialgleichungen 
1. Ordnung mehrfach durchgeführt worden ist®). Es wird nur das Ergebnis mitgeteilt, welches 
für das hier beschriebene Verfahren günstiger ausfällt als für die Adam ssche Interpolations- 
formel. 


Teil II (Differentialgleichungen 2. Ordnung) geht nach Darstellung des Verfahrens und 
der Anfangsiteration ausführlicher auf die Fehlerfortpflanzung ein. Soweit die Gedanken- 
gänge denen der Arbeit von G. Schulz’) ähnlich verlaufen, können sie hier kurzgefaßt 
werden. Die dortigen Ergebnisse werden dahin ergänzt, daß auch die Fehler der Anfangs- 
werte abgeschätzt, die Interpolationsformel anstatt der Störmerschen Extrapolationsformel 
zugrunde gelegt und in der Differentialgleichung 7 = fit, y, „) eine Abhängigkeit der Funk- 
tion f auch von % zugelassen wird. 


I. Verfahren für Differentialgleichungen 1. Ordnung. 


1. Näherungsaufbau nach einer Interpolationsformel. 
Zu einer Diflerentialgleichung 1. Ordnung 


ı) F.Bashforth, J.C. Adams: An attempt to the test the theories of capillary action. Cambridge 188, 
S. 15 bis 21. 

2),C, Störmer: Sur les trajeetoires des eorpuseules eleetrises dans l’espace sous l’action du magnetisme 
terrestre avee application aux aurores boreales. Arch. Sei. phys. nat. Bd. 24 (1907), S. 221 bis 247. 

3) Vgl. etwa Günther Sehulz: Formelsammlung zur praktischen Mathematik. Sammlung Göschen Bd. 1110. 
Berlin und Leipzig 1937, S. 116 bis 126. 

4, E. J. Nyström: Über die numerische Integration von Differentialgleichungen. Acta Soc. Sei. fennicae 
B1. 50 (192%), Nr. 13, 8. 29 bis #3. — A. Bennet, W. E. Milne, H. Bateman: Numerical integration of differential 
equations. Bull. nat. Res. Coune. Nr. 92. Wash. 1933, S. 28, 71 bis 07. — W. E. Milne: Amerie. Math. Monthly Bd. 33 
(1926), S. 455 bis 460 und Bd. 40 (1933), S. 322 bis 327. — E. Lindelöf: Remarques sur !’intögration numerique des equati- 
ons differentielles ordinaires. Acta Soe. Sci fennicae A Bd. 2 (1938), Nr. 13. — Sch. Mikeladze: Über die Integration 
von Differentialgleiehungen mit Hilfe der Differenzenmethode. Bull. Acad. Sei. USSR, Ser. Math. Nr. 5 bis 6 (1939), 
S. 627 bis 642 (Russisch). 

5) Die Arbeit entstand im Anschluß an größere numerische Rechnungen, welche am Institut für Praktische 
Mathematik der Technisehen Hochschule Darmstadt durchgeführt wurden. 

6) R. v. Mises: Zur numerischen Integration von Differentialgleichungen. Z. angew. Math. Mech. Bd. 10 (1930), 
S. sı bis 92. — G. Sehulz: Interpolationsverfahren zur numerischen Integration gewöhnlicher Differentialgleichungen. 
Z. angew. Math. Mech. Bd. 12 (1992), 8. 4 bis 5%. — W. Tollmien: Über die Fehlerabschätzung beim Adamsschen 
Verfahren zur Integration gewöhnlicher Differentialgleiehungen. Z. angew. Math. Mech. Bi. 18 (1938), S. 83 bis 90, 

?7)G. Sehulz: Fehlerabsehätzung für das Störmersche Integrationsverfahren. Z. angew. Math. Mech. Bd. 14 
(1934), S. 224 bis 234. 


2 


& 


Nr. Febr. 198 Collatz u. Zurmühl,. Differentialgleichungen 1. und 2. Ordnung 43 


welehe der Lipschitzbedingung 


zenügt (wobe: M durch den größten Betrag der Ableitung x - ol innerhalb des in Betracht 


kommenden Bereiches ersetzbar ist), seien bereits eine Anzahl von Näherungswerten y; an 
(den gleiehabständigen Stellen #; mit der Spanne Yt=h und damit auch die Ableitungen 
n—2,n —1,n) bekannt. 


Zur Neurechnung einer Näherung y,4+, für den exakten Wert y(f,+,) bedienen wir uns 
der Gleichung 


welche die ersten Glieder der integrierten Stirlingschen Interpolationsformel *) 


benutzt und mit der Simpsonschen Regel 


übereinstimmt. Dabei bedeuten die aufsteigende Differenzen, also I f„=fu fun... Im 
Differenzenchema, Tafel 1, in welchem zweckmäßig die mit » multiplizierten Funktionswerte 
und deren Differenzen tabuliert werden, benutzt ıman «ie umrandeten Werte, versehen mit 
den am Kopf des Schemas vermerkten Faktoren 2 und 1/3. 


Tafel 1. Rechenschema zur Interpolationsformel (1.3). 
2 
3 
Y hy=hf 
hf 
KV 
hlu-ı 


Ha 


In (1.3) ist zunächst die zweite Differenz I ?f,+, der neuen Zeile a+1 noch nicht be- 
kannt. Sie läßt sich jedoch, wenn das Verfahren erst einmal in Gang gekommen ist, aus der 
Reihe der bisher bekannten 1? f-Werte durch Extrapolation leicht schätzen®) (etwa durch 
Aufzeichnen der Y*f oder der über f). Ausgehend von einem so geschätzten 
at man die gesuchte Näherung in fortschreitender Verbesserung nach der Iterationsvor- 


S) Vgl. etwa G. Schulz: Formelsammlung®), S. 856. - 


‚ »,;, Hat man keinen Anhaltspunkt zum Sehätzen des neuen Wertes r? /nzı1, so kauın man statt dessen den Wert 
(r? fun +73 fn) benutzen. Die Formel geht dann in die nach r3 f„ abgebrochene Extrapolationsforme] von Nyström 
über. Vgl. E.J. Nyström, a.a. 0.), 8. 32 bis 33, Gl. 
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fit, +1» ı) 
Die Iteration kommt bei genügend kleiner Schrittweite A nach ein bis zwei Durchrechnungen 
zum Stehen. 
2. Iterative Bestimmung der Ausgangswerte. 


Zur Ermittlung der zum Beginn der Rechnung erforderlichen Ausgangswerte bedient 
man sich mit Vorteil eines stufenweisen Aufbaus von Rohwerten mit anschließender iterativer 
Verbesserung, etwa nach folgendem Schema’): 


A. Anfbau von Rohwerten: 


1. Vorläufige Rohwerte auf Zeile 1: 
(y)=M+hfo; daraus: 
—h- 


| 
| 

2. Verbesserte Rohwerte auf Zeile 1: | 
| 
| 
| 


f+ ; daraus: fl". 


> 


3. Rohwerte auf Zeile 2: 


4. Verbesserte Rohwerte auf Zeile I und 2 


+hl2 ” +3 
5. Rohwerte auf Zeile 9: 
Ausgehend von diesen Werten der Näherungsstufe r 1 folgt nun: 


B. Iterative Verbesserung: 


| 


Dabei wechseln ab: 
1. Verbessern der drei y-Werte nach (1.6), 
2. Verbessern der Funktionswerte und der Difterenzen, 


bis sich die Werte innerhalb der mitgeführten Stellenzahl nicht mehr ändern. — Für geringere 
Genauigkeitsansprüche und genügend kleine Schrittweite genügt die Anfangsiteration mit zwei 


10) Zur Bereehnung von yı vgl. G. Schulz: Formelsammlung), S. 87; ya und 43 werden im wesentlichen nach 
Formel (1.3) bestimmt. 


2 
| 
| 
| 
| XUM 


ach 
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Werten y, und »,, also eine Rechnung innerhalb der ge- Tafel 2. Schablone zur Berechnung der 

. Ausgangswerte nach (1.6) und zur fort- 
strichelten Umrandung. Auf A 3 folgt dann unmittelbar laufenden Rechnung nach (1.3). 
die iterative Verbesserung nach B. 


Die Anfangsiteration kann durch Schablonen nach y Irw 
Art von Tafel 2 erleichtert werden ''), wobei die im Dif- 7:77 Ei 2 
ferenzenschema zu verwendenden Werte in (hier schraf- 77 
fiert gezeichneten) Schablonenfenstern, versehen mit den = VZZA3 
zugehörigen Faktoren, sichtbar werden. Der obere Teil 7 2- % a7 

der Schablone dient lediglich zur Anfangsreehnung von — 

der untere dagegen außer zur Anfangsbereehnung von 9, 


und 9, auch zur fortlaufenden Rechnung nach Gl. (13). 


3. Bemessung der Schrittgröße. 
Die Schrittgröße h ist so klein zu halten, daß 
1. die Iteration nach (1.3) genügend rasch konvergiert, und daß 


+» 


2. das in der ausführlicheren Stirlingformel (1.4) auf 1 ?f folgende Glied plan 
auf die Näherung „,„;, möglichst ohne Einfluß bleibt. 


Maßgebend für die Bemessung der Schrittweite ist die im folgenden als „Schrittkennzalıl* 
bezeichnete Größe of 
k=|xh mit 

(vgl. 6. und 7.). Sie soll für mittlere Genauigkeitsansprüche den Wert von etwa k=0,05.... 0,10 
einhalten, wobei sich für die Anfangsiteration mit Rücksicht auf rasche Konvergenz eher der 
kleinere Wert empfiehlt. Bei veränderlichem # ist der Schritt h entsprechend abzustufen. 

Zur Schrittverdopplung brauchen die Werte y; und die mit doppelter Spanne h=2h 
versehenen Werte f;h der Zeilen » — 4, n—2 und » lediglich noch einmal angeschrieben 
und die f;h mit neuen Differenzen versehen zu werden, worauf die Rechnung ihren Fortgang 
nimmt. 

Zur Schritthalbierung ist eine neue Anfangsiteration erforderlich. Von der zuletzt be- 
rechneten »-ten Zeile aus kann man zwei Funktionswerte der Zeilen n —'/ und » +": 
zwischenschalten nach 


1 


Die mit halber Spanne h=h/2 versehenen Werte kA; i=n —1,n n, + '/ oder neu 
i=0,1,2,3) mit ihren Differenzen dienen von y, aus zur Berechnung von ,. Y,, 4, gemäß 
(1.6) mit anschließender iterativer Verbesserung der Werte. 

4. Ein Beispiel. 

Als Beispiel zur Anfangsiteration wählen wir die Differentialgleichung „=t-+y mit der 
Anfangsbedingung y—0 für = 0. Die Schrittgröße sei h=0,1, die Schrittkennzahl ist da- 
mit k= xh]—=0,l, also - namentlich für die Anfangsiteration — reichlich bemessen. Nach 


Tafel 3. Stufenweiser Aufbau dreier Rohwerte Y,- 


0,0 ° 0 

41) 0,1 0 0.010 0000 10 0000 
0.0 0 0 

4A 2) 0.1 0,005 000 0.010 5000 10 5000 

453) 0.2 0.021 000 0.022 1000 11 6000 1 1000 
0.0 0 
0.1 0.005 158 0.010 5158 105158 

14% 0.2 0,021 367 0.022 1367 116209 11051 
45) 03 0.049 800 0,034 9800 12 8433 12224 1173 


11) Die Differenzen sind hier, ähnlich wie im Zahlenbeispiel Tafel 3 und 4, nicht auf Lücke geschrieben, sondern 
je um eine halbe Zeile gesenkt. 


- 
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Aufbau der Rohwerte in Tafel 3 nach den Gleichungen unter 2 A erfolgt die iterative Ver- 
besserung nach (1.6) in Tafel 4. Bei Anwendung der Schablone Tafel 2 würden die einzelnen 


Tafel 4. Iteration (der Ausgangswerte Y,. Ys. nach Gl. (1.6). 


0,0 0 

r 0,1 0.005 170 0.010 5170 10 5170 

0.2 0.021 400 22 1400 116230 1 1060 
03 0.040 851 349851 12 8451 12283 1161 
00 

0.1 0.005 171 00105171 10 5171 

0,2 0.021408 22 1403 116232 11061 
0.3 0.040 858 319858 128455 1 2223 1162 
0.0 0 
0.1 0.005 171 00105171 105171 

0.2 0.021.103 22 1403 116232 11061 
03 0,040 850 349859 12 8456 12224 1169 


Näherungsstufen nicht untereinander geschrieben, sondern (mit Bleistift) einkorrigiert werden, 
bis die Werte feststehen. Bei halber Schrittweite h = 0,05 verläuft die Iteration entsprechend 
rascher. 


5. Glätten bei Kurvenaufrauhung. 


Indem die Formel (1.3) bei der Berechnung von „+, unter Umgehen von y,„ auf 
Ya, zurückgreift, besteht eine unmittelbare Bindung von einem „-Wert nur zum zweit- 
folgenden, während ein Zusammenhang zwischen benachbarten y-Werten nur mittelbar über 
die Differentialgleichung hergestellt wird. Das kann —- z. B. dureh 
gewisse Unebenheiten bei empirischen Funktionsgesetzen — dazu 
führen, daß sich im Verlaufe der Rechnung der Näherungsverlauf in 
zwei getrennte Kurvenzüge nach Art von Bild 1 aufspaltet. Die 
Aufrauhung der y-Werte pflanzt sich über die Differentialgleichung 
auf die Funktionswerte f und damit — verstärkt — auf die Diffe- 
renzen fort, bei denen sie sich in einem Schwanken etwa der | °f- 
Werte zuerst bemerkbar macht. 

In diesem Falle kann man die aufgespaltenen y-Kurven durch 
Einschalten einer Anfangsiteration wieder aneinander binden. Man 
rechnet zunächst, ausgehend von einem Wert 7,, welcher um drei 
Bild. Aufgerauhte „Kurve.  Sehritte hinter dem letzten bekannten Wert zurückliegt, 


1 


und erhält mit 22= y, — y, geglättete y-Werte auf einem mittleren (striehpunktierten) Kurven- 

vgl. Bild 1. Nach Verbessern von Funktionswerten und Differenzen erfolgt Anfangsiteration 
nach (1.6), bis die Werte feststehen. Die Rechnung nimmt ihren Fortgang, bis eine erneut 
auftretende Aufrauhung eine Glättung notwendig macht. 

6. Konvergenzbeschleunigung. 

Will man die Schrittweite % nieht so klein halten, daß die Iteration nach (1.5) von 


selbst genügend rasch konvergiert, so kann die Konvergenz durch eine besondere Rechnung 
beschleunigt werden. Man erhält für die Unterschiede 


der ». Näherungsstufe gegenüber der endgültigen Näherung durch Taylorentwieklung des 
Funktionswertes von fu +1 aus die Näherung 


+ l ) 3 h . . . . . . . . . . . 1.S). 
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Die Iteration nach (1.5) konvergiert also, solange k=|zh <3, und zwar um so rascher, je 


kleiner die Schrittkennzalhl % ist. 
Schreibt man nun (1.8) für » und » —- I an und subtrahiert, so erhält man nach kurzer 
Umreehnung 
l- „ah 
womit aus dem Unterschied zweier aufeinander folgender Näherungsstufen die Korrektur 
‚ und damit nach (1.7) die gesuchte Näherung 1 + ı berechnet werden kann. 

Für Systeme von Differentialgleichungen tritt an die Stelle der einen Gl. (1.9) ein 
Gleiehungssysten, welches hier nicht angeschrieben werden soll. Vielmehr sei auf den prak- 
tisch oft wiehtigeren Fall hingewiesen, daß nur eine einzige Veränderliche des Systems schlecht 
konvergiert. Um nun nicht wegen dieser Veränderlichen den Schritt verkleinern zu müssen, 
kann man für sie die oben beschriebene Konvergenzbesehleunigung nach (1.9) und (1.7) an- 
wenden. 

Für die Anfangsiteration nach (1.6) ergibt sich in ähnlicher Rechnung wie oben ein 
Gleiehungssystem für die Korrekturen —  - Für die verkürzte Anfangsiteration nach 
(1.6) innerhalb der gestrichelten Umrandung lauten «die Gleiehungen: 


2 
1 zahle’ = „NM 


7. Fehlerabschätzung zum Interpolationsverfahren nach (1.3). 
Für das Interpolationsverfahren nach Adams sind Fehlerabschätzungen schon wieder- 
holt durchgeführt worden. Auch für die Formel (1.3) erhält man wie bei G. Schulz*) für 
den Fehler = eine Schranke von der Gestalt 
wobei €, und €, von der Schrittanzahl » unabhängige Konstanten sind und z von der Schritt- 
kennzahl k=|xh| abhängt. An die Stelle der für das Adams-Verfahren beim Näherungs- 
grad r=3 (Berücksichtigung der Glieder bis [ ’f ein- 
schließlich) geltenden Formel rn 


—/S — 
24 (2? — 
tritt hier die Gleichung 
14: 
3(2?—1) 
(1.13). 


Bild 2 zeigt die Abhängigkeit der Größe 2 von k für die 
Adamssche Interpolationsformel und für Formel (1.3). Die 
für das Anwachsen der Fehlersehranke (1.10) entscheidende A 
Größe z fällt bei gleichem k für die Formel (1.3) kleiner, 
also günstiger aus. Für kleine Werte von k='zh!| gilt 
für (1.11) angenähert ” 


zzl-+I«h|, 


womit die Schrittkennzalil # einen gewissen Anhalt für das 


Anwachsen der Fehler gibt. Bild 2. Größe z in Abhängigkeit von «h. 


II. Verfahren für Differentialgleichungen 2. Ordnung. 


8. Interpolationsformel zum Störmer-Verfahren. 
Bei der Differentialgleichung 2. Ordnung 


gehorche f der Lipschitzbedingung 


Ns | | 
| 
/ | | | 
| | | 
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wobei M und Z, durch die größten Beträge der Ableitungen f, und f, innerhalb des in Betracht 
kommenden Bereiches ersetzbar sind. Aufbauend auf einer Reihe bereits bekannter Näherungs- 
werte und und damit bekannter Funktionswerte 7; = f(t;, y;, 9) = f; an gleichabständigen 
Stellen #; - 2,n —1,n) berechnen wir neue Näherungen y,„., und nach den 
Formeln 


Die zweite der beiden Gleichungen ist Formel (1.3), angewandt auf hy, die erste benutzt die 
Anfangsglieder der zweimal integrierten Stirlingformel '?) 
Die praktische Anordnung der Rechnung ist aus dem Schema Tafel 5 ersichtlich, in welchem 
die umrandeten Werte, versehen mit den Faktoren 1 und "fie bzw. 2 und "/s, benutzt werden. 


(2.4). 


Tafel 5. Rechenschema zur Interpolationsformel (2.3). 


Faktoren für y: 1 1/12 
t Faktoren für hy: 2 1/3 
In-ı In-ı l h 2}, 
In hin 
| In ri 


Zur Neurechnung der Näherungen y„;, und Ay„;, kann der noch unbekannte Wert 
I ®f,+, wieder aus der Reihe der schon bekannten | ?f-Werte durch Extrapolation geschätzt 
werden ®). Ausgehend von dem geschätzten Wert I ?f,”. , liefern die Gln. (2.3) die gesuchten 
Näherungen in iterativer Verbesserung nach 

\ 
] -(v-+- 1) _ 3 glr) (2.5). 
mit 

Für den in den Anwendungen wichtigen Sonderfall fehlender Abhängigkeit der Differential- 
gleichung von jy, also für 7 =f(f, y), brauchen die h y- Werte nicht berechnet, oder, wenn ihre 
Kenntnis erwünscht, nieht mit iteriert zu werden. 

Für die Bemessung der Schrittgröße h sind wiederum einerseits genügend rasche Kon- 
vergenz der Iterationen, andererseits die in den Gl. (2.4) und (1.4) auf I ?f folgenden Glieder 

50 und "fu +, maßgebend, welche auf die Näherungen und 

möglichst ohne Einfluß bleiben sollen. Über eine die Schrittweite bestimmende Schrittkenn- 
zahl # (ähnlich wie in 3) siehe in 13, letzter Absatz. 

9, Iterative Bestimmung der Ausgangswerte. 

Die Ermittlung der zum Beginn der Rechnung erforderlichen Ausgangswerte kann ähnlich 
wie unter I in stufenweisem Aufbau von Rohwerten mit anschließender iterativer Verbesserung 
erfolgen. Im Falle fehlender Abhängigkeit der Differentialgleiehung von ;,, also für y=fit; y), 
erübrigen sich Rohwerte und Iterationen für y. 


12) Vgl. G. Schulz: Formelsammlung 3), S. 88. 


13) Nimmt man anstatt 72 den Wert so hat man die nach abgebrochenen Extrapolat ions- 


formeln von Störmer ?, und Nyström ®). 


4 
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A. Aufbau dreier Rohwerte: 


1. Vorläufige Rohwerte auf Zeile 1: 
hıy)=hin + daraus: 
2. Verbesserte Rohwerte auf Zeile 1: 


daraus: 
3. Rohwerte auf Zeile 2: 


daraus: 


4, Verbesserte Rohwerte auf Zeile 1 und 2: 
1 
Un +h? m 57 2"), 


1 
daraus: 
5. Rohwerte auf Zeile 3: 


B. Iterative Verbesserung: 


hy | 


v v v v 1 


2.6). 


| 
| 
| 
| 
| 
3 l 
| | 
3 
2 | 
| 
N 
- 
ıl- 
re 
ler N 
+1 
In | 
Ing 
_ 
4 


50 Collatz u Zurmühl, Differentialgleichungen 1. und 2. Ordnung Fu aus 


Es wechseln ab: 
1. Verbessern der Werte und hy; nach 9; 
2. Verbessern der Funktionswerte flt;, und deren Differenzen, 
bis sich die Werte nicht mehr ändern. 
Für geringere Genauigkeitsansprüche und genügend kleine Schrittweite A" genügt 


wieder die Anfangsiteration mit zwei Zeilen 1 und 2, Tafel 6. Schablone zur Anfangsiteration nach 


also die Reehnung innerhalb der gestrichelten Um- ) nd zur fortlaufenden Rechnung der 
° y-Werte nach (2.3). 


randung. 

Die Anfangsrechnung kann wiederum durch An- y 
wendung von Schablonen erleichtert werden. Die T- TA Te: za 
Schablone zur Berechnung von y'') zeigt Tafel 6. Der [= 777 
untere Teil der Schablone dient auch zur fortlaufen- 
den Rechnung. Zur Berechnung der hy-Werte be- "7077: 7- 
nutzt man die in Tafel 2 für Differentialgleichungen [24 [ I pm: | 
1. Ordnung angegebene Schablone, wobei man y durch 


hy und hf durch h?’f zu ersetzen hat. 


10. Konvergenz der Iteration für die Ausgangswerte. 


Drückt man in den Formeln (2.6) alle Differenzen dureh die Funktionswerte f,. fi. fa fs 
aus, so erhält man Formeln der Gestalt 


3 
Mi 
wobei die o;; und r;, die Werte haben: 
360 60 120.45 s 24 24 24 
3 22 2 28 4 1 
5, W=elss 5° 
10 4020 ss 


Schreibt man (2.7) für » und » - I an und zieht die entstehenden Gleichungen voneinander 


ab, so erhält man 
3 


=] 


9; 


Die Differenzen der Funktionswerte auf der rechten Seite werden mit Hilfe der Lipschitz- 
bedingung nach (2.2) abgeschätzt: 


(28). 


«(v-+1]) Y (—]) «(—]) 


Mit der Abkürzung 
= 


ergibt sich nach ee der beiden Gln. (2.5) mit M bzw. L und Addition die Beziehung 


3 
k=1 


tz- 


Ing 
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Mit einem Gleichheitszeichen werden daraus Gleichungen für obere Schranken V, der 


Größen v,''. Diese Gleichungen lassen sich bequem in Matrizenform schreiben. Dazu faßt 


man die V‘' als Komponenten eines Vektors dv" auf. Mit der Matrix 


lautet dann Gl. (2.9): 
A 


oder, indem man diese Beziehung rekursiv anwendet, 


= 
Nun konvergiert I’ ve +" —(_V’W)v“” nach bekannten Sätzen, wenn die Matrix A charakteri- 
v—( 
stische Zahlen hat, deren sämtliche Beträge kleiner sind als 1. Ist das erfüllt, so konvergiert 
also auch die Anfangsiteration nach Gl. (2.6). Zur Untersuchung der charakteristischen Zahlen 
der Matrix A schreiben wir 


unter Einführen eines Parameters p. Die größte charakteristische Zahl s einer zweiten Matrix 


| £ 


7 ist dann eine Funktion von p, welche 
\ | punktweise berechenbar und in Bild 3 
\ dargestellt ist. Da der Betrag der cha- 
rakteristischen Zahl von W kleiner als 
\ 1 sein soll, hat man die Bedingung 


s-Mh’<1. 

> Als Konvergenzgrenze erhält man damit 

für die beiden Größen Mh? und Lh die 

Beziehungen 


1 
2 1: 
—s(p)' . . (2.13), 


0 7 2 
welche in Parameterform die fast gerad- 
linige Grenze des Konvergenzbereiches 
in einer Ebene mit den Koordinaten Ih 
und Mh? angeben, siehe Bild 4. 


Bild 3. (oben). Verlauf s = s(p). 


Bild 4 (rechts). Konvergenzgrenzen 
fürdie Anfangsiteration nach GI. (2.6). 


0 7 


11. Fehler der Ausgangswerte. 
Für die zweite Ableitung j der exakten Lösung y(f) schreiben wir die Newtonsche 
Interpolationsformel mit absteigenden Differenzen '') 


mit dem Restglie«dl — Dia — 2) (u 


3) afıa) 
(e), 


R= 


1 
wo statt die Veränderliche — (t -#,) eingeführt und eine Zwischenstelle aus dem 
Intervall ist. 


Eine Integration liefert 


2 3 a a 
t+tuh u 
mit S-\ Rdt=h\Rdu. Das Integral kann abgeschätzt werden durch 


u 
5 


14) (1, Schulz: Formelsammlung 3), S. 74. 


h 
r 
2 | | 
| 
F | | 
| | | | \ 
er 
— 
8). 
ıaAA 
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Die Ausrechnung ergibt 
19 
Imax 


a 
Seh f max 


h? 


Eine nochmalige Integration liefert 
3 u* u? 


tuh 
mit 
to 
Bei Ausrechnung erhält man die Werte 


- 4! 50: 

Die Näherungen %,, 9, genügen nun den Gleichungen, die aus (2.6) durch 
Fortlassen der Klammern (»v) und (» +1) hervorgehen. Zieht man diese von -(2.14) bzw..(2.15) 


ab, so erhält man für die Fehler «, bzw. £,: 
3 


3 
k=1 
Schätzt man die Differenzen der Funktionswerte mittels der Lipschitzkonstanten nach 
(2.2) ab, so folgt 


3 
N + |Sil. 


Mit Gleichheitszeichen ergeben sich 6 Gleichungen für 6 unbekannte obere Schranken Y,; 
und Z; der Fehler und nämlich: 
3 
3 (i = 1, 2, 3) 


Faßt man die Unbekannten als Komponenten eines Vektors x auf 
so läßt sich das Gleichungssystem (2.16) in der Gestalt schreiben 
wobei die Matrix 


und der Vektor 
(D,, U,, U,, S,, S,) (r,, 
Die Lösung von (2.17) lautet in Matrizenform mit der Einheitsmatrix €: 


ist. 


m 
| 
30° 
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Haben alle Elemente von A höchstens den Betrag a, so sind alle Elemente von A? dem 
Betrage nach höchstens’ gleich 6a? und, die von A“ höchstens gleich 46a, Die Auf- 


summierung der Reihe (2.18) liefert dann, wenn man die obengenannte Abschätzung durch- 
führt, für die oberen Schranken Y; und 2, der Feilkeleäge: 


(2.19). 


12. Konvergenz der Iteration nach Gl. (2.3). 
Schreibt man die Gl. (2.3) in der Gestalt 


h 


für die Indizes » und » - I an und zieht die so entstehenden Gleichungen voneinander ab, 
so erhält man 


+1) (r) (—1) 
— Ya +ıT fazı) 


’+)) ) h 


Schätzt man die rechten Seiten mit Hilfe der Lipschitzkonstanten ab, so erhält man 


für die Änderungsbeträge 


_ 
die Ungleichungen 
h? 
v) v—| 


h 


i> 


Mit der Abkürzung folgt 


+1 <( 


Da die nicht sind, nehmen sie für 


mindestens so stark wie die Glieder einer geometrischen Reihe ab. Ist 
diese Bedingung erfüllt, so ist die Konvergenz der Iterationen gesichert. 
u} In Bild 5 ist in einer Ebene mit den Koordinaten Lh und Mh? das 
ne Konvergenzgebiet eingezeichnet. Der Vergleich zeigt, daß das Gebiet viel 
Bild 5. Konver- größer ist als das Konvergenzgebiet für die Anfangsiteration. Man führt 
genzgreuze fürdie daher die Anfangsiteration vorteilhaft oft mit kleinerer Schrittweite als die 


Iterationen nach 
G1. (25). eigentliche Rechnung durch. 


13. Fehlerabschätzung für die Formel (2.3). 
Den Ausgangspunkt für die Fehlerabschätzung bildet wieder die Newtonsche Inter- 
polationsformel für aber diesmal mit aufsteigenden Differenzen geschrieben (mit !=t,+,+ uh). 


1 1 2 


mit dem Newtonschen Restglied 


+1 2) 3 
ur 


wo r eine Zwischenstelle aus dem Intervall tu+, ist. 


ha 
4 
1 
MR? 
5 
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Eine Integration dieser Gleichung liefert 


t u 
mit S=h\Rdu. Dieses Integral kann abgeschätzt werden dureh 
; 
wobei Tin) = \ du gesetzt ist. 
—1 
Es ist T(-2)= du=1. Insbesondere folgt für oder —2: 
mit 
S, > max * 
Weitere Integrationen liefern aus (2.21) 
1 
= \ ydt+ \ yıat + U, | 
n 27 
mit (2.2). 


\Sdu+h \Sdu > Imax 


-1 


Die Näherungen y,„, 7. genügen nun den Gleichungen 


72 


1 10 


Subtraktion von (2.21) bzw. (2.22) liefert für die Fehler e, und £,: 


2 
Sntı Belfttun+ 1-0) — 
(= 


oder wenn die Differenzen der Funktionswerte mittels der Lipschitzkonstanten abgeschätzt 
werden 


U, 


.n-ı sh Be(M \enzı-el + 


YıJM 
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Bedeuten Y, und Z,„ obere Schranken für und '{„|, so kann man für Y,„ und die 
Gleichungen verwenden 


V° Da, (M 1-2 L -o) U, | 
(2.23). 
Mit diesen Formeln ist im Prinzip die Fehlerabschätzung geleistet. Hat man die Fehler- 
schranken Y, Z bis Y„, Z,„ berechnet, so liefern diese Formeln (2.23) zwei Gleichungen zur 
Berechnung von und Zu+:- 


Diese schrittweise Berechnung der Fehler ist oft etwas mühsam, und man erhält keinen 
Ueberblick über das Wachstum der Fehler. Man kann die Größen Y, und Z, durch ein 
Potenzgesetz abschätzen, da (2.23) lineare Differenzengleichungen für Y,, Z, sind. 


Eine Sonderlösung der inhomogenen Gleichungen ist 
'S, U, S, U, 


Die Lösungen der den (2.23) entsprechenden homogenen Gleichungen werden durch den Ansatz 
gewonnen 


Dabei ist z Wurzel der Gleichung 


2 —22+1—- Mh — Lh? N 
D 


Diese Gleichung hat mehrere Wurzeln z. 
Entscheidend für das Fehlerwachstum ist die 
Wurzel z vom größten Betrage. In Bild 6 
ist diese Wurzel in Abhängigkeit von den 
Größen Lh und Mh? aufgetragen; dabei sind 
die Kurven z2=konst. gerade Linien, denn 


aus Gl. (2.24) folgt: 


Bild 6. Geraden <—=konst. in der Lh- M h?-Ebene. 


0==0 


wobei sich das Glied A? LM heraushebt und eine in Ah und Mh? lineare Gleichung entsteht. 


Zur praktischen Fehlerabschätzung wird es meist vorzuziehen sein, die Fehler schritt- 
weise aus (2.23) zu berechnen, da bei der Abschätzung der Fehler nach einem Potenzgesetz 
die Fehlerschranken, besonders bei großer Schrittanzahl », die wirklichen Fehler in vielen 
Fällen oft stark überschreiten. 


Um über ein längeres Intervall, in welchem sich die Werte Lh und Mh? ändern, an- 
nähernd mit gleicher Genauigkeit, d. h. mit gleicher Stärke der Fortpflanzung der Fehler zu 
rechnen, kann man sinnvoll die Schrittweite so abstufen, daß man sich in Bild 6 annähernd 
auf der gleichen Fehlergeraden z = konst. bewegt, wobei der Zahlenwert von z sich nach den 
Genauigkeitsansprüchen richtet. Der Wert von z— 1 entspricht dabei wieder angenähert der 
Schrittkennzahl k=|zh|, welche in Nr. 3 für Differentialgleichungen 1. Ordnung eingeführt 
und für mittlere Verhältnisse mit 0,05 bis 0,10 angegeben worden war. Für kleine Werte 
von z—1 gilt die Näherung: 
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Bemerkung zur Kirchhoffschen Platten- 
strömung. Es ist bekannt. daß eine Unstetig- 
keitsfläche, wie sie z. B. bei der bekannten Kirch- 
hoffschen Plattenströmung auftritt, unter Berück- 
sichtigung der Flüssigkeitsreibung nur als Anfangs- 
zustand möglich ist. und im weiteren Verlauf rasch 
in der Weise vernichtet wird. daß sich die Wirbel. 
welche auf der Trennungsfläche liegen, in die Flüs- 
siekeit hinein verbreiten: die Trennungsfläche er- 
weist sich also als stark labilt). Man sieht, daß 
einer solchen Betrachtung die folgende Voraus- 
setzung (welche auch im Wesen der Flüssigkeits- 
reibung liegt) zugrunde liegt: 


(V): Die beiden vorher getrennten 
Flüssigkeiten affizieren sich über 
die Trennungsfläche hinweg. 


Für wirkliche Flüssiekeiten ist diese Voraus- 
setzung stets erfüllt. da sie mit der Brownschen 
Molekularbewegung ausgestattet sind. Für ideale 
jedoch nieht. Will man also hier in Übereinstim- 
mung mit der Theorie die scharfe Trennungsfläche 
behalten, so hat man dafür künstlich Sorge zu tra- 
een, indem man z. B. die Trennungsfläche als 
Membran gestaltet, oder die ruhende Flüssigkeit 
dureh ein anderes Medium wie Luft ersetzt und so 
die Gefahr einer Affizierung nahezu ausschaltet 
freier Helmholtzscher Strahl); in jedem Fall 
erhält man dann eine ideale Strömung, die physi- 
kalisch. d. h. im Rahmen der Theorie der idealen 
Flüssigkeiten. möglich ist: doch scheint sich, so- 
lange man in diesem Rahmen bleibt, die Instabili- 
tät der Trennungs- bzw. Begrenzungstläche nicht 
zu ergeben. Indessen können wir (ähnlich wie dies 
Prandtl in anderer Hinsicht bei Aufstellung 
seiner Tragflügeltheorie getan hat) durch einfache 
Hinzunahme der Voraussetzung (N) diesen Rahmen 
erweitern (und so in gewisser Weise der Reibung 
Rechnung tragen, ohne daß es jedoch nötig wäre. 
von der Reibung zu sprechen!) und uns fragen, ob 
wir vielleicht in der solehermaßen ergänzten Theo- 
rie der idealen Flüssigkeiten besagte Instabilität 
einsehen können. Die folgende heuristische Be- 
trachtung will versuchen, auf diese Frage eine be- 
jahende Antwort zu geben. 


dn__ 


Bild 1. 


Betrachten wir das in Bild 1 abgegrenzte Flüssig- 
keitselement: solange die Voraussetzung (V) nicht 
zutrifft, zerfällt es in zwei Teilchen. ein ruhendes 
und ein mit der Geschwindigkeit bewegtes, (die 
nichts miteinander zu tun haben: ihre Begrenzungs- 
linien parallel zur Trennungslinie sind mit 0, 1. 2 
bezeichnet. die diesbezüglichen Drucke seien Po. 
Pi. Pa. wo also P,—=P, ist: ds sei ihre Länge 
parallel der Trennungslinie, dr ihre Höhe. In der 


1) Vgl. Cl. Schaefer: Einführung in die Theoretische 
Physik, 1. Bd.3. Aufl. Berlin u. Leipzig 1929, S. 937 bis 943. 


Normalenrichtung haben wir dann für das bewegte 
Teilchen mit r als Krümmungsradius die Gleichung 


(Ps 


und für das ruhende Teilchen die Identität 


Aus diesen beiden Gleichungen ersieht man, dab 
beim Durchtritt durch die Trennungslinie außer der 
Gieschwindigkeit auch das normale Druckgefälle 
einen Sprung erleidet. 

Lassen wir jetzt die Voraussetzung (V) zu, so 
bilden die beiden Flüssigkeitsteilchen von vorhin 
ein Ganzes, welches die mittlere Geschwindigkeit 5 
hat. An Stelle von (1) erhalten wir dann 


(U/2)? 
r 


— p) ds . (3). 


Weren D#, = P, ist hier gemäß (1) die linke Seite 
aber doppelt so groß wie die rechte! Ohne eine 
entsprechende Fremd- oder Haltekraft von der 


Größe — (92 — ds könnte also die neue 


Strömung nicht bestehen, d. h. die Instabilität der 
Trennungsfläche wäre erwiesen. Genauer ergäbe 
sich: dem zu großen Druckgefälle kann das Teil- 
chen offenbar nur dadurch entgehen, daß es seinen 
Krümmungsradius entsprechend verkleinert. Auch 
sieht man, daß dieses Resultat nicht an die Beson- 
derheit der Kirchhoffschen Plattenströmung 
gebunden ist. und z. B. auch für Trennungsflächen 
gilt. welche nicht ein Totwasser abgrenzen. 

In dieser Betrachtungsweise ist eine Strömung 
mit Totwasser nur stabil, wenn die Trennungslinie 
an jeder Stelle einen unendlich großen Krümmungs- 
radius besitzt, also geradlinig ist. Aber auch dieser 
Fall ist i.a. nicht möglich: ist nämlich 2 (z) die kom- 


plexe Stromfunktion, also F (2) — die kom- 
dz 


plexe Geschwindigkeit, so würde F(z) an der Be- 
grenzung eine Konstante € und daher überall = € 
sein, d. h. es läge eine reine Parallelströmung vor. 


Bild 2. 


Bild >. 


(43083) 


2 veranschaulicht einen solchen Fall. Liegt 
jedoch eine Strömung geren eine Schneide wie in 
Bild 3 vor, so kann die Trennungslinie nicht gerad- 
linie sein, ist also (nach unten) gekrümmt: alsdann 
wird sie an jeder Stelle ihren Krümmungsradius 
dauernd zu verkleinern suchen, und sich daher in 
bekannter Weise zu einem Wirbel aufspulen. 
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Kurze Bemerkung zur graphischen 
Lösung gewöhnlicher linearer Differen- 
tialgleichungen 1. Ordnung. 1)ie Praxis ver- 
langt gelegentlich die Kenntnis einer Funktion, die 
einer linearen Differentialgleichung 1. Ordnung ge- 
nügt. Zwar gehört die Auflösung einer solchen 
Differentialgleichung längst zum Allgemeingut nicht 
nur des Mathematikers, sondern auch des Techni- 
kers, es sei aber doch gestattet, auf die folgende 
graphische Lösungsmethode hinzuweisen. Sie ist 
besonders dann geeignet, wenn es darauf ankommt. 
für beliebige Stellen © die Funktion y(r) nume- 
risch zu kennen, falls beliebige Anfangswerte u 

für beliebige Stellen = .r, vorgeschrieben sind. 
Auf ein solches Problem führt z. B. die Ermittlung 
der Bahngeschwindigkeit V «eines Flugzeugs als 
Funktion der Flughöhe h, das unter konstantem 
Gleitwinkel (speziell im senkrechten Sturzflug) aus 
einer vorgegebenen Ausgangshöhe h — h, mit be- 
kannter Anfangsgeschwindigkeit 7 — V, unter Be- 
rücksichtigurg der mit der Höhe variablen Luft- 
(lichte eine gewisse Höhendifferenz durchstürzt !). 

Sei die lineare Differentialgleichung 1. Ordnung 


mit im Intervall a&r= 5 stetigen Koeftizienten 
und f, (x) vorgelegt. Man setze 


x 
(t) dt 
ea 


— 


und — was wegen F, (#)=#=0 
möglich ist — . 
hit) 


a 


di=F,(z) 


Dann ist die aus der Gleichung 


Y 
F,(z) 


Y (Ko ) 


F, (+) = F, (z.) 


— (8) 
zu ermittelnde Funktion die eindeutig be- 
stimmte Lösung der Differentialgleichung (1), die 
für — den Wert y= Yo = Y (X,) annimmt. Mit 


vi) „ 
=-2 


stellt aber die aus (3) hervorgehende Gleichung 
.. (4 


eine Geradenschar im reehtwinkligen 3-System 
dar: Parameter ist dabei :r. 

Man denke sich die zu gewissen Parameterwerten 
zehörigen Individuen dieser Geradensehar gezeich- 
net. Wünscht man nun den Wert v=y, der Lö- 
sune Y von (1) für .r, zu kennen, falls 
(ler Anfangswert 4 = für — .r, vorgeschrieben 
ist, so gehe man vom Punkt y, der Abszissenachse 
(y-Achse) senkrecht hoch bis zum Punkt P, auf der 
Geraden mit dem Parameterwert ,, dann von P, 
horizontal weiter bis zum Punkt P, auf der Ge- 
raden mit dem Parameterwert 7, und hat in der 
Abszisse von P, den gesuchten Wert y, gefunden 
(siehe Bild 1). 

Um nieht nach dem Parameter 7 interpolieren zu 
müssen und um dureh das Zeichnen von vielen Ge- 
raden die Übersichtlichkeit nicht leiden zu lassen. 
ist es vorteilhaft, das v,2-Achsenkreuz überhaupt 


1) Siehe hierzu H.Behrbohn® Graphische Konstruktion 
von Sturzfluggeschwindigkeiten. Luftf.-Forschg. Bd. 17 
(1940) Lfg. 6, S. 167 bis 171, wo für den genannten speziellen 
Fall aus der Praxis die im folgenden skizzierte Methode 
mit Erfolg verwendet wurde. 


\ 


Bild 1 


mit keiner fertig gezeichneten Schargeraden zu be- 
lasten, sondern dafür lieber die jeweils gebrauch- 
ten Geraden ad hoc rasch zu konstrnieren. 

Ein Mittel hierzu ist beispielsweise folgendes: Die 
(reraden (4) schneiden auf der y-Achse den Ab- 
schnitt Y(z)=F,(2)F, (x), auf der z-Achse den 
Achsenabschnitt Z (2) =F, (“) ab. Es genügt dann 
vollkommen. sowohl der z-Achse als auch der Y- 
Achse jeweils eine z-Achse zu überlagern, über 
der die entsprechenden Achsenabschnittskurven Y 

und Z— aufgetragen werden (siehe 
Bild 2). Statt der Abschnitte auf der v- und z- 
Achse können natürlich auch irgendwelche andere 
für das jeweils vorliegende Problem evtl. besser 
zeeignete Bestimmungsstücke der Schargeraden als 
Funktion des Scharparameters angegeben werden. 


Von Nutzen kann auch die möglicherweise vor- 
handene Einhüllende E der Geradenschar (4) sein. 
Durch Differentiation nach dem Scharparameter 
entsteht aus (4) die weitere Gleichung 


1F,(x) 
2 


so daß sich E unter Beachtung von (2) in der Pa- 
rameterdarstellung 

(2) 


2=Fı(2)+ 


| 
Fol) 


; | 57 
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ergibt; Y und = sind die Koordinaten des Berüh- 
rungspunktes B von E mit «der zum entsprechen- 
(len Parameterwert © gehörigen Schargeraden. Liegt 
E im 4.3-System gezeichnet vor, so hat man eine 
eute Kontrolle für das Zeichnen der Geraden. da 
diese ja E in B berühren müssen. 
Des weiteren liefert die Kurve E in der Abszisse 
Y/g ihres Schnittpunktes mit der Geraden 2 vonst 
c einen analytischen Extremwert derjenigen Inte- 
eralkurve 4 von (1). die durch die An- 
fangsbedingung = Yo gerade den durch 
F 
lest; aus 
0, 4. R. 


bestimmten fest- 


= € 


In (x) 
nach (1) für Y Up. 


folet nämlich fo 

Die oben erwähnte Aufgabe aus der Flugmechanik 
führt zur Differentialgleichung 


3 
y=-ay+ 


falls mit. die Luftdichte, mit y das Quadrat der 
Flugseschwindigkeit und mit «@ und ‚geeignete 
Konstante (in die der aerodynamische Widerstand 
und die Flächenbelastung des Flugzeugs eingehen) 
bezeichnet werden. Als funktionale Bindung zwi- 
schen Flughöhe und Luftdichte ist dabei ein gegen- 
über «der polytropischen Atmosphäre angenäherter 
Exponentialansatz getroffen worden. 


Für dieses Beispiel wird (mit @ — 0) 


=er®%, 


(in der Bezeichnungsweise von Jahnke-Emde, Funk- 
tionentafeln. 2. Aufl.. Berlin 1933. S. 79). 


Die Geradenschar (4) lautet demnach 
Sie hat die Abschnittsfunktionen 
\ 
Z(r)—=BEi 


und die Einhüllende E 


z=ßEi(l(ar) — 


as 


Zu einer Anfangsluftdichte 7 — 7,. die zu einer 
bekannten Gleitflugausgangshöhe gehört. und zu 
einem vorgegebenen Geschwindigrkeitsquadrat u 

Yo im Anfangszustand gehört dann jeweils ein 
eindeutig bestimmtes Geschwindigkeitsquadrat v7), 
das mit abnehmender Flughöhe (d. h. mit zuneh- 
mender Luftdichte «*) zunimmt, bis ein durch 

zegebenes Maximum erreicht wird; an- 
schließend fällt annähernd linear mit 


Bild 2 gibt die Verhältnisse dieses speziellen Bei- 
spiels im Intervall 1&.r=4 wieder. Dabei ist die 
Achsenabschnittsfunktion Zr) ersetzt durch die 
Abschnittsfunktion Z = 2 (4 = 0,2) der Schar- 
geraden mit der zur 2-Achse parallelen Geraden 
v—02. Für 2 — 1,0 sei v= 0,41. Bild 2 gibt dann 
beispielsweise für 7 — 2,0 den Wert u — 0.45. 
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Zur Theorie der Wirbeilbewegung. 

In einem Aufsatz mit derselben Überschrift ') 
hat A. Kneschke ein Minimalprinzip angege- 
ben, aus dem die Bewegungsgleichungen isolierter 
Potentialwirbel in ebenen, einfach zusammenhän- 
genden Gebieten einer idealen Flüssigkeit ab- 
geleitet werden können. Dabei geht Kneschke so 
vor, daß er zunächst diese Bewegungsgleichungen 
potentialtheoretisch berechnet und dann ein Vari- 
ationsprinzip konstruiert, das mit diesen Bewe- 
eungsgleichungen äquivalent ist. Im folgenden 
soll nun gezeigt werden, wie man dieses Minimal- 
prinzip mit Hilfe des Hamiltonschen Prinzips 
direkt aufstellen und damit auch physikalisch 
deuten kann. Ferner wird dieses Prinzip auf das 
Wirbelband einer tragenden Linie angewandt, wo- 
bei sich die Prandtische Gleichung für den Ab- 
wind ergibt. 


1. Isolierte Wirbelfäden. 


Wenn sich in einem unendlich ausgedehnten ebe- 
nen Flüssigkeitsgebiet isolierte Wirbel befinden, 
so ist die Eigengeschwindiekeit 1; 
des ö-ten Wirbels an der Stelle z=ir; 
vleich der von den anderen Wirbeln 7, induzier- 
ten Geschwindigkeit d; = mit 
ist die Geschwindigkeit. die der Wirbel 7, an der 
Stelle x; induziert: sie steht senkrecht zur Ver- 
bindungslinie r;, der beiden Wirbel 7‘, und 
Denkt man sieh den Wirbel 7’, dureh einen ruhen- 
den starren Körper in einer Strömung mit der 
Geschwindigkeit d; und der Zirkulation 7, ersetzt, 
so wird auf diesen Körper von der Flüssigkeit 
nach dem Kutta-Joukowskischen Satz eine Kraft 
senkrecht zu db; ausgeübt. In diesem Sinne kann 
man daher sagen, daß ein Wirbel 7, auf einen an- 
deren 7; eine Kraft in Richtung der Verbindungs- 
linie der beiden Wirbel ausübt, die folgende 
Größe hat: 


1 
ik 


[din = dir 


File, 


wobei |dir dir. 
—o 


und IK 
o ist die Diehte der Flüssigkeit. Zwei Wirbel mit 
eleichem Drehungssinn (7, und 7‘, gleiches Vor- 
zeichen) stoßen sich ab, St; und r;,,, haben dieselbe 
Richtung (s. Bild 1): bei entgegengesetztem Dre- 
hungssinn ziehen sie sich an, die Richtung von $;, 
ist entgegengesetzt der von Tj;. 


Bild 1. 


Mit Hilfe der so eingeführten Kräfte kann man 
nun die kinetische Energie E des durch die Wirbel 
verursachten Strömungsfeldes leicht berechnen. 
Zunächst ist es naheliegend, sich alle Wirbel im 
Unendlichen zu denken und nun jeden einzelnen 
Wirbel an die vorgeschriebene Stelle r,; zu brin- 
ven. Um die Arbeit. die dabei infolge der gegen- 
seitigen Abstoßung bzy. Anziehung der Wirbel zu 
leisten ist, ist dann % kinetische Energie des 


1) A. Knesehke: Zur Theorie der Wirbelbewegung. 
Z. angew. Math. Mech. Bd. 18 (1938) S. 343. 
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Feldes vergrößert worden. Bei dieser Rechnung 
muß man aber mit Grenzwerten arbeiten, die 
‚dureh die folgende Betrachtung umgangen werden 
können. Es sei E die kinetische Energie der 
Flüssigkeit. wenn sieh die Wirbel 7’, an den Stel- 
len tr; befinden und die gegenseitigen Abstände T;,, 
haben. Wir denken uns nun, ohn® die Zirkulation 
zu Ändern. einige Wirbel verschoben, so dab die 
Abstände vj, und die kinetische Energie E* be- 
tragen. Die Differenz E*—E entspricht der beim 
Versehieben der Wirbel von einer äußeren Kraft 
zu leistenden Arbeit 4A: 


* 
Un 
A=E*—-E= \ ik [si o 
Un 


ler — 5, 


Daraus folgt nun: 


E=,_ + Eu 
ET 

wobei 


Da bereits di» Energie eines Potentialwirbeis un- 
endlich groß ist, sind Ex und damit auch E eben- 
talls unendlich groß. Jedoch hängt Ex nicht von 
der Lage der Wirbel ab, sondern ist eine orts- 
unabhängige Konstante. Die Abhängigkeit der 
kinetischen Energie von der gegenseitigen Lage 
der Wirbel, d. h. also von den r;,. wird somit 
ganz durch das erste Glied erfaßt. Dieses ist 
sleiehlautend mit der Kirchhoffschen Bahnfunk- 
tion, die somit. wie bereits Kneschke erwähnt, bis 
auf eine additive Konstante gleich der kinetischen 
Energie des Strömungsfeldes ist. 


Um das Hamiltonsche Prinzip anwenden zu kön- 
nen, müssen wir nun noch das Potential der ge- 
venseitigen Kräfte der Wirbel aufstellen. Dabei 
ist zu beachten, daß die Kraft auf einen Wirbel 
nicht der Beschleunigung. sondern seiner Ge- 
schwindigkeit proportional ist. Man kann aber 
auch in derartigen Fällen das erweiterte Hamil- 
tonsche Prinzip anwenden ?). wenn man das 
Kräftepotential U so aufstellt. daß 


d 
 dt\dy; 


Y;= 


daher setzen: 


Dann lautet das Hamiltonsche Prinzip, angewandt 
auf die Wirbelbewegung: 


mub 


(3). 


2) Handbuch der Physik V, Berlin 1927. 
») 8; ist die Resultierende der Kräfte, die die Wirbel 
7; auf den Wirbel 7; ausüben. 


E)dt | 
to 


— 


ET 

Das ist aber das von Kneschke angegebene» Mini- 
malprinzip, das man jetzt so formulieren kann: 
Für die wirkliche Bewegung der freien Wirbel ist 
der zeitliche Mittelwert der Differenz aus dem 
(oben definierten) Potential der gegenseitigen An- 
ziehungskräfte der Wirbel und der kinetischen 
Energie des Strömungsfeldes stationär. 


2. Der Abwind hinter einer tragen- 
den Linie. 

Um das obige Minimalprinzip auf die Tragtlügel- 
theorie anwenden zu können, müssen wir es zuU- 
nächst auf die Bewegung einer Wirbelschieht über- 
tragen. Die Koordinaten der Schicht seien 7 und 

die Zirkulation y (7.{) und s die Bogenlänge 
(s. Bild 2). Für ein infinitesimales Stück der 


Bild 2. 


Wirbelschieht =y (y.2) ds an der Stelle 
benutzen wir nun die Gleichungen (2). (3) und (4). 
wobei durch d 7’ usw. zu ersetzen ist. Dann er- 
gibt sich: 


1 


(s) 


stearfı 


ds | ar 


Dabei ist dZ/’ konstant zu halten. Die zu dieser 


Variationsaufeabe 


ö\Fy,2; y,2)dt=0 


zugehörigen Euler-Gleichungen 


d 


dt und 


Fy Fj 


lauten: 


= \ ds 
(s) 

1 \ — 2) 
2) 


s) 


ds 


Aus diesen Gleichungen läßt sich jetzt als Son- 
derfali der Abwind hinter einer tragenden Linie 
berechnen. Weit hinter einer tragenden Linie 
können wir nämlich die Strömung in einer Quer- 
ebene senkrecht zur Flugriehtung als zweidimen- 
sional betrachten. Das abgehende Wirbelband soll 
noch eben sein: Z=0, 20, d. h. wir setzen vor- 
aus, daß die Eigenbewegung der freien Wirbel- 
fäden klein ist gegenüber der Fluggeschwindig- 
keit. Wie üblich sei 


59 
tı 
| 
(4). 
2n 
| | | 
= 79:5) al ds 
18). 
= 
mit 
0 
U=, | 
> 
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8 
| dy 


(lie Zirkulation der Wirbelschieht an der Stelle 
vn. und b die Spannweite. Wir betrachten nun 
an irgendeiner Stelle y des Wirbelbandes den zu- 
e N dy und fragen 
nach seinen Bewegungsgleichungen. Für diesen 
Fall ergeben die Gleichungen (6): 


dy | 


Die zweite Gleichung sagt aus. dab die Wirbel- 
fäden keine Querbewegung in Y-Richtung aus- 
tühren, die erste gibt die Abwärtsgeschwindigkeit. 
Berücksichtigt man, daß der Abwind w (y) an der 
tragenden Linie selbst. wo die freien Wirbel nur 
in einer Richtung ins Unendliche verlaufen. rur 
halb so groß ist wie der weit hinten. so findet 
man die bekannte Formel für den Abwind am Ort 
der tragenden Linie: 


gehörigen Wirbelfaden = 


+ 


1 
on 


— 


y=0 


Damit ist das Variationsprinzip (5) bestätigt. 


Ein anderes Minimalprinzip für die tragende 
Linie hat F. Ziller*) von der Abwindgleichung 
ausgehend, aufgestellt. Dabei wird aber nicht, wie 
hier, die Bahn eines abgehenden Wirbelfadens 
selbst variiert, sondern der effektive Anstellwinkel 
an einer Stelle des Wirbelbandes. Ziller berech- 
net damit Auftriebsverteilungen, indem er die 
Variationsaufgabe direkt mit einem Ritz-Ansatz 
löst. Physikalisch erscheint es nun wahrschein- 
lich. daß sein Minimalprinzip auch mit dem Ha- 
miltonschen Prinzip in Beziehung steht: es ist je- 
doch zu undurchsichtig. um diesen Zusammenhang 
erkennen zu lassen. Andererseits liegt mathema- 
tisch kein zwingender Grund für einen derartigen 
Zusammenhang vor, da es durchaus möglich ist. 
daß zwei verschiedene Variationsaufgaben zu ein 
und derselben Euler-Gleichung führen, ohne daß 
zwischen den variierten Ausdrücken selbst eine 
innere Verwandtschaft besteht. 

Göttingen. K. Wieghardt. 312 


4) F. Ziller: Be 
endlicher Spannweite. 


Beitrag zur Theorie des Tragflügels von 
Ing.-Arch. Bd. 11 (1940) S. 239. 


BUCHBESPRECHUNGEN 


Prof. Dr. HUGO DINGLER, Max Planck und 
die Begründung dersogenannten mo- 
dernen theoretischen Physik 32 8. 


Berlin 1939, Ahnenerbe-Stiftung-Verlag. Preis 
0,50 M. 


Das Heft enthält zwei bereits in einer Fachzeit- 
schrift veröffentlichte Aufsätze von Dingler: 1. „Zur 
Entstehung der sogenannten modernen theoreti- 
schen Physik“ und 2. „Determinismus oder In- 
determinismus?* In beiden wendet sich der Verf. 
zeren die erkenntnistheoretischen und methodo- 
iorischen Anschauungen Plancks. So bekämpft er 
in 1. vornehmlich das von Planck der Physik ge- 
setzte Ziel: die Loslösung des physikalischen 
Weltbildes von der Individualität des bildenden 
Geistes, um allmählich zur Erkenntnis der realen. 
dd. h. von den menschlichen Sinnen unabhängigen 
Außenwelt zu gelangen. Dingler stellt dieser „ge- 
spenstischen metaphysischen Überwelt“ die „pri- 
märe Lebenswelt“ «egenüber, die sich ihm ergibt 
aus „jenem frischen, noch von keinen Zweifeln an- 
zekränkelten Erleben des Alltags. das wir alle 
kennen und haben. und in dem auch die Sinnes- 
wahrnehmungen nicht nur subjektiv sind, sondern 
jeder unmittelbar weiß, daß sein Nebenmensch im 
wesentlichen genau die gleichen hat“. Es geht 
also wieder einmal um das Wirklichkeitsproblem, 
‚las die Philosophen seit Jahrtausenden beschäftigt. 
ohne daß bisher eine zweifelsfreie Lösung geglückt 
wäre. — Auf ein Mißverständnis auf S. 7 sei hin- 
zewiesen. Dingler beziehtigt da Planck einer An- 
schauung, die dieser indessen als Folgerung des 
von ihm bekämpften Machschen Sensualismus ent- 
schieden ablehnt mit den Worten: „Ich glaube 
nicht, daß ein richtiger Physiker jemals auf solche 
seltsamen Gedankengänge verfallen würde.“ Ding- 
lers scharfe Kritik: „Welch ein katastrophales 
Programm für eine Physik“ ist also an die falsche 
Adresse gerichtet. Im 2. Aufsatz ergibt sich be- 
züglich der Frage der Determiniertheit eine ge- 
wisse Übereinstimmung zwischen Dingler und 
Planck, da beide «lie Indeterminiertheit des 


menschlichen Willens 
„zegenwärtiger Wille“ 


proklamieren. soweit er 
(Planck) bzw. „aktiver 
Wille“ (Dingler) ist. In der Tat beruht auch (das 
„eindeutig-methodische System“  Dinglers. das 
(dieser ausführlich in seiner „Methode der Physik“. 
dargestellt hat, auf der Tätigkeit dieses aktiven 
Willens: vgl. die Besprechung in dieser Zeitschrift 
Bd. 19 (1939). S. 59. Immerhin findet er schließlich 
„bei Planck zuletzt doch wieder eine bedauerliche 
Unsicherheit und Verwirrung im Methodischen. 
welche eindringlich auf die Notwendigkeit hin- 
weist. hier zunächst einmal Ordnung zu schaffen“, 
Solange freilich der von Dingler gewünschte. auf 
Grund seines e.-m. Systems durchzuführende 
lückenlose, allen Forschungsergebnissen Rechnung 
tragende Aufbau der Physik noch nicht vorliegt. 
wird das physikalische Weltbild wohl auch weiter- 
hin durch jene von Dingler verworfenen, bzw. nur 
als zeitweiliger Notbehelf aufgefaßten Theorien be- 
stimmt werden, ohne daß übrigens deswegen 
alle Physiker auch mit «der erkenntnistheoretischen 
Stellunz Plancks einverstanden zu sein brauchten. 


Berlin. E. Mosch. 274 


J. P. AREND, Die wirkliche Mechanik 
des Geschehens im Lebenden und 
Toten. 456 8. Zürich und Leipzig 1941, Rascher 
Verlag. Preis geb. 10 M. 


Der Verf. glaubt. das „kausalitätserschließende 
bBezugssystem” gefunden und „das höchste Ziel 
der Wissenschaft, die Erkenntnis des Wirklichen, 
das Wahre“ erreicht zu haben. Das Buch setzt 
sich mit der gegenwärtigen Wissenschaft in allen 
ihren Zweigen auseinander. Ein kurzer Bericht, 
wie er hier nur möglich wäre. ist wegen der neu- 
artigen Vorstellungen und Gedankengänge des 
Verf. leider nicht möglich: der Ref. muß sich daher 
auch jeden Urteils über die Beweiskraft der Aus- 
führungen sowie über Bedeutung und Tragfähig- 
keit der Ergebnisse enthalten. 
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Dr. ©. KLOUCEK, Das Prinzip der fort- 
seleiteten Verformung als Weg zur 
Ausschaltung der Unbekannten aus 
dem Formänderungsverfahren. 391 S. 
Berlin 1941, Verlag Wilhelm Ernst & Sohn. Preis 
brosch. 18 M. 


Der Verfasser beschäftigt sich in diesem Buche 
mit dem Spannungsnachweis in Rahmen, deren 
Knoten durch geeignete Randbedingungen am 
Tragwerk unverschieblich sind oder als unverschieb- 
lich angenommen werden. Der Verschiebungszu- 
stand ist dann allein durch die Knotendrehwinkel 
bestimmt. Diese werden aus den Bedingungen des 
Gleichgewichts der Stabendmomente am Knoten 
berechnet, die sich als Funktion der Stabbelastung 
und der Knotendrehwinkel anschreiben lassen. Die 
linearen Gleichungen können in diesem Falle leicht 
dureh Iteration gelöst werden. Dasselbe gilt vom 
Gaußschen Algorithmus, wenn der Ansatz die drei- 
gliedrige Form erhält. In diesem Falle bestehen 
zwischen den unbekannten Knotendrehwinkeln ein- 
fache Kennbeziehungen, mit denen die Lösung 
schnell für jedes einzelne Belastungsglied ange- 
veben wird, um daraus das endgültige Ergebnis 
abzuleiten. Diese Rechenvorschrift läßt sich oft in 
allgemeinen Fällen als Näherung verwenden. Sie 
ist bei der Bedeutung der Rahmenberechnung für 
die Baustatik in zahlreichen Abwandlungen vorge- 
schlagen und in der Praxis des Bauingenieurs je 
nach Geschmack und Erziehung verwendet worden. 
Hierzu zählt auch die Rechenvorschrift des Ver- 
fassers. Dabei werden zunächst die Knotendreh- 
winkel aus dem Gleichgewicht der Stabend- 
momente eines Knotens bei voller Einspannung 
der angeschlossenen Stäbe bestimmt und die Fort- 
pflanzung dieser Verformung auf die benachbarten 
Knoten verfolgt. Die endgültigen Knotendreh- 
winkel ergeben sich durch Superposition und damit 
die Stabendmomente, welche die Gleichgewichts- 
bedingungen erfüllen. Ich vermag die Erwartung 
(les Verfassers, die dieser an die Brauchbarkeit der 
Lösung stellt, nicht zu teilen. Die übliche Rechen- 
vorschrift ist übersichtlich und die Auflösung der 
linearen Gleichungen gerade hier kein Problem. 
Die theoretischen Betrachtungen des Verfassers 
sind ausführlich erläutert und durch Tabellen und 
zahlreiche Beispiele zur statischen Untersuchung 
von Baukonstruktionen ergänzt. 


Dresden. K. Beyer. 330 


Prof. Dr. AUGUST THUM und Dr.-Ing. OTTO 
PETRI, Steifigkeit und Verformung 
von Kastenquerschnitten mit Ver- 
rippungen, Wanddurcehbrechungen u. 
bei Behinderung der Querschnitts- 
verwölbung (VDI-Forschungsheft 409, Beilage 
zu „Forschung auf dem Gebiete des Ingenieur- 
wesens“, Ausgabe B, Bd. 12, Juli/August 1941). 
28 S. m. 72 Bildern und 8 Zahlentaf. Berlin 1941. 
VDI-Verlag G.m.b.H. Preis 5 M. 


Ausgehend von dem Gedanken, daß der unge- 
heuere, Verbrauch an Werkstoffen in der heutigen 
Zeit zu äußerster Sparsamkeit zwingt, weisen die 
Verfasser vor allem auf die weit bedeutsamere 
Forderung hin, daß nur der Leichtbau auf vielen 
Gebieten des Maschinenwesens den höheren Ge- 
schwindigkeiten und Leistungen, die als Folge 
einer natürlichen Entwicklung anzusehen sind, ge- 
recht werden kann. Für den Leichtbau genügt 
jedoch die Güte des Werkstoffes allein noch nicht. 
Das Wesen dieser Bauweise muß vielmehr darin 
erblickt werden, Konstruktionen mit größtmög- 
licher Steifigkeit zu schaffen, was nur durch gün- 
stige Formgebung und Kenntnis der Bean- 
spruchungsart erreicht werden kann, um unnötige 
Werkstoffanhäufungen zu vermeiden. 


Nach einer mathematischen Entwicklung des Be- 
eriffes Steitigkeit bei Zug-, Druck-, Biege- und Ver- 
drehungsbeanspruchung werden bei gleichem Ma- 
terialverbrauch die Starrheiten von 5 der gebräuch- 
lichsten Hohlquerschnitte und einem Querschnitt 
offener Bauart miteinander verglichen. Für die sich 
stets schwierig gestaltende rechnerische Ermittlung 
der Schubspannungen bei Verdrehung von Hohl- 
querschnitten wird ein Näherungsverfahren ange- 
zeben. Die Ergebnisse der sich anschließenden 
Biege- und Verdrehungsversuche, für die ein Probe- 
körper mit kastenförmigem Hohlquerschnitt als 
Grundform, einmal ohne und dann mit 6mal verschie- 
denartiger Innenverrippung, gewählt worden ist. 
erweisen recht befriedigend die Brauchbarkeit des 
rechnerischen Verfahrens. Um keine Lücke zu 
lassen. wurden die Versuche auch auf die Ermitt- 
lung der Spannungen bei den unvermeidlichen 
Wanddurchbrechungen «les Hohlkörpers ausgedehnt 
und darüber hinaus selbst noch der Einfluß (der 
Behinderung «der  Querschnittsverwölbungen bei 
Verdrehung untersucht. 

Bei Besprechung dieser vorzüglichen Forschungs- 
arbeit, die zu klaren Erkenntnissen führte, darf 
nicht übersehen werden, daß das Ziel nur durch 
eine wohlüberlegte Versuchsanordnung und durch 
mühsame, sorgfältige Feinmessungen erreicht wor- 
den ist. 
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Prof. Dr. Dr.-Ing. AUG. FÖPPL ($) und Dr. 
LUDWIG FÖPPL, o. Prof. a.d. Techn. Hochschule 
München, Drang und Zwang.eine höhere 
Festigkeitslehre für Ingenieure. 
Erster Band, 3. Aufl. XII + 346 S. m. 70 Abb. 
München u. Berlin 1941, Verlag R. Oldenbourg. 
Preis geb. 15,70 M. 

Die vorliegende dritte Auflage des bekannten 
Buches weist gegenüber der 1924 erschienenen 
zweiten nur geringe Änderungen auf. Durch- 
gehend ist kurz auf neuere, den betreffenden Ge- 
genstand behandelnde Arbeiten hingewiesen. An 
größeren Änderungen ist zu erwähnen, daß dem 
Abschnitt über Plattenbiegung ein ausführlicher 
Bericht über die Arbeit von Lewe: „Die Lösung 
des Pilzdeckenproblemes durch Fouriersche Rei- 
hen“ angefügt ist und daß in dem letzten, die 
Scheiben behandelnden Abschnitt die Paragraphen 
fortgelassen sind, die sich mit Näherungslösungen 
beschäftigten. Als Grund für diese Kürzung gibt 
Verf. an, daß nach seiner Ansicht Näherungs- 
lösungen von Scheibenproblemen nicht mehr die 
gleiche praktische Bedeutung zukomme wie früher, 
nachdem man heute derartige Probleme mit den 
Methoden der Spannungsoptik behandeln könne. 

Wahrscheinlich infolge der wenigen Änderun- 
gen ist es möglich geworden, den Preis dieser 
Neuauflage trotz guter Ausstattung im Vergleich 
zu den heute sonst für derartige Werke gleichen 
Umfanges üblichen Preisen außerordentlich nie- 
drig zu halten. Das ist sehr zu begrüßen und wird 
zur Verbreitung des Buches beitragen. Für die 
Neuauflage wäre aber doch eine Überarbeitung 
des Ganzen sehr wünschenswert gewesen. Dabei 
hätten die späteren Zusätze wirklich eingearbeitet 
werden müssen und vor allem wäre manches, 
z. B. bei der Behandlung und Anwendung der bei- 
den Minimalsätze, schärfer zu fassen gewesen. 
Auch sähe man gern weitere neue Arbeiten be- 
rücksichtigt, wie etwa die in dieser Zeitschrift er- 
schienenen von C. Weber über die Eingrenzung 
elastischer Größen mittels der beiden Minimal- 
prinzipe. 

Doch auch in der vorliegenden Form wird das 
Buch theoretisch interessierten Ingenieuren noch 
gute Dienste leisten können. 

Dresden. 
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Dr. HEINRICH BLASIUS, Studienrat an der In- 
eenieur-Schule zu Hamburg, Wärmelehre, 
Physikalische Grundlagen vom ttech- 
nisehen Standpunkt. 3. neubearb. Aufl. 
VIII + 299 S. m. 126 Fig., 67 Aufg. u. 15 Tabellen. 
Hamburg 1941, Verlag von Boysen & Maasch. 
Preis geb. 8 M. 

Fast alle neueren Bücher über Technische 
Wärmelehre knüpfen eng an die von Zeuner be- 
eründete und von Mollier fortgeführte Darstel- 
ung dieses Gebietes an. Es ist ein unbestreit- 
bares Verdienst des Verfassers, daß er sich be- 
müht, die Grundlagen der Wärmelehre dem An- 
fänger auch auf anderen Wegen zugänglich zu 
machen. Die Form wurde in erster Linie durch 
die Bedürfnisse von Ingenieurschulen gegeben. 
Ausgehend von eingehend  durchgerechneten 
Zahlenbeispielen werden die einzelnen Gebiete 
dann unter sparsamster Verwendung mathema- 
tischer Hilfsmittel allgemein behandelt. Trotz 
dieser elementaren Darstellungsweise, die den 
mathematisch wenig vorgebildeten Leser mit gro- 
Bem Geschiek an das Wesen der Vorgänge heran- 
führt, wird überall wissenschaftliche Genauigkeit 
angestrebt. Das Buch geht inhaltlich erheblich 
über die klassische Thermodynamik der Gase und 
Dämpfe hinaus und bringt im gleicher Darstel- 
lungsart Abschnitte über feuchte Luft, Misch- 
dampf, Lösungen usw. Freilich vermißt man ge- 
rade hier das Fehlen von Schaubildern, die dem 
Anfänger auch dort einen klaren Überblick geben, 
wo er dem verwickelten Rechnungsvorgang nicht 
folgen kann. 

Bedauerlich ist, daß der Verfasser auch in der 
neuen Auflage von genormten oder festliegenden 
Bezeichnungen abweicht, so daß der Gebrauch des 
Buches neben den üblichen Taschenbüchern min- 
lestens sehr erschwert wird. Besonders unglück- 
lich ist die Bezeichnung „Wärmestrom“ statt En- 
thalpie oder Wärmeinhalt. Dazu werden Ver- 
wechslungen noch dadurch erleichtert, daß die 
„Innere Energie“ mit „Energieinhalt“ bezeichnet 
wird. Eine allgemeine Einführung der üblichen 
Bezeichnungen könnte den Wert des besonders 
für Ingenieurschulen sehr geeigneten Buches 
sicher noch erhöhen. 
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OTTO SCHAUB, Beiträge zur Strömungs- 
lehre. 47 8. Biel (Schweiz) 1940, Selbstverlag 
0. Schaub-Scheurmann. Preis brosch. 3 fres. 


Der Verfasser versucht, für die Helmholtz- 
Kirchhoffschen Strömungen mit freier Strahlgrenze 
eine Lanze zu brechen, weil nach seiner Meinung 
die darauf aufbauende Widerstandstheorie durch 
die 40 Jahre jüngere Prandtische Grenzschichts- 
theorie zu Unrecht verdrängt würde. Die Strö- 
mungen mit und ohne freie Strahlgrenzen unter- 
scheidet er als „offene“ und „geschlossene“ Strö- 
mungen, die als stationäre Strömungen gleich- 
berechtigte Lösungen der hydrodynamischen 
Gleichungen darstellen. Wieweit man sie als 
stationäre Strömungen beobachten kann, hängt 
natürlich davon ab, welchen Einfluß die vernach- 
lässigte Reibung besitzt und ob Stabilität der 
Lösung gegenüber kleinen Strömungen vorhanden 
ist. Im Gegensatz zu allen eingehenden theoreti- 
schen und experimentellen Untersuchungen glaubt 
der Verfasser, daß die „offene“ Strömung, also 
diejenige mit freien Strahlgrenzen, sich unmittel- 
bar durch den Beschleunigungsvorgang ausbildet. 
Er beweist dies für den zweidimensionalen Fall 
eines Spaltes in ebener Platte, bei dem er eine 
gleichmäßige Geschwindigkeitsverteilung im Spalt 
oder Strahl statt der ungleichmäßigen Geschwin- 
digkeitsverteilung im Spalt bei der geschlossenen 
Potentialströmung erwartet, mit folgenden Worten: 


„Die gegebenen Anfangsbedingungen sprechen 
nun keineswegs zugunsten der Entstehung einer 
geschlossenen Strömung. Denn der im Augenblick 
des ersten Bewegungsbeginnes über die ganze 
Durchflußöffnung gleichförmig verteilte Druck- 
sprung ruft!) einer ebensolchen Beschleunigung 
der Massen links der Spaltöffnung, womit eine 
unveränderliche Geschwindigkeit eingeleitet wir. 
Man braucht sich nur den Vorgang so vorzustellen, 
daß eine die Öffnung abschließende Membrane 
plötzlich entfernt wird, um zu erkennen, daß jede 
andere Druckverteilung rechts der Membrane gegen 
das hydrostatische Grundgesetz der Druckfort- 
pflanzung verstoßen würde. 

Und weil im Augenblick des Bewegungsbeginnes 
die Zentrifugalkräfte, die bei geschlossener Strö- 
mung die Ursache der unterschiedlichen Druck- 
verteilung im Öffnungsquerschnitt sind, fehlen, bzw. 
verschwindend klein sind im Verhältnis zum an- 
treibenden Druckgefälle, so ist überhaupt gar 
keine Möglichkeit zur Ingangbringung einer un- 
gleichförmigen Strömung gegeben.“ 

Um die Tatsache, daß jeder Stromfaden wegen 
der Kontinuität an allen Orten gleichzeitig be- 
schleunigt werden muß, wobei die einzelnen Strom- 
fäden bei ihrer verschiedenen Geometrie eine 
unterschiedliche Aufteilung des Druckgefälles ver- 
langen, und daß die ungleichen Druckgradienten 
auch eine ungleichmäßige Geschwindigkeitsver- 
teilung erzeugen. kümmert sich der Verfasser 
offenbar nicht. Denselben Mangel an Beweiskraft 
weisen auch seine übrigen Ausführungen, bei denen 
er z. B. das gelegentliche Auftreten geschlossener 
Strömungen aus anfänglich offenen Strömungen 
etwa im Venturirohr umständlich erklären muß, an 
den entscheidenden Stellen auf. Man stellt daher 
mit einer gewissen Befriedigung fest, daß der- 
artige Abhandlungen offenbar nur noch im Selbst- 
verlag erscheinen können. 


Braunschweig. 


1) Der kleine Fehler im Satzbau steht im Original, ist 
aber natürlich belanglos. 
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Dr. CLEMENS SCHAEFER, o. ö. Prof. d. Physik 
an der Univ. Breslau, Einführung in die 
Maxwellsche Theorie der Elektrizi- 
tät und des Magnetismus. 4. verb. Aufl. 
142 S. m. 33 Textfig. Berlin 1941, Verlag Walter 


de Gruyter & Co. Preis geb. 6 M. 

Verf. behandelt von dem Coulombschen Gesetz 
ausgehend zunächst die Gesetze der Elektrostatik 
und Magnetostatik. In den folgenden Kapiteln 
werden die Gesetzmäßigkeiten über das Magnet- 
feld des elektrischen Stromes und die Faradaysche 
Induktion erörtert, wobei Verf. auch Anwendun- 
gen stationärer und quasistationärer Felder be- 
rücksichtigt. Das Schlußkapitel befaßt sich mit 
der Ausbreitung von elektrischen Wellen in Iso- 
latoren und Leitern. Dabei wird auch der Gültig- 
keitsbereich der Maxwellschen Theorie im Zu- 
sammenhang mit der neueren Quantentheorie in 
den Kreis der Betrachtungen gezogen. 

Es wäre zu begrüßen, wenn Verf. in der näch- 
sten Auflage auch das Pohl-Miesche-Maßsystem 
der Elektrodynamik heranziehen würde. 

Die Darstellung ist gegenüber früheren Auf- 
lagen verbessert und klarer gestaltet. Ich emp- 
fehle das ausgezeichnete Buch allen Physikstudie- 
renden deshalb besonders, weil es in didaktischer 
Hinsicht kaum zu übertreffen ist. 
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Prof. G. MAHLER, Physikalische For- 
melsammlung. Neubearbeitet von Prof. 
K. Mahler, Studienrat a. d. Schubartoberschule 


‚Aalen in Württemberg. 7. verb. Aufl. (Sammlung 
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Göschen Bd. 136). 152 S. m. 69 Fig. Berlin 1941, 
Verlag Walter de Gruyter & Co. Preis geb. 1,62 M. 

Die bekannte und vielbenutzte Formelsammlung, 
die alle wichtigeren Sätze und Formeln der Phy- 
sik enthält, soweit ihr Verständnis nicht die Hilfs- 
mittel der höheren Mathematik voraussetzt, ist in 
einer neuen Auflage erschienen. Die Änderungen 
gegen die vorhergehende beziehen sich im wesent- 
lichen auf die Anpassung der Bezeichnungen an 
die Vorschriften des Normenausschusses. Sach- 
liche Erweiterungen haben insbesondere die über 
die Fluglehre und die Spektroskopie handelnden 
Abschnitte erfahren, womit den Fortselfritten der 
Wissenschaft Rechnung getragen ist. 
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Dr. KARL KOMMERELL, em. o. Prof. d. Math. 
a. d. Universität Tübingen, Vorlesungen über 
analytische Geometrie der Ebene. 
VII + 368 S. m. 92 Fig. Leipzig 1941, K. F. Koeh- 
ler Verlag. Preis 20 M. 

Seine „Vorlesungen über analytische Geometrie 
ıles Raumes“ t) vervollständigt der Verfasser durch 
‚ten vorliegenden Band zu einem geschlossenen 
I,ehrbuch: er hat in ihm zugleich mit der Darbie- 
tung des notwendigen Stoffes als Aufgabe die Her- 
anführung des Anfängers an die Gedanken, die in 
ler Geometrie und in ihrer analytischen Behand- 
lung wirksam werden. Von den elementaren Metho- 
dien der Parallel- und Polarkoordinaten ausgehen 
werden die Grundbegriffe der metrischen Geometrie 
entwickelt und sehr bald — unter Benutzung der 
komplexen Zahlen — die metrischen Transforma- 
tionsgruppen eingeführt. Diese geben Anlaß zur 


Aufstellung der projektiven Gruppe und mit ihr 
zusammen zur frühzeitigen Berücksichtigung des 
Erlanger Programmes von F. Klein. Nunmehr kann 
der Verfasser den Aufbau des geometrischen Lehr- 


gebäudes gleichzeitig von der geometrischen und von 
der analytischen Seite her gestalten. Zunächst leitet 
das Polarsystem des Kreises zur Einführung der 
uneigentlichen Elemente und die Zentralprojektion 
zur Einführung der homogenen Koordinaten. Die 
Entwickelung der Grundbegriffe der projektiven 
Geometrie im Verein mit der Einführung der Linien- 
koordinaten ermöglicht die Feststellung des Duali- 
tätsprinzips, die projektive Erzeugung der Kurven 
zweiter Ordnung und zweiter Klasse, sowie die 
Ableitung der Sätze von Pascal und Brianchon. Die 
Einordnung der Kreisgeometrie durch die unendlich- 
fernen imaginären Kreispunkte und die Formel von 
Laguerre führen zur Cayley-Kleinschen Maßbestim- 
mung, in der der erste Abschnitt seinen Abschluß 
findet. 

Der zweite Abschnitt ist den Eigenschaften der 
Kegelschnitte gewidmet, die in homogenen Parallel- 
koordinaten zum Ausdruck gelangen. Eine sorg- 
fältige Untersuchung der gemeinsamen Punkte 
zweier Kurven zweiter Ordnung führt zum Kegel- 
schnittbüschel und zur Kegelschnittschar, der die 
projektive Fassung der Brennpunktseigenschaften 
angeschlossen wird. Der dritte Abschnitt führt die 
Dreieckskoordinaten ein und wendet sie an auf die 
projektive Theorie der Kegelschnitte, sowie der 
Kollineationen und Korrelationen, dabei die schon 
vorher angebahnte Benutzung der Invariantentheorie 
wesentlich vertiefend. Ein Anhang entwickelt die 
notwendigen Sätze über zwei- und dreigliedrige 
Determinanten. 

Dem Buch sind dieselben Vorzüge nachzurühmen 
wie den Vorlesungen über analytische Geometrie 
des Raumes: wie diese entspricht es in voilsten 
Maße der Absicht des Verfassers, „der Geometrie 
zu dienen“. 
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1) Vgl. die Besprechung in Z. Math. Mech. 


angew. 
Bd 20 (1940) S. 294. 


Dr.-Ing. W. MEYER ZUR CAPELLEN, Mathe- 
matische Instrumente. (Mathematik für 
Physiker und Ingenieure, herausgeg. von Prof. Dr. 
A. Kratzer, Bd. 1.) X + 247 S. m. 210 Fig. Leipzig 
1941, Akademische Verlagsgesellschaft Becker & 
Erler Kom.-Ges. Preis geb. 30 M. 


Die Anwendung mathematischer Methoden dringt 
in immer weitere Gebiete vor. Dabei handelt es 
sich nicht nur um mathematische Formulierungen 
und daran anschließende theoretische Untersuchun- 
gen, sondern oft um zahlenmäßige Auswertung der 
sich ergebenden Formeln. Da für diese oft mathe- 
matische Instrumente benutzt werden können, be- 
steht heute allgemein ein großes Interesse für diese 
Apparate. Die Veröffentlichungen darüber sind 
aber sehr zerstreut und teilweise schwer zugäng- 
lich. Es ist daher zu begrüßen, daß Verf. diese zu- 
sammenfassende Übersicht gegeben hat. 


Er faßt sein Thema im weitesten Sinne und 
rechnet zu den mathematischen Instrumenten nicht 
nur, wie man es heute meist tut, Instrumente zum 
Integrieren, Differenzieren, Rektifizieren, Analysie- 
ren, Synthetisieren, Umzeichnen, Auftragen und 
Kurvenzeichnen, sondern behandelt außerdem 
Rechenschieber (in der Hauptsache „System Darn- 
stadt“) und sonstige logarithmische Rechengeräte. 
Rechenmaschinen und Vorrichtungen zur Lösung 
von Gleichungen. Dabei werden nicht nur mecha- 
nisch, sondern auch elektrisch oder elektrooptisch 
arbeitende Apparate besprochen. 


Die Darstellung ist getriebetechnisch orientiert. 
Als Einführung gibt Verf. daher einige Abschnitte 
über die verschiedenen Rechengetriebe, die ja 
heute in vielen Maschinen — ich erwähne nur die 
Kommandogeräte der Flugabwehr — verwendet 
werden. Besonders hervorgehoben zu werden ver- 
dient der Abschnitt über Rechenmaschinen, der 
sich durch zahlreiche instruktive Lichtbilder aus- 
zeichnet, die den inneren Aufbau und die Wir- 
kungsweise (Schaltwerk, Zehnerübertragung, Lösch- 
vorrichtung, Schlittentransport usw.) der verschie- 
denen Systeme erläutern. Diese Bilder stammen 
fast alle von Herrn F. Hecht (Zella-Mehlis) und 
auch Aufbau und Text dieses Abschnittes lehnen 
sich eng an ein Manuskript des Herrn Hecht 
an. Während bei den Rechenmaschinen die Ent- 
wicklung zu einem gewissen Abschluß  ge- 
kommen zu sein scheint, kann bei den Ge- 
räten zum Lösen von Gleichungen nur von einem 
Anfang gesprochen werden. Die eigentlichen 
mathematischen Instrumente dagegen, insbesondere 
Planimeter, harmonische Analysatoren und Integra- 
toren stehen mitten in einer lebhaften Entwicklung. 
Zahlreiche Neukonstruktionen sind in den letzten 
Jahren geschaffen, an vielen wird noch gearbeitet. 
Vor allem konzentriert sich heute das Interesse auf 
Apparate zur Integration von Differentialgleichun- 
gen, auf die Verf. vielleicht etwas ausführlicher 
hätte eingehen sollen. 


Natürlich behandelt Verf. allgemein nur eine Aus- 
wahl dieser Instrumente, nämlich solche, die er für 
wichtig oder entwicklungsfähig hält. Um ein weiter- 
führendes Studium einzelner Fragen zu erleichtern, 
gibt er am Schluß ein rund 270 Nummern um- 
fassendes Verzeichnis des neueren Schrifttums. 


Das vorliegende Buch ist der erste Band einer 
neuen Reihe von Lehrbüchern, die unter dem ge- 
meinsamen Titel „Mathematik für Physiker uni 
Ingenieure“ erscheinen sollen. Mir scheint mit 
dieser Darstellung der mathematischen Instru- 
mente ein guter Anfang gemacht zu sein. Man 
kann nur wünschen, daß die folgenden Bände, vor 
allem was Ausstattung und Druck anlangt, ebenso 
vortrefflich ausfallen, wie der vorliegende, 
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64 Nachrichten 


Z. angew. Math. Mech. 
Bd. 22 Nr.1 Febr. 1942 


Dr. LOTHAR SCHRUTKA, o. Prof. d. Math. a. d. 
Techn. Hochschule Wien, Leitfaden der In- 
terpolation. VII+80 S. Wien 1941, Springer- 
Verlag. Preis steif geh. 4,80 M. 

Dieser ganz elementar geschriebene Leitfaden 
behandelt in den ersten 4 Abschnitten (Steigungen. 
Steigungen mit gleichen Argumenten, Interpolation 
von Polynomen und parabolische Interpolation als 
Näherungsverfahren) die Interpolation für beliebige 
Argumente. In dem nächsten Abschnitt über para- 
bolische Quadratur, in dem in der Hauptsache, die 
Sehnentrapez-, die Tangententrapez- und die Simp- 
sonsche Regel mit Fehlerabschätzung abgeleitet 
werden, geht Verf. dann zu gleichabständigen Ar- 
gumenten über. Ausführlich wird in einem weiteren 
Abschnitt die Interpolation bei solchen gleich- 
abständigen Argumenten behandelt. Der letzte 
Abschnitt endlich beschäftigt sich mit der Mouton- 
schen Aufgabe, also der Untertafelung. An Stoff 
wird etwas weniger gebracht als das, was man 
über diesen Gegenstand in den gebräuchlichen 
Lehrbüchern der praktischen Analysis findet, abge- 
sehen vielleicht von dem letzten Abschnitt, der 
sich in diesen Büchern oft nicht findet. Besonders 
sei auf die vom Verf. gebrachten Restabschätzun- 
gen verwiesen. Die Darstellung ist sehr geschickt, 
so daß das Buch als erste Emführung durchaus 
empfohlen werden kann. 
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WALTHER LIETZMANN, Lustiges und 
Merkwürdiges von Zahlen u. Formen. 
5. durchges. u. erg. Aufl. VI+ 296 S. m. 197 Fig. 
u. 16 Taf. Breslau 1941, Verlag Ferdinand Hirt. 


Preis geb. 8 M. 

Das jetzt in Neuauflage erschienene unterhaltende 
Buch bringt, wie die früheren Auflagen, eine große 
Fülle anregenden und interessanten Stoffes, der aus 


Anekdoten, Gedichten, Stücken aus Romanen uni 
Essays, Spielen usw. in bunter Folge zusammen- 


getragen ist, so daß man sich gern müßige Stun- 
den damit vertreibt. Von den früheren Auflagen 
unterscheidet sich die Neuauflage hauptsächlich 
wohl dadurch, daß ein Säuberungsprozeß an 
Material und Autorennachweisen durchgeführt 
wurde und die ausgemerzten Bestandteile einwand- 
freien Platz gemacht haben. 
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Ferner sind bei der Schriftleitung folgende Bücher 
eingegangen (ausführliche Besprechung bleibt vor- 
behalten): » 


Dr. habil. ALBERT KOCHENDÖRFER, Dozent 
für Physik a. d. Techn. Hochschule Stuttgart, Pla- 
stische Eigenschaften von Kristal- 
len und metallischen Werkstoffen. 
(Reine und angewandte Metallkunde in Einzeldar- 
stellungen, herausgegeben von W. Köster, Nr. 7.) 
XI + 312 S. m. 91 Abb. Berlin 1941, Verlag Julius 
Springer. Preis geb. 28,50 M. 


Dr. ROBERT HAUSSNER, oö. Prof. a. d. Univ. 
Jena, Sammlung von Aufgaben und 
Beispielen zur analytischen Geome- 
trie der Ebene mit den vollständigen 
Lösungen. 2. unveränderte Aufl. (Sammlung 
Göschen, Bd. 256). 139 S. m. 22 Fig. im Text. 
Berlin 1941, Verlag Walter de Gruyter & Co. 
Preis geb. 1,62 M. 


Prof. Dr. A. WITTING, Differentialrech- 
nung. 3. neubearb. Aufl. (Sammlung Göschen, 
Bd. 87). 201 S. m. 95 Fig. und 200 Beispielen. 
Berlin 1941, Verlag Walter de Gruyter & Co. Preis 
geb. 1,62 M. 


Dr. phil., Dr.-Ing. e. h., Dr. phys. e. h. R. W. 
POHL, o. ö. Prof. a. d. Univ. Göttingen, Ein- 
führungindie Mechanik, Akustik und 
Wärmelehre X-+322 S. m. 527 Abb. Berlin 
1941, Verlag J. Springer. Preis geb. 15,80 M. 


NACHRICHTEN 


Vortragsreihe. 
Neuere Festigkeitsprobleme des Ingenieurs. 


Im Außeninstitut der Technischen Hochschule 
Berlin veranstaltete dieses in Gemeinschaft mit 
dem Verband Deutscher Elektrotechniker, Bezirk 
Berlin und der Berliner Mathematischen Gesell- 
schaft im Herbst 1941 eine Vortragsreihe: Neuere 
Festigkeitsprobleme des Ingenieurs. Die Reihe be- 
stand aus 8 Vorträgen. Es sprachen: 

1. und 2. Dr.-Ing. habil. K. Marguerre, 
Adlershof: Einführung in die Elastomechanik. Ein- 
führung der Grundbegriffe, der Balken unter Bie- 
gung und Torsion, Formänderungsarbeit, die ener- 
getische Methode in der Balkentheorie, Knickung 
des Stabes. 

3. Prof. Dr.-Ing. G. Mesmer, Darmstadt: Ex- 
perimentelle Verfahren zur Spannungsbestimmung. 
Messung der Spannungsänderung an der Oberfläche, 
Messung der absoluten Spannung an der Öber- 
‚fläche, Messung an Modellen. 

4. Prof. Dr.-Ing. K. Klotter, Berlin: Schwin- 
gungserscheinungen im Bau- und Maschinenwesen. 
Unterschied zwischen der statischen und dynami- 
schen Betrachtungsweise, die Grundbegriffe der 
Schwingungslehre, erzwungene Schwingungen, Rolle 
der Dämpfung. 

5. Prof. Dr.-Ing. R. Grammel, Stuttgart: Prak- 
tische Verfahren zur Lösung technischer Eigenwert- 


probleme. Eigenwertprobleme mit endlich vielen 
Freiheitsgraden, Eigenwertprobleme von kontinuier- 
lichen Systemen, Methoden von Rayleigh, Ritz und 
Galerkin. Besonders dankenswert, daß der Vor- 
tragende auch sein eigenes Verfahren darstellte. 

6. Dr.-Ing. K. Marguerre, Adlershof: Knick- 
und Beulvorgänge. Stabilitätsprobleme, Energie- 
methode, Plattenknickung. Geschlossene Lösungen 
bei der Rechteckplatte. 

7. Dr.-Ing. habil. W. Flügge, Adlershof: Festig- 
keit von Schalen. Drehschalen, Kesselboden, un- 
symmetrische Belastung, Drehschalen mit nega- 
tiver Krümmung, windschiefe Vierecke, Behälter- 
theorie, allgemeine Biegetheorie der Drehschalen. 

8. Prof. Dr. G. Hamel, Berlin: Ausblick auf 
neuere Ansätze. Nichtlineare Beziehungen zwischen 
Formänderung und Spannung. Prinzip der virtu- 
ellen Arbeiten, Minimalprinzipe, Zusammenhang mit 
den entsprechenden Sätzen der Allgemeinen Me- 
ehanik. 

Berlin. 
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Persönliches. 


Anläßlich seines 65. Geburtstages wurde dem 
o. Professor für Strömungslehre an der Techn. 
Hochschule Berlin Dr.-Ing. H. Föttinger die 
Goethe-Medaille für Kunst und Wissenschaft ver- 
liehen. 
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